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Relaxationsmethode* 
Von Miroslay Kotal, Prag 
Mit 9 Textabbildungen 


Zusammenfassung. Die Relaxationsmethode ist eine Abart der Iterationsmethoden zur Lésung 
des Gleichungssystems. Sie wird als eine Netzmethode fiir Berechnung der elektromagnetischen 
Felder dargestellt. In diesem Falle mu Laplacesche oder Poissonsche Differentialgleichung in 
eine Differenzengleichung transformiert werden. Singulare Punkte in dem Potentialfeld verur- 
sachen bei den Berechnungen groBe Schwierigkeiten. Als Beispiel wird eine Lésung des Stirn- 
feldes eines Generators eingefihrt. 


J. Einleitung 


Kine Untersuchung der elektromagnetischen und elektrotechnischen Felder ist 
wegen der verschiedenen Randbedingungen sehr schwierig. Bei experimentaler 
Untersuchung dieser Felder werden entweder direkte Methoden oder passende Modelle 
benutzt. Von den graphischen Methoden ist die Lehmannsche Methode die bekannteste. 
Manchmal kommen auch die analytischen Methoden in Betracht. Fiir die unkompli- 
zierten Randbedingungen wird die Methode der konformen Abbildung (Schwarzsches 
Theorem) angewandt. Die Netzmethode, die als ein spezieller Fall einer Relaxations- 
methode dargestellt werden kann, ist wenig bekannt. Diese Methode, die von 
Southwell! erarbeitet wurde, ist eine Abart der GauB-Seidelschen Iterationsmethode 
fiir die Lésung der simultanen Gleichungen. Eine ahnliche Methode von Liebmann 
beruht auch auf Iterationsverfahren. Diese Methoden werden auch zur Lésung der 
Differentialgleichungen, die mit gewisser Annaherung als Differenzengleichungen 
ausgedriickt werden kénnen, beniitzt. In unserem Falle fiihren wir nur das Rechen- 
verfahren der Lésung einer partiellen Differentialgleichung eines Potentialfeldes ein, 
die in allgemeiner Form lautet: 


7) Oy a) oy 
6x (4 se.) ae ay (B =) i f(x, y). (1) 
Wenn die Koeffizienten A, B als eine Funktion von x, y dargestellt werden kénnen, 
erhalt man eine Gleichung mit veranderlichen Koeffizienten. Falls A = B= 1 ist, 
hat man: @ @2 
gat + oye =, 9). (2) 
Diese Gleichung mit f(x, y) auf der rechten Seite wird auch Poissonsche Gleichung 
genannt. Wenn wir f(x,y) auf Null bringen, bekommen wir eine Laplacesche 
Gleichung: ep . a 
ax? dy? 0 (3) 
oder: Vig 0: 


* Diese Benennung ist meistens in der Literatur benutzt, obwohl sie sich nur auf Berechnung 


eines Gleichungssystems bezieht. 
1 R. V. Southwell: Relaxation Methods in Theoretical Physics. Oxford University Press. 


1949. 
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Man nennt die Methode Relaxationsmethode, weil das ganze Rechenverfahren 
als Berechnung der Fachtriger betrachtet werden kann. In diesem Falle halten wir 
die Reste, die sog. Residuen, fiir ungewiinschte auBerliche Krafte, die Schritt fir 
Schritt auf erforderliche Kleinheit herabgedriickt werden (eine Relaxation). Bei der 
Losung der elektromagnetischen Felder? bedeuten die Residuen die magnetischen 
Quellen der betreffenden Punkte des Netzes (das heiBt die rechten Seiten der 
Poissonschen Gleichungen). Alle Berechnungen werden unter einer Voraussetzung 
eines ebenen Feldes durchgefiihrt. Das dritte Glied tritt also in der Gl. (3) nicht auf. 

Das gleiche Rechenverfahren kann auch bei biharmonischer Gleichung angewendet 


werden: ey dtp at 
bat + aan age + oye = He, 9). (4) 


Il. Uher die Wahl des Netzes und des entsprechenden Bereiches 


Man iiberdeckt das Feld, in dem die Gl. (3) gelost werden soll, mit einem quadra- 
tischen Netz von der Maschenweite a. Manchmal wird auch das Feld in Drei- oder 
Sechsecke geteilt. In speziellen Fallen kénnen auch die Parallelogramme beniitzt 
werden. In einzelnen Netzpunkten werden die numerischen Werte mit vorgeschrie- 
bener Genauigkeit, die der Gl. (2) und (3) gentigen, berechnet. Den Bereich wahlt 
man auf solche Weise, dai sein Rand durch Netzpunkte gebildet werden soll, oder 
daB die Randpunkte sich immer zwischen inneren und auferen Netzpunkten befinden. 


III. Rechnerische Voraussetzungen 


Ks werden im Felde, das der Gl. (3) geniigt, ein beliebiger Punkt 0 und vier Nachbar- 
punkte 1, 2, 3, 4 angegeben (Abb. 1). Mit y,; bezeichnen wir ein zugehériges Potential 
des betreffenden Punktes 7. Wir gehen vom Taylorschen Satz der Reihenentwicklung 
aus. Man erhalt fiir das Potential im Punkt 1 


4 
Oe a pee BE 
ga Wag se vo=(1 + Gy eager Om ae (5) 
Voraussetzung fiir die Giltigkeit des Taylorschen Satzes ist, daB 
alle vorkommenden Ableitungen existieren und stetig sind. Ahn- 
licherweise erhalt man auch fiir andere Netzpunkte: 


Ota F) 2 92 \ 
ita oy ere As 
Ps =e yo=(1 ty dy x 2! dy? ++} Pos (6) 
e 7) 
ere oie deals 7) Qs? on? 
Abb. 1. Va tine er (1 =aGaspa oh sayungs alee (7) 
Relaxationsfigur a 
Ca eae 2 72 
FAs oy bi 0 Oni 
Pa ere é Po (1 Ay dy | Oy dy? Sora Sie eral .) Po- (8) 


Nach Vernachlassigung der Glieder von dritter und héherer Ordnung und nach 


der Durchfiihrung der Elimination von £. und a aus den Gl. (5) bis (8) bekommen 
wir fir die Poissonsche Gleichung . 


2 (L— Pe | fe Vo.) 2 Y2— Po , Pa Po 
i Qs I = 


a, + as a dg + Oy Gs azn 
OD oy 
ma ane Bye ane f(a, Y)o- (9) 


* A. Weigand: Die angeniherte Berechnung ebener und rotationssymmetrischer Potential- 
felder mit Hilfe des Differenzenverfahrens. Berlin. 1953. 
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Wenn wir 
Cy = 45-4) 
=, =D 
einsetzen, wird die Gleichung des Potentials im Punkte 0 von der Form sein: 


(Pi + 3) 6? + (Yo + Ya) a2 — f(x, Y)y a? B? 


Po iv 2 (a + b?) (10) 
Im Falle des quadratischen Netzes ist a gleich 6b 
2 
oe en a ae Y)o (11) 
und fiir die Laplacesche Gleichung?: 
++ 
= Vite ie ep: : Ps + Ya ; (12) 


Da die Glieder von zweiter und héherer Ordnung in Gl. (5), (6), (7), (8) vernach- 
lassigt werden, erweist die Endgleichung (10), (11), (12) einen gewissen Fehler, der 


Abb. 2. Diagonalrelaxationsfigur Abb. 3. Dreiecksnetz Abb. 4. Sechsecksnetz 


annahernd der Ordnung a? ist. Nach einer ausftihrlichen Errechnung*: > bekommen 
wir aus der Gl. (5) bis (8), wenn gleichzeitig a, =a, =a; =a,=a eingesetzt wird: 


4 


Si —ane= 2/2 (28 + OH) 4 (Oe 4 i) 4 lM ag 
daraus folgt: : 
worin die Restglieder aN 

Pe tet (Bi A) as) 


sind. Wenn wir in Nachbarpunkten die Ableitungen von dritter Ordnung und im 
Mittelpunkt die Ableitungen von vierter Ordnung in das Rechenverfahren einfihren, 
erhalten wir folgenden Ausdruck fiir die Restglieder: 


4 2 
Fee Ol 1), (16) 
; A o4 o* 
worin M, absoluter Maximalbetrag der Ableitungen von vierter Ordnung na und ae 


ist. 


8 Hugh McKibbin: Relaxation Methods I., Il. Electr. Rev. 167, 87 (1955). 
4 ZL, Collatz: Numerische Behandlung von Differentialgleichungen. Berlin: Springer-Ver- 


lag. 1955. 
aa L. V. Kantorovitsch-V. I. Krylov: Die Methoden hodherer Analysis. Moskau. 1951 
(in russischer Sprache). 
7s 
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Entsprechend stellen wir auch die Taylor-Entwicklung fiir die Relaxationsfigur 
nach Abb. 2: 


4 
Dy Pi 4% =P $92 + %s + Hs — 4% = 


to 
Nes PY an a (a2 O*H “ee 
4\a1 da heer dak ot O oat ay? | dy? as 
a Aq, - 68 os 
4 (a B+ Soe + 1B te + ee) + of. (17) 


Man kann beide Ausdriicke (13), (17) unter gewissen Voraussetzungen zusammen- 
fassen, was zur Verdichtung der Netzpunkte und dadurch zur gréBeren Genauigkeit 
der Berechnung fihrt. Fiir Dreiecksnetz ergibt sich (Abb. 3): 


3a? 
3'9.—690= 21 Ago + mT Ste es 
Se 
33 po BH ee ea 27 Ayo : 
+-(43 ox ioe Oat Oy? + 15 16 dx? dy? + 76 me) 4 atieiels (18) 
daraus 
6 
3 
DV: — 6 = FV Ago + Ro (19) 
= 
worin 


9 4 
ewe ege ny ee (20) 


ist. In erster N&iherung kann man also die Laplacesche 
Gleichung durch die lineare Beziehung ersetzen: 


| 

| 

by! gg = a oe a Pa (21) 
Abb. 5. Relaxationsfigur Und die Poissonsche Gleichung: 

zur Erhdhung der Be- 


Soe 2 
rechnungsgenauigkeit % = Pi + P2 + Ys ¢ Pat ste f(x, ue aay (22) 


Man kann auch ein Sechsecksnetz verwenden (Abb. 4): 


3 
So 3D) = 7 VAG + Ry, (23) 
worin 
a@ (3 9 dy f 
Poe ae jah a oe ae (24) 
ist. Fir praktische Anwendung nimmt diese Gleichung folgende Gestalt an: 
1 
Po = 3 (Pi + P2 + Ys); (25) 
1 2 
Po= x (M1 + G2 +%s) — Mx, Wo- (26) 


Weitere Erhéhung der Berechnungsgenauigkeit® erhilt man aus der Relaxations- 
figur Abb. 5. Aus der Gl. (5) bis (8) folgt: 


o 1 o* 

a+ — 29 = 0(58) 12 a(S) to, 
‘oa 1 o4 

oto. — 290 = 0 (52) 12 as( Pe) 4 


° H. Motz-L. Worthy: Calculation of the Magnetic Field in Dynamo-Electric Machines 
by Southwell’s Relaxation Method. J. Inst. Electr. Engr. 1945, 522. 


(27) 
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Ahnlicherweise kénnen wir schreiben: 


Ce 4 ot 
pr + oy — 24 = 40 (52) 9) at (te) ae) 
4 


O Po HDq Ce) 
Px + Py — 2g = 4a?( | 3 at et) oe 
Unter Voraussetzung, daB Glieder von héherer Ordnung als a! vernachlassigt werden, 
a4 4 
und nach Elimination ss e 0 erhalten wir: 
x oy 


4 1 
BV? Po = 3 (V1 + G2 + Ps + Ys — 490) — Gq (GY +2 +437 +P —490)- (29) 


Daraus kénnen wir den Wert gy, ermitteln: 


4 1 
Po = 5 (V1 + Ya + Ys + Pa) — GO (GL + Px + Gs + Hy’). (30) 


IV. Matrizenrechnung und Relaxationsmethode 


Der Funktionswert g kann als eine lineare Kombination der inneren Werte 9; 
und der Randwerte ,, dargestellt werden’: ®: 
a Tis 9; +2 Vin Pe + Hi = 9, (31) 
; ; 
worin gy; absolutes Glied der Gleichung ist. 


T,; und U,,; sind die numerischen Koeffizienten. Diese sind von Netzdimensionen 
und der Anordnung der Netzpunkte abhangig. Sie erweisen aber eine Abhangigkeit 
weder von den Randwerten noch von den berechneten Werten. Nach dem bekannten 
Maxwell-Satz iiber reziproke Verschiebung gilt: 7;; = 7';. Diese Gleichheit, die 
immer fiir eine harmonische Funktion erfillt wird, ist notwendige, aber nicht hin- 
reichende Bedingung fiir Konvergenz des Relaxationsverfahrens. Randwerte 9, 
k6nnen als stabil betrachtet werden. Gl. (31) driickt eine lineare Transformation aus. 

Es sei gy; die Lésung der partiellen Differentialgleichung (3) im Punkte 7. Dann 
kann die Funktion ¢;” nach »-ten Iterationen als gesuchte Losung und einem ge- 
wissen Fehler y,, sog. Residuum, ausgedriickt werden: 


9 =O +4; y=1,2...n, (32) 
Pe = Ve + Ve = Pr pimple Dine i (33) 
Dieser Fehler y; verschwindet beim Randwert , = 0. 


Zur Vereinfachung dieses Verfahrens setzen wir die Randwerte », gleich Null. 
Ahnliche Voraussetzungen macht Poincaré bei der ,,Ausfegemethode“ (,,balayage“). 
Fiir den betrachteten Fall kann man schreiben: 


n 
ah Tin PO + Vi = pi. (34) 

k=1 
Der Grundgedanke des Relaxationsverfahrens besteht darin, durch systematische 
Korrekturen Ag,®, Av”... 4g, der Versuchswerte die Residuen_schrittweise 


immer kleiner zu machen. Diese Korrekturen ergeben die neuen Residuen: 
n 
YO =Ol) +S Te Ape-?: (35) 
k=1 


7 BE. Stiefel: Uber einige Methoden der Relaxationsrechnung. Z. angew. Math. Physik 8, 
1 (1952). 

8 G. Temple: Relaxation Methods Appl. to Linear Systems. Proc. Roy. Soc. London, Ser. A 
169, 476. 
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Die Niaherungswerte ¢,”) und Residuen yp,” fassen wir zu Vektoren zusammen. 
Als Skalarprodukt dieser zwei Vektoren bekommen wir folgenden Ausdruck : 


n n n nN 
0 yw) = S79 yO =D GOD Tin Gi +L PO vie (36) 
(=a iad) eee fad 
Indem wir 
ly tg 
Qe) => = Tn VP OX ae Pp ; (37) 
(i, Via 


setzen, das heiBt eine quadratische Funktion, worin >’ 7, py, eine positiv 
(ik) 
definite quadratische Form ist. Durch Differenzieren dieses Ausdruckes folgt: 


4) f (v) (38) 
09, 7 Tin Pe” + Ve 
was eigentlich i-tes Residuum ist: 
; 7) : 
yp = Par (39) 
Weiterhin kénnen wir schreiben: 
y; = grad Q, (40) 


das bedeutet, daB die 7-te Komponente von y gleich der partiellen Ableitung von Q 
nach , ist. Dadurch ist eine Auflésung der gegebenen Gleichung vollstandig aqui- 
valent mit der Bestimmung des Minimums der quadratischen Funktion. In der 
Mechanik stellt diese Funktion eine Potentialenergie des Systems dar. 

Diese Auffassung fiihrt zur geometrischen Bedeutung der Relaxation. Zu jedem 
Naherungsvektor y;”) wird ein Gewichtsvektor p; angeordnet wie folgt: 

Ag™=Ap®  p, * 9, (41) 

worin A ein Parameter ist. Weitere Behandlung kann nach Ritz durchgefiihrt werden. 

Bei einer Blockrelaxation® ist A der negativ genommenen Summe der Residuen, 
die durch Anzahl der gekreuzten Strecken des Randgebietes dividiert wird, gleich. 
Die Anderung der Residuen in mehreren Punkten verhalt sich auf ahnliche Weise 
wie die Residuen in einzelnen Punkten. 

Bei einer Methode des stirksten Abstieges wird Gewichtsvektor p, entgegengesetzt 


gleich dem Residuenvektor y,; gewahlt. Aus Matrizenanalysis folgt, daB in diesem 
Falle Korrekturfaktor 4 der reziproke Rayleighsche Quotient des Residuenvektors ist: 


Rope (42) 
a WV; 


V. Potentialtheorie 


Kine harmonische Funktion ist eine beliebige, zweimal stetig differenzierbare 


Losung der Gleichung © Bip | sigh aay 
| 
da? y3 Gat ae (43) 
bzw. a ap 
Ou a oye v. (3) 


Bei der Losung dieser Gleichung kénnen folgende Randwertaufgaben betrachtet 
werden!?: 


®° H. McKibbin-Webb: Relaxation Methods I., II. Electr. Rev. 1953, 67, 190. 


*0 Ph. Frank-R. Mises: Die Differential- und Integralgleichungen der Mechanik und Physik, 
Bd. I. 1930. 


Relaxationsmethode 99 
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a) Dirichletsches Problem (erste Randwertaufgabe): gegeben y am ganzen Umfang 
eines geschlossenen Bereiches, gesucht » im Innern oder im AuBenraum des Bereiches. 


b) Neumannsches Problem (zweite Randwertaufgabe): gegeben 


oy 
on’ 


das heiBt die Normalkomponente am Umfang des Bereiches, gesucht @. Das Rechen- 
verfahren entspricht in diesem Falle aber nicht dem des Dirichletschen Problems. 


c) Problem der Warmeleitung. 


d) In der Elektrotechnik kann das vierte Randwertproblem betrachtet werden. 
Dieses wird als Refraktion oder Zwischenebene charakterisiert. 


VI. Randbedingungen fiir die Relaxationsfiguren 


a) Lineare Interpolation 


Unter einer Voraussetzung, da die Funktion am Rande als linear betrachtet 
werden kann, setzen wir ein: 


a = A, (44) 

Ag = Pa, (45) 

a3 = 7%, (46) 

4 = 17 4. (47) 
Aus Gl. (9) folgt: 


_ fi | Pe | Ps | Pa 
l= g 


1 1 
(pt+r)p ' t+q)q ' (ptnrr vol, =7|=0.- (48) 


Fiir Poissonsche Gleichung ist an der rechten Seite anstatt Null folgender Ausdruck: 


Cnty 
co 


Wenn p, g oder r wechselweise gleich Eins gelegt werden, konnen die gesuchten Aus- 
driicke fiir die Figuren mit den im Verhialtnis 7, p, q verkiirzten Relaxationsstrahlen 
hergestellt werden. 

b) Parabolische Interpolation 


Fiir den einzigen Strahl einer Relaxationsfigur kann eine parabolische Interpolation 
verwendet werden, indem durch die Punkte 1, 2, 3,4 eine Parabel, die folgender 
Gleichung gentigt, gelegt wird (Abb. 6): 


tg} Te On@ 
ames ee des) | peek (49) 


Nach Bestimmung der Differenzen 


0 
: a od, sf ae ne 
oe? 
oF = G+ 91-202 (51) 


bekommt man fiir den Punkt 3: 


] 
G2 = G2 +P (Gs — 91) + yPP(Me + 91 — 292): (52) 
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i | “A : pa 3 \4 
oebehae: seNae 2 a 
Abb. 6. Parabolische Interpolation Abb. 7. Hilfsbild bei einer Singularitat 


VII. Rechenverfahren bei der Singularitat 


Wenn die gesuchte Funktion in einem Punkte, z. B. im Punkte B (Abb. 7) unbe- 
stimmten Wert erreicht, kann’ die obige regelmaBige Relaxationsrechnung nicht 
benutzt werden!!. 

Wir fiihren zuerst die Laplacesche Gleichung (3) in die Polarkoordinaten tier; 
Punkt B wird als Mittelpunkt betrachtet: 


eh 0 Or 
a (53) 
Nach Einsetzung 
y (r, 3) = RB (r) 6 (8) (54) 
bekommt man: 
oe 
CR aR av} 
eho 75 pr rote E(r) GH = 9. (55) 
Wenn wir 
do 1 
t= — capt BiB): pee 
schreiben, hat man: 
CR dR 
aa, te Capac: Le: (57) 
Aus der Erlosung der Gleichung 
6” + nv20=0 (56a) 
folgt: 
6 = C, (r) cosn & + Cy (r) sinn B. (58) 
Ks gilt = = 0 fir @=0 und = «. Daraus bekommt man 
= Circo 2.0, (59) 
Weiterhin 
ka 
Me Ne a (be ==G, zs 3 Se) 24s. 2) (60) 
und schlieBlich erhilt man: 
k=-+0o 
p= 2 A,r cosnd. (61) 
Auf einem Kreise, der mit kleinem Radius @ um den Punkt B gezogen wird, gilt: 
op 
poe” (2) =o. (62) 


4 H. Motz: The Treatment of Singularities of Partial Differential Equations by Relaxation 
Methods. Quart. Appl. Math. 4, 371 (1947). 
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Kine entsprechende Reihenentwicklung liefert: 
A~=Q"Ay, (b=1,2...). (63) 


Wenn g@ zu Null strebt, verschwinden die negativen Exponenten und gm kann 
durch folgenden Ausdruck dargestellt werden: 


ze ee, Lox: 
g=A,+A,r% cos — 9 +A,r * cos “* 9 + Azar * cos 4... (64) 
Falls « = 22, hat man (Abb. 8): 
1 3 
os dp — 
yg=A,+A,r? cos-7 + A,r cos? + Agr? Coss +... (65) 


Abb. 8. Quadratisches Netz mit Singularitét im Punkte B 


Man legt einen Einheitskreis durch die Punkte 1, 2, 3, 4. Die Werte der Konstanten 
der Reihe (65) sind wie folgt: 


Ay = 0°25 (9; + 2 + Ps + Pa), (66) 
H 135° 1 45° 
A, = (2 — Ps) 9g COS 5 (Ys — Pi) 3% COS Be? 
A, = 0°191 (pz — ys) — 0°462 (yp, — Ps); (67) 
ahnlicherweise : 
A, = 0°354 (1 — 2 — Y3 + Pa); (68) 
A, = — 0°191 (p, — G4) + 0°462 (p3 — 2). (69) 


Die Werte in den Punkten 1’, 2’, 3’, 4’ werden durch folgende Gleichungen be- 
stimmt: 
Py = 0°593 py + 0°236 gp, + Oly; + 0°073 gy, 
Py = 0236 y, + 0°457 9, + 0-21 q3 +0:1 os 
YQ =O19, +0219, + .0:457 p, + 0°236 gy, 
Gy = 0073.9, + Oleg, +0236 y, + 0°593 gy. 


In anderen Netzpunkten berechnet man die numerischen Werte nach den Aus- 
driicken (12), (48). Residuen in den speziellen Punkten werden nur addiert oder 
subtrahiert. 


| 
(70) 


VIII. Beispiel der Lésung 


Als Beispiel wurde das Stirnfeld eines dreiphasigen Generators gewahlt (Abb. 9). 
Es handelt sich um den Fall, in welchem nur der Stator mit Strom durchflossen 
wird. Wir betrachten einen Augenblick, wenn der Momentanwert in einer Phase sein 
Maximum erreicht und in anderen Phasen halbe Stréme flieBen. Die zugehdrige 
MMK im Stirnfelde werden der Wicklungsanordnung und der Verbreitung der MMK 
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Abb. 9. Stirnfeld eines dreiphasigen Generators 


in den Nuten entnommen. Beim Herunterdriicken der Residuen wurde auch die zwei- 
malige Blockrelaxation benutzt. 


1X. Zusammenfassung 


Aus der oben angefiihrten Analysis ist sichtbar, daB die Relaxationsmethode eine 
geniigend prazise mathematische Methode zur Losung von Differentialgleichungen ist. 
Diese Methode wurde in unserem Falle zur Loésung einer harmonischen Gleichung 
angewandt. Es ist ihr Vorzug, da zum Unterschied von anderen mathematischen 
Methoden keine prazisen Rechenmaschinen erforderlich sind und die Rechen- 
operationen mit Rechenschieber durchgefiihrt werden kénnen. 


(Eingegangen am 20. Dezember 1956) 
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Uber die richtige Farbenwiedergabe 
Von F. Lachner, Wien 
Mit 5 Textabbildungen 


Zusammenfassung. In der vorliegenden Arbeit wird versucht, vom physiologisch-physikali- 
schen Standpunkt aus zu zeigen, da den iiblichen Farbenreproduktionsverfahren Grenzen gesetzt 
sind und da insbesondere eine richtige Wiedergabe des reinen Spektrums nicht méglich ist. 
Es wird gezeigt, daB eine Verbesserung durch verzweigte Verfahren méglich ist, so daB wenigstens 
nicht zu satte Farben richtig wiedergegeben werden, was bei additiver Wiedergabe leichter 
gelingt als bei subtraktiver, bei der noch eine zweite Matrix zu realisieren ist und auBerdem ein 
trichromatisches Wiedergabewei8 nétig wire. Angedeutet wird ferner die Méglichkeit der An- 
wendung von mehr als drei Wiedergabegrundfarben, die auch eine hinreichende Wiedergabe des 
reinen Spektrums erméglichen wiirde. 


Mit den im Handel befindlichen Erzeugnissen werden oft zufriedenstellende farbige 
Bilder erzielt. Trotz sorgfaltigster Durchfiihrung der farbenphotographischen Ver- 
fahren gibt es aber in gewissen schwierigeren Fallen véllige Versager, insbesondere 
bei der Farbaufnahme eines Spektrums. Die vorliegende Abhandlung befaft sich mit 
den Ursachen dieser Versager und mit den naturbedingten Grenzen der Leistungs- 
fahigkeit der tiblichen Dreifarbenverfahren. Es wird ferner angegeben, durch welche 
Abanderungen eine Verbesserung mdglich ist. Die Abhandlung soll das Problem 
der richtigen Farbenwiedergabe vom Standpunkt der physiologisch-physikalischen 
Farbenlehre behandeln und in den Kreisen der reinen Techniker, bei denen viel irrige 
Meinungen verbreitet sind, aufklirend wirken. 

Fiir die Beschaftigung mit dem Problem der naturgetreu erscheinenden Farben- 
wiedergabe geniigen nicht technische und physikalische Kenntnisse allein; es sind 
vielmehr auch solche aus dem Gebiet der physiologischen Optik erforderlich. (Analog 
muB z. B. der technische Spezialist fiir Elektroakustik auch Kenntnisse aus dem 
Gebiet der physiologischen Akustik haben.) Vor der Inangriffnahme des Problems 
ist es daher nétig, etwas iiber die beim technischen Unterricht kaum angedeutete 
physiologische Farbenlehre anzugeben. Wie in jeder Naturwissenschaft, bilden Versuche 
die Grundlage. Vorerst die Versuche iiber die Mischung farbiger Lichter, die sog. 
additive Farbmischung. Man fand als Versuchsergebnis, dali es zwischen vier 
farbigen Lichtern stets eine Mischungsbeziehung gibt (sofern nicht der wenig wahr- 
scheinliche Sonderfall einer Mischungsbeziehung innerhalb einer Gruppe von drei 
der vier gegebenen Farben vorliegt). Jede Mischungsbeziehung kann durch eine 
Farbgleichung dargestellt werden, die sich wie eine lineare Vektorgleichung verhilt. 
In der Farbgleichung wird jede Farbe durch ein Symbol, zweckmaBig in gleicher Weise 
wie ein Vektor, bezeichnet; also mit einem fett gedruckten oder tiberquerten Buch- 
staben oder durch Anwendung von Frakturbuchstaben!. 


1 Die in der Tensorrechnung vielfach wbliche, falschlich oft symbolfrei genannte Bezeich- 


nungsweise, bei der ein Vektor F in gleicher Weise wie seine Koordinate F, nur durch einen 
Buchstabenindex bezeichnet wird, hat sich in der Farbenlehre als wenig zweckmaBig erwiesen. 


Man hat oft mit Farbtripeln zu tun, wobei man den k-ten Vektor des Tripels mit /,, und seine 
i-te Koordinate mit F',; bezeichnet. Doppelindizes kennzeichnen aber in der erwaéhnten Tensor- 
schreibweise nur die Konstanten einer Transformationsmatrix. Es hat sich aber oft als zweck- 
maBig erwiesen, fiir die Matrix ein Symbol ohne Indizes zu beniitzen. Zum Unterschied von 
einem Vektor wahlt man oft groBe griechische Buchstaben. Sehr zweckmafig und auch leicht 
bei der Schreibmaschine sowie bei der Setzmaschine anwendbar hat sich die Bezeichnungs- 
weise nach F. Jung erwiesen, der Vektoren durch einfache Uberquerung und Dyaden bzw. Matrizen 
durch Doppeliiberquerung bezeichnet (Triaden durch dreifache usw.). Fir die Setzmaschine wiire 
auch die folgende Schreibweise méglich: Vektoren durch fette (oder halbfette) kleine und Matrizen 
durch fette groBe Buchstaben. In der vorliegenden Abhandlung wird vorwiegend die Jungsche 


Bezeichnungsweise beniitzt. 
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In einer Farbgleichung steht vor jedem Farbsymbol im allgemeinen noch ein 
Zahlenfaktor (Skalar). Beim Mischungsergebnis zweier gleicher Farblichter ist es 
die Zahl 2. Der Faktor wichst linear mit der Energiezunahme des Farblichtes. Nach 
dem vorhin iiber die Ergebnisse der Farblichtermischungen Angegebenen hat man 
also allgemein eine implizite Farbgleichung 


Op. =O (1) 


wobei iiber den Index k von 1 bis 4 summiert wird. (Sofern im Einzelfall nicht aus- 
driicklich das Gegenteil angegeben, ist ein Buchstabenindex, der in eiem Produkt 
zweimal vorkommt, Summationsindex.) In (1) sind die Konstan- 
ten a, zum Teil negativ. 


Die Glieder von (1) mit negativen Konstanten werden nun 
auf die rechte Gleichungsseite gebracht (wobei sie positiv werden). 
Es ergeben sich nun zwei Méglichkeiten: 1. Auf einer Gleichungs- 
seite steht ein Glied und auf der anderen drei Glieder. In diesem 
Fall kann die eine Farbe durch entsprechend dosierte Mischung 
der drei anderen nachgebildet werden. 2. Je zwei Glieder stehen 
auf jeder Gleichungsseite, was bedeutet, daB eine Mischung zweier 
Farben einer Mischung der anderen beiden Farben gleich aussehend 
gemacht werden kann. In diesem Falle kann also keine der vier 

6(2,-%,1%) Farben einer Mischung der drei anderen gleich aussehend gemacht 
Abb.1.Grund-und Werden. Die weniger bekannte Tatsache, daB auch der Fall 2 auf- 
AufriB eines Farb-  treten kann, ist eine der Ursachen dafiir, daB die Dreifarbenreproduk- 

rake baec tion nicht in allen Fallen eine richtige Wiedergabe erméglicht. Prak- 

tisch brauchbar kann also nur der Fall 1 sein, und zwar nur dann, 
wenn die wiederzugebende Farbe allein auf einer Gleichungsseite steht und in der 
Gleichung kein negatives Vorzeichen vorkommt. 


Man wihlt nun drei bei den weiteren Untersuchungen unverandert belassene sog. 
Bezugsfarben B; zwischen denen es keine Mischungsbeziehung gibt. Es sei nun F 


irgendeine zu messende Farbe. Dann kann die (1) entsprechende Mischungsbeziehung 
zwischen diesen vier Farben auch in der Form 


B= FYB, (2) 
geschrieben werden, wobei der Summationsindex 7 von | bis 3 lauft. Die (positiven 


oder negativen) Konstanten F;” sind die ,,FarbmaBzahlen“ (im B-System), wobei 
im Falle des Normsystems der obere Index weggelassen wird. Man hat dann 


P= Py Bok Be BeBe (2) 
Diese Gleichung kann man als Vektorgleichung auffassen und mit Hilfe der dar- 


stellenden Geometrie graphisch darstellen, indem man den B; ein im allgemeinen 
schiefwinkeliges Vektordreikant zuordnet. Abb. 1 zeigt ein Beispiel, bei dem ein 


rechtwinkeliges Normbezugstripel N; gewaihlt wurde. Auf diese Weise kann jeder 
Farbe ein Farbvektor in einem ,,Farbraum‘ zugeordnet werden. Die F, nennt man 
auch Farbkoordinaten oder Farbwerte. 


Wird die Helligkeit (Leuchtdichte in Kerzen pro Quadratmeter) einer Farbe, 
nicht aber die ,,Farbart‘‘ geindert, dann Andert sich nur die Lange, nicht aber auch 
die Richtung des Farbvektors. Damit wird auch die Bedeutung eines Zahlenfaktors 
vor einem Farbvektor verstindlich. 

Es werden nun die einzelnen Farben eines beziiglich der Wellenlangenskala igo- 
energetischen Normspektrums gemessen, dessen hellste Stelle (bei 555 nm) die Leucht- 
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dichte 1 K/m2 hat. Alle seine Wellen gemischt ergeben das Normweils Hy. Als Bezugs- 
farben bzw. Eichfarben wahlt man zweckmaig die Spektrallichter Hg 436, Hg 546, 
Cd 644 nm, die unmittelbar gemischt, Normweif ergeben wiirden. Als Ergebnis er- 


halt man die S;? bzw. S (A). Die Spitzen dieser Farbvektoren bilden die Spektral- 
raumkurve. Sie liegt auf einer (allgemeinen) Kegelfliche, die durch die Spektral- 
vektoren gebildet wird. Hat man nun ein Farblicht mit der spektralen Strahldichte- 
verteilung (Knergieverteilung) F (2), so ist es nun einleuchtend, daB dazu der Farb- 


vektor = * 
F=O|FSds (3) 


gehort, wobei C eine Konstante ist. Da der Kegel iiberall nach innen konkav ist, 
mu8 der Farbvektor jeder reellen Farbe innerhalb dieses Kegels liegen. 

La8t man nun ein Farblicht auf drei Photozellen mit verschiedener spektraler 
Empfindlichkeitsverteilung P,, (A) fallen, so erhalt man drei elektrische Spannungen 
(bzw. Stréme), die unter gewissen Voraussetzungen Farbkoordinaten darstellen kénnen, 
aber stets positiv sind. Das Gebiet der reellen Farben, also der gesamte Spektralkegel, 
mu innerhalb des durch die entsprechenden Bezugsvektoren A, gebildeten Dreikants 
liegen. Eine einfache Uberlegung ergibt, daB 


P,Ap= 8 = 8/7 B; (4) 
sein mu, wenn man durch die Photozellen wirkliche Farbraumkoordinaten gewinnen 


will. Haben die Photozellen (bzw. die mit den entsprechenden MeBfiltern versehenen 
Photozellen) von dieser ,, Luther-Bedingung“< 


(nach Prof. Luther) abweichende Kurven, so sind sie fiir eine solche Farbmessung 
unbrauchbar. Gleich aussehende Farben mit verschiedener spektraler Energie- 
verteilung, sog. bedingt-gleiche Farben, hatten verschiedene Farbmafzahlen. 

Beim Farbensinn sind nun tatsachlich drei Sehstoffe (und ein weiterer fiir das 
Dammerungssehen) vorhanden, die sich wie die drei Photozellen verhalten und die 
spektralen Grundfarbenempfindlichkeiten G; (A) haben, wobei_ selbstverstandlich 
G; = 8,° ist. Dazu gehéren die physiologischen Grundfarbenvektoren G;. Auch die 
drei Aufnahmeschichten der Dreifarbenreproduktion verhalten sich wie die drei 
Photozellen. Der Photozellenspannung entspricht die Exposition bzw. die Transparenz 
des positiven Dreifarbenauszugs, nimlich 7'p, (sofern mit dem richtigen Kontrast 
kopiert wurde). Die spektrale Empfindlichkeitsverteilung der (mit den Dreifarb- 
auszugs-Aufnahmefiltern versehenen) Aufnahmeschichten mu8 also ebenfalls S et) 
sein, das heiBt es miissen die Aufnahmefilter die Luther-Bedingung erfiillen. 


Die Farbaufnahme muB also stets eine einwandfreie Farbmessung darstellen, 
was durch die Erfiillung der Luther-Bedingung erméglicht wird. Man hat 
8,0 A, =8; Ny. (6) 


Die Normbezugsvektoren N, liegen beim Internationalen CIE-System im Gebiet 
der virtuellen Farben und die S;-Kurven (Abb. 2) sind tabellarisch gegeben. Da die 
eine S,-Kurve ein Nebenmaximum hat und sich ein Nebenmaximum bei Filterkurven 


schwer realisieren 148t, macht man das Tripel der A, verschieden von dem Tripel 
der W;.: Die P, = S,4-Kurven sind dann. lineare Kombinationen der S;-Kurven, 
zweckmiBig nach Abb. 3. Die Ordinaten sind das Produkt aus Schichtempfindlich- 
keit und Aufnahmefiltertransparenz. Die Konstanten der linearen Kombination 


(bzw. der Luther-Bedingung) bilden eine Matrix a, die kontragredient ist zur Matrix A, 
die aus den Koordinaten der A,, gebildet wird. 
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Bei der farbenphotographischen Aufnahme der Farbe F hat man F,4 = Tr 


und bei der Wiedergabe F;” = 7'p;. Man muB die reellen Wiedergabegrundfarben 
(die meist einem Drittel des Spektrums entsprechen) durch die den Wiedergabe- 
Positivtransparenzen entsprechenden Dosierungen 7'p;’ steuern. Man hat also 


P2PporAa =F es, (7) 
Das Tripel der A, ist aus den vorhin 
angegebenen Griinden virtuell und das 


CIE-NORM 
(1BK) 
S} SA) 
S.A) 


5 7102 mp 4 £ 6 72.07 nm 
$3=% Sahih $422 (nm) 
Abb. 2. Normspektralwertkurven Abb. 3. Aufnahmegrundfarben (spektrale 


Empfindlichkeitswertverteilung) 


der B, reell. Man hat daher eine Transformation von dem einen Bezugssystem in das 
andere. Fir die Koordinatentransformation hat man die Matrix V anzuwenden. Hs ist 


ge = Vie PF Gs (8) 
wobei 


V=(4 B4), (9) 


ist, nach bekannten Regeln der Tensorrechnung. Fiir die Matrizen bzw. Dyaden 
wurde hier zum Unterschied von Vektoren die doppelte Uberquerung gewahlt. (Auch 
groBe griechische Buchstaben werden beniitzt.) Die Gl. (8) k6nnte man auch in der 
Form 


F? — y FA (8’) 
schreiben. Es darf daher 7'p;’ nicht mit 7'p;, tibereinstimmen, wie dies beim gewodhn- 
lichen unkorrigierten additiven Verfahren der Fall ist. 

Bei dem durch ein sog. ,,verzweigtes Verfahren“ korrigierten additiven Verfahren 
sind die korrigierten Positivtransparenzen 7'p;’ lineare Kombinationen (mit teilweise 
negativen Konstanten) der unkorrigierten Positivtransparenzen 7'p,. Beim analogen 
korrigierten Farbfernseher sind die Steuerspannungen der Wiedergabegrund- 
farben lineare Kombinationen der drei Photozellenspannungen. Es kénnen alle 
Farben richtig wiedergegeben werden, die sich innerhalb des Dreikants der Wieder- 
gabegrundfarben befinden, die man daher méglichst satt wihlen wird. Man wird also 
moglichst reine Spektrallichter 436, 546 (oder besser Tl 535), 644 nm wihlen. 
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Hatte man fiir das subtraktive Verfahren ideale, nur ein Drittel des Spektrums 
absorbierende Zweizonenfarben Gelb, Purpur, Eisblau (Cyanblau), die in den anderen 
beiden Dritteln des Spektrums keine Nebenabsorptionen haben, dann ware es dem 


additiven Verfahren aquivalent. Die blauviolette Wiedergabegrundfarbe B, sollte 
nur durch den Gelbfarbstoff allein gesteuert werden, was aber nicht der Fall ist. 
Wegen der Nebenabsorption des Purpurfarbstoffes im Blaudrittel trigt derselbe 
auch zu einer zusitzlichen Blaulichtsteuerung bei. Es ist nun naheliegend, die Gelb- 
konzentration proportional der Purpurkonzentration durch ein sog. ,,Maskenverfahren “ 
herabzusetzen. Analog bei der Kompensation anderer Nebenabsorptionen. 

Bei Anwendung eines trichromatischen Wei kann das subtraktive Verfahren 
exakt auf ein aquivalentes additives zuriickgefiihrt werden. Man geht dabei aus von 
den fiir das additive Verfahren korrigierten drei Einzeldiapositiven, deren Transparenzen 
die T'p;’ sind. Thre Schwarzungen sind die S,’ = — lg T'p,’, die zu dem Hilfsvektor 


S’ = —lg 7’ pa To) 


zusammengefaBt werden koénnen. (Die Definition der Funktion eines Vektors ist 
hier zu beachten!) Man stellt nun ein neues Tripel von Diapositiven her, dessen 
Schwarzungsvektor 


S’=MS' : (11) 


ist. Dazu kann das sog. Maskenverfahren nach E. Albert beniitzt werden. In jedem 
dieser Diapositive wird das Silber durch je einen Farbstoff ersetzt, so daB z. B. die 
Flachenkonzentration des Gelbfarbstoffes c, = H, 8S,” ist. Analog fiir Purpur und 
Kisblau. Man hat also 


¢=HS". (12) 


Nur die Glieder in der Hauptdiagonale der Matrix von H sind von Null verschieden. 

Der k-te Farbstoff hat bei c, — 1 fiir das 7-te Grundfarbenlicht die Dichte D,,;. 
Die aufeinander liegenden drei Farbstoffschichten haben fiir das i-te Farblicht die 
Gesamtdichte 


Dy = Dui Cn = Deir €x- (13) 
Der Gesamtdichtevektor ist ie 
D=D,¢. (14) 
Fiir 
D=S' (15) 


wird das subtraktive Verfahren dem additiven aquivalent. 
Durch Zusammenfassung von (11) bis (15) erhalt man 


Ks mu also 223 
D,HM =I, (16) 


also gleich dem Identitatsaffinor (Idemfaktor) sein. Daraus ergibt sich 
M = (0, = 40,284 17) 


womit man also die Nebenabsorptionen bzw. Nebendichtenkonstanten (D,,, fiir + + /) 
kompensieren kann. Beim subtraktiven Verfahren hat man also die 
Matrix V und die Matrix M zu realisieren! (Auf diesen Umstand wurde vom 
Verfasser erstmalig klar hingewiesen.) 
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Es wird nun untersucht, ob V und I nicht durch eine einzige Matrix ersetzt werden 
konnen. Ks gilt ¥ 
Fa=AoF (18) 
und % 
F=4A,F* (19) 
und analog fiir B. Man hat also 


js Al ape sy 1 —B, antilog (—D) = B, antilog (—D,H M-— log (V- A, F)). 


In der Klammer hat man 


D, H Mog (V A,“ F). 


Ware nicht log zwischen M und J, so konnte man als Maskierungsmatrix die Produkt- 
matrix (M V) realisieren und wiirde sich die getrennte Anwendung eines verzweigten 


Verfahrens zur V-Realisierung ersparen. Oder man kénnte sich die Maskierung er- 
sparen und dafiir die Produktmatrix bei einem solchen Korrekturverfahren anwenden, 
wie es auch beim additiven Verfahren zur Korrektur bentitzt wird. 


Liegen jedoch alle wiederzugebenden Farben F=F,+AF in der Umgebung 


einer Farbe Fy, so liBt sich die erwihnte Produktmatrix realisieren. Fiir jede 
Koordinate hat man den 


log (C + X) = log (C (1 + x)) = log C + log (1 +2) =logC +2 


zu bilden, wobei x = X/C klein ist. Bentitzt man natiirliche (Nepersche) Logarithmen 
und bezieht alle Dichten (auch D) auf diese, so hat man fiir die Koordinaten Ausdriicke 
der Form K+ 2. Man hat dann fir log (V A, Af) = = log (V F4) die Form 
Ke W-AF*. Die Matrix W unterscheidet sich von der V-Matrix dadurch, daB 


jede Zeile mit je einem Faktor multipliziert ist. Hier hat man das Produkt M-W 
zu realisieren. 


Anwendung eines Bildtransformators 

Die mit der Strahlenteilungskamera (oder nacheinander) aufgenommenen und 
(mit dem resultierenden Gamma G = 1 im Kontakt) kopierten Diapositive S werden 
von drei Photozellen PZ wie bei der Bildtelegraphie abgetastet. Der Vektor der 
drei Photozellenspannungen ist dann F“, der in das elektronische Rechengerat hinein- 
gesteuert wird. Das Gerat liefert drei Spannungen zum Steuern dreier Kopierbelich- 
tungslampen LQ, so da ihre Helligkeiten 7” = antilog (— 8’) entsprechen. Der 
Rechenautomat hat also 


T” = antilog (— M (— log (V F4)) (20) 


zu bilden, wobei die Minuszeichen wegfallen. Man erhalt so vorerst fiir das subtraktive 
Verfahren korrigierte Negative. Kine Kontaktkopierung derselben liefert die nun 
voll korrigierten Diapositive, deren Silber dann durch Farbstoff ersetzt wird. Bei diesem 
Kontaktkopieren muf aber im allgemeinen mit dem resultierenden Gamma gleich 
Kins kopiert werden. Werden davon abweichende Gammawerte beniitzt, ist dies 


gleichbedeutend mit einer Anderung von H. 


Die Abb. 4 zeigt das Blockschema des zum Rildiransfomaator: gehorenden elek- 
tronischen Rechengerates. Die bei der Durchleuchtung und bildpunktweisen Abtastung 


Uber die richtige Farbenwiedergabe 109 


der unkorrigierten Auszugsdiapositive erhaltenen Photozellenspannungen (denen der 
Hilfsvektor F“ zugeordnet ist) werden links in das die Matrix V realisierende Gerat 
Abb. 5 hineingesteuert, das die Farbraumkoordinatenumwandlung bewirkt. Die 
Konstanten V,;, kénnen durch neun Justierdrehknépfe eingestellt werden. Dieses 


(3 =F Fexp(U (log (V FAW=4(T 4) =O 
Abb. 4. Bildtransformator 
V-Gerat liefert Spannungen, die fiir die Kopierlampensteuerung zum Kopieren der 


fiir das additive Wiedergabeverfahren korrigierten Ausziige dienen. Daf die Hellig- 
keit einer Lampe im allgemeinen nicht der Lampenspannung U proportional ist, 


V-Gerat Viz Vig V3 


Abb. 5. Matrixnachbilder 


muB noch durch ein eigenes Gerat f (F”) gebildet werden, das weiter unten gezeichnet 
ist. Erst dieses f-Gerat steuert dann die Belichtungslampen, die mit dem tblichen 
Symbol fiir Gliihlampen eingezeichnet sind (trigheitslosere Glimmlampen waren 
hier vorzuziehen). Hat man bereits fiir das additive Verfahren korrigierte Diapositive, 
so kann man diese statt der unkorrigierten durch Photozellen abtasten und die Photo- 
zellenspannungen bei den Klemmen F® hineinsteuern und sie statt der aus dem 
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V-Geriit entnommenen dem Logarithmiergeriit log zufiihren. Nach dem log-Gerat 


kommt wieder ein Matrizennachbilder, und zwar fiir M (mit neun Justierdrehknépfen 
fiir die M,,). Dann kommt ein Antilogarithmiergeriit (antilog), das Spannungen 


liefert, zu denen der Hilfsvektor 7’ gehort. Zum Schlu8 kommt wieder ein /-Gerat, 
das Belichtungslampen steuert fiir die fiir das subtraktive Verfahren korrigierten 
Ausziige. 18 Zahlenwerte sind justierbar, da ja zwei Matrizen zu realisieren sind. 
Je nach den Aufnahmegrundfarben, den Wiedergabegrundfarben und den Farb- 
stoffen des subtraktiven Verfahrens kann nun der Bildtransformator richtig justiert 
werden. Weitere Hilfseinrichtungen kénnen das unmittelbare Abtasten der ur- 
spriinglichen Auszugsnegative erméglichen. Ebenso wiire es auch méglich, unmittel- 
bar durch den Bildtransformator die korrigierten Diapositive (statt der entsprechenden 
Negative) zu belichten. (Erster Bildtransformator von Siemens & Halske, F. Fischer 
1936.) 


Auch fiir den Dreifarbendruck (mit und ohne zusatzlichen Graudruck), der in 
gewissem Sinne ein gemischt additiv-subtraktives Verfahren darstellt, gibt es Bild- 
transformatoren, die die so lastige Retusche iiberfliissig machen (so daB sie hochstens 
fiir die letzte Feinbearbeitung in Frage kommt). Der Bildtransformator leistet mehr 
als das Maskierungsverfahren (das erst in letzter Zeit groBere Verbreitung finden 
konnte). 


Es wurden insbesondere vom Verfasser mehr oder weniger automatisch arbeitende 
Mehrfarbenverfahren entwickelt, die eine wesentliche Erweiterung des Bereiches der 
richtig wiedergebbaren Farben zu noch hoheren Sattigungen hin ermdglichen, so daB 
bereits eine hinreichend richtige Wiedergabe des gesamten sichtbaren reinen Spektrums 
ermoglicht wurde. 


Wie die obigen mathematischen Untersuchungen ergeben haben, ist eine richtige 
subtraktive Wiedergabe bei Anwendung eines trichromatischen Betrachtungsweih 
(bei nicht zu satten Farben) tatsichlich moglich. Bei emem gewohnlichen WeiB wire 
ein ganz wesentlich komplizierteres Rechengeréit beim Bildtransformator nétig, das 
wenig Aussicht auf Kinfiihrung in die Praxis hatte. Der einfachere Bildtransformator 
ermoglicht bei zweckmaBiger, Kompromissen entsprechender Justierung auch beim 
gewohnlichen Weil, insbesondere der. Xenonlampe, eine meist hinreichend gute 
Wiedergabe, zumindest eines Vollfarbenspektrums (von Oberflichenfarben). 


Bei den tiblichen Mehrschichtenfarbfilmen treten noch andere Fehlerquellen auf, 
die einen Bildtransformator sehr komplizieren wiirden. Eingebaute Farbmasken 
sind hier nur halbe Mafinahmen. (S. Osterreichische Photozeitung 1956/11, S. 10, 12.) 


Anmerkung bei der Korrektur: Beim Neuentwurf der Farbnormen ist die 
Méglichkeit vorgesehen, beim Farbvektor statt der Uberquerung die Einklammerung 
anzuwenden. Fir die Matrix kame dann die doppelte Klammer (oder eine eckige) 
in Betracht. Z. B. ware dann (Ff) ein Farbvektor und ((V)) eine Matrix. Fiir die 
Spektralkurve wird statt des Buchstabens S der Buchstabe P vorgeschlagen, da 8 
fiir die Stahldichte vorgesehen ist. 


Literatur 


W. Schultze: Farbenphotographie und Farbenfilm (enthalt ein reichhaltiges Literatur- 
verzeichnis). 

OTF, 1950, Hefte 5/6 und 7/8, Dreifarbenfernseher mit richtiger Wiedergabe. 

Photographische Korrespondenz: 1954, H. 10 und 11 ,,Korrigiertes subtraktives Farbfilm- 
verfahren“; 1955, H.9 ,,Farbenreproduktionsverfahren mit m Aufnahme- und n Wiedergabe- 
grundfarben‘; 1956, H. 3 additives Verfahren. 
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Zu den Tanischen Grenzschichten 
Von H. Gértler und H. Witting, Freiburg i. Br. 
Mit 8 Textabbildungen 


Zusammenfassung. Eine neue Methode! zur Berechnung stationirer, ebener und inkompres- 
sibler laminarer Grenzschichten bei allgemeinem auBerem Druckverlauf mit Hilfe einer formal 
exakten Lésung dieses allgemeinen Problems der Grenzschichttheorie in Gestalt einer unend- 
lichen Reihe in besonders geeigneten Variablen wird dargestellt. Es wird gezeigt, daB diese Methode 
sich erneut besonders gut bewahrt in Anwendung auf eine neuerdings von I. Tani untersuchte 
Klasse von Grenzschichtstromungen zu auBeren Geschwindigkeitsverteilungen der Gestalt 
U(x) = U) + u, x”. Nachdem bisher nur der Fall n = 1 gelést worden war (L. Howarth), 
ist die Abhangigkeit des Grenzschichtverlaufs und insbesondere der Ablésestelle von dem Expo- 
nenten von besonderem praktischem Interesse bei der Beurteilung des Einflusses von Druck- 
anstiegen auf den Strémungsverlauf. 

Tani benutzte in enger Anlehnung an Howarth eine Reihenentwicklung vom Blasiusschen 
Typus (das hei8t nach Potenzen der iiberstrémten Wandbogenlange) und eine auf plausiblen An- 
nahmen beruhende approximative Fortsetzungsmethode dort, wo die verfiigbaren Glieder der 
abgebrochenen Reihe wegen zu langsamer Konvergenz keine brauchbare Approximation mehr 
hefern. 

Es wird gezeigt, dafi auch im Falle der Tanischen Grenzschichten die Methode der neuen 
Reihe der Methode der Blasiusschen Reihe tiberlegen ist. Sie liefert mit wesentlich weniger 
Gliedern eine ebenso weitreichende gute Approximation, wie sie von Tani erzielt worden 
ist. Dort, wo auch die nach der neuen Reihe (beim gegenwirtigen Stand der Vertafelung der 
universellen Koeffizientenfunktionen) erzielbare Approximation stromabwirts nicht mehr aus- 
reicht, wird eine sorgfaltige Fortsetzungsrechnung mit der von den Verfassern frither entwickelten 
Differenzenmethode bis zur jeweiligen Ablésungsstelle angeschlossen. Es ergibt sich Ablésung 
der Grenzschicht durchweg etwas weiter stromabwarts als nach dem Howarth-Tanischen Fort- 
setzungsverfahren, jedoch betragt die Abweichung nur wenige Prozent der von der Vorderkante 
an bis zur Ablésestelle insgesamt tiberstromten Wandbogenlange, und sie ist von keinem praktischen 
Interesse. 


I. Einleitung 


Fiir das allgemeine Problem der Berechnung laminarer Grenzschichten bei beliebig 
vorgegebenem auferem Druckverlauf wurde kiirzlich eine neue Methode entwickelt?. 
Sie beruht auf einer formal exakten Loésung dieses allgemeinen Randwertproblems 
der Grenzschichttheorie — zunachst fiir ebene und inkompressible Stromungen — 
in Gestalt einer unendlichen Reihe. Die durch die Wahl besonders geeigneter Variablen 
erzielten wesentlichen Vorteile der Methode wurden an einer gréBeren Zahl von 
Erprobungsbeispielen bestatigt!,*. Insbesondere ist der Arbeitsaufwand bei der 
Anwendung dieser Methode nicht gré8er als derjenige der — wegen ihrer Kinfachheit 
in der Praxis bisher bevorzugten — Parametermethoden vom Typus des Pohlhausen- 
Verfahrens. Diese Verfahren sind jedoch nicht streng begriindet, vom Ansatz her 
ungenau und erfahrungsgemaB unzuverlassig. 

Das neue Verfahren ist nicht nur anwendbar auf umstro6mte Konturen mit runder 
Nase, sondern allgemeiner auf Profile mit beliebigem Keilwinkel an der Vorderkante. 
Insbesondere sind auch — Vorderkantenwinkel Null — spitze Profile oder, was grenz- 
schichttheoretisch auf dasselbe hinauskommt, Platten mit beliebigem auBerem Druck- 
verlauf der Berechnung zuginglich. Neuerdings ist das Verfahren auch auf Probleme 
der Grenzschichtbeeinflussung durch Absaugung bei allgemeiner kontinuierlicher 
Verteilung der Absaugegeschwindigkeit langs der Wand erweitert worden’. 


1H. Gortler: A new series for the calculation of steady laminar boundary-layer flows. J. 
Math. Mechan. 6, Nr. 1 (1957). 
2 H. Gortler und H. Witting: Einige laminare Grenzschichtstromungen, berechnet mittels 
einer neuen Reihenmethode. Zschr. angew. Math. Physik, erscheint in Kirze. 
3 H. Gortler: On the calculation of steady laminar boundary-layer flows with continuous 
suction. J. Math. Mechan. 6, Nr. 2 (1957). ‘ 
8* 
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Die vorliegende Untersuchung beschiftigt sich mit einer gewissen Klasse von 
Grenzschichten an lings angestromten Platten oder Profilen mit spitzer Vorderkante. 
Neben der seit den Anfangen der Grenzschichttheorie bekannten exakten Blasius- 
schen Lisung fiir die Platte bei konstantem auBerem Druck war bis in neuester Zeit 
nur eine einzige weitere Plattenstrémung exakt behandelt worden, naimlich die Grenz- 
schicht bei linear mit der tiberstr6mten Tiefe anwachsendem Druck. L. Howarth’ 
hat diese Strémung mit Hilfe einer unendlichen Reihe nach Potenzen der tiberstromten 
Tiefe berechnet. (Neben der Lésung allgemeinerer Plattenstrémungsprobleme in ?* 
und? ist in1 dieser Fall des linearen Druckanstiegs nach der Methode der neuen 
Reihe erneut berechnet worden, weil die Howarthschen Ergebnisse eine willkommene 
Vergleichsméglichkeit boten.) Die Howarthsche Reihe ist vom gleichen Typus wie 
die in der Literatur wohlbekannte Blasiussche Reihe fiir Grenzschichten an Profilen 
mit runder Nase. Erst neuerdings ist uns bekanntgeworden, dai A. N. Tifford 
in einer nicht im Druck erschienenen Untersuchung® unter anderem die Howarthsche 
Reihe verallgemeinert hat fiir Plattenstromungen mit allgemeinem auBerem Druck- 
verlauf. Veréffentlicht wurde jedoch vor kurzem eine Verallgemeinerung der Howarth- 
schen Reihe auf eine Klasse einfacher nichtlinearer Druckverteilungen, die von 
I. Tani® 7 entwickelt und der numerischen Auswertung zuganglich gemacht worden 
ist. Diese einparametrige Schar von Lésungen hat uns eine sehr willkommene Még- 
lichkeit geboten, die Leistungsfihigkeit der Methode der neuen Reihe zu vergleichen 
mit jener der von Tani verwendeten Reihe vom Blasiusschen Typus (Entwicklung 
nach Potenzen der tiberstr6mten Tiefe). 

Die Klasse der Tanischen Grenzschichten stellt auf den ersten Blick eine fiir die 
Anwendung der Methode der neuen Reihe besonders ungiinstige Klasse von Strémungen 
dar, wie weiter unten im einzelnen ausgefiihrt werden wird. Es hat sich aber gezeigt, 
daB sich die neue Reihe bei der Anwendung auf diese Stromungsklasse durch be- 
sonders giinstiges Konvergenzverhalten ausweist, so da der angedeutete Nachteil 
nicht zum Tragen kommt. Die neue Reihe erweist sich damit der Reihe vom Blasius- 
schen Typus erneut eindeutig tiberlegen. 


II. Kurze Darstellung der Methode der neuen Reihe 


Um die Gesichtspunkte des Vergleiches verstandlich zu machen, mu8 hier zunachst 
kurz die Methode der neuen Reihe beschrieben werden. Wir beschrinken uns dabei 
auf die wesentlichen Ziige. 

Das Randwertproblem der — ebenen und stationiren — Strémung in einer 
laminaren Grenzschicht einer inkompressiblen Fliissigkeit ist durch Bewegungs- 
und Kontinuitaétsgleichung 


U Ug + 0 Uy = U(x) U(x) + ¥ Uy, | 


te + 0, = 0 j iD) 
und die Randbedingungen 
== Oy Tity vay 0, 
essa)! eee 
u—>U(x) fir y\Re->oo ae (2) 


* L. Howarth: On the solution of the laminar boundary layer equations. Proc. Roy. Soc. 
London, Ser. A 164, 547—579 (1938). 

BANE Tifford: Heat transfer and frictional effects in laminar boundary layers. Part 4: 
Universal series solutions. Wright Air Development Center Technical Report 53-288 (August 1954). 

° I. Tani: On the solution of the laminar boundary layer equations. J. Physic. Soc. J apan 4, 
149—154 (1949). 

* I. Tani: On the solution of the laminar boundary layer equations. ,,50 Jahre Grenzschicht- 


forschung“‘, herausgegeben von H. Gértler und W. Tollmien. Braunschweig: Verlag Friedr. 
Vieweg u. Sohn. 1955. 
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im wesentlichen gegeben. Es kommt im allgemeinen Fall eine Anfangsbedingung 
(Kinlaufbedingung im Anfangsquerschnitt « = 0) hinzu, die jedoch bei Problemen 
mit singulirer Vorderkante x = 0, y = 0 durch die Vorschrift des (bekannten) Ver- 
haltens der Lésung in der Umgebung von Vorderkanten (Staulinien, Plattenspitzen) 
beim Losungsansatz ersetzt wird. In (1) und (2) bedeuten a die Wandbogenlange, 
gemessen vom eben genannten vorderen singuliéren Punkt des Stromungsfeldes, 
y den senkrechten Abstand von der Wand, wu und v die Geschwindigkeitskomponenten 
in der Grenzschicht in x- bzw. in y-Richtung, U(x) die vorgegebene auBere Geschwin- 
digkeitsverteilung, » die konstante kinematische Ziahigkeit der Fliissigkeit und Re 
schlieBlich eine kennzeichnende Reynoldssche Zahl der Strémung. 

Fiir die Methode der neuen Reihe ist die Einfiithrung der folgenden neuen unab- 
hangigen Koordinaten grundlegend: 


: x ls 
|U(2)de, 4 =—@)¥ . Uteyas} (3) 
0 / 


Ist (x, y) die vermdge y, =u, y, = — »,. v(x, 0) = 0 definierte Stromfunktion 
der Grenzschichtstr6mung, durch deren Einfiihrung die zweite Gl. (1) (Kontinuitats- 
gleichung) befriedigt wird, so wird ferner als neue abhaingige Variable die Funktion 


/ ch we 
F(E, 0) = yla, y) | {2>( U(«) aa| (4) 
/ 0 
benutzt. In den Variablen (3), (4) geht das Grenzschichtproblem (1), (2) iiber in 
Fung +F Fy, + B(E){1 —Fy3} = 28 (FFs, — Fe Fyn} G) 
mit den Randbedingungen 
F(é, 0) = 0, Pgs te; 0) = 0, F, ig, co) =1 (0 = E(ap)). (6) 


In dieser Formulierung erscheinen die Vorgaben des jeweiligen speziellen Problems 
(auBere Geschwindigkeit und Stoffkonstanten) in (6) tiberhaupt nicht — die Rand- 
bedingungen sind ,,universell‘‘ — und in der Differentialgleichung (5) erscheinen 


sie nur in der Koeffizientenfunktion 
uv 


B(é) = 2 U"(x) | U(x) da / U%(a). (7) 
0 
Wegen ihrer zentralen Bedeutung fiir das Verfahren wurde diese Funktion als ,, Haupt- 
funktion der Grenzschicht“ bezeichnet. 
Es wird vorausgesetzt, daB die Hauptfunktion eine im Grundintervall analytische 
Funktion von & ist, welche eine Reihenentwicklung 


B(é) = Bu é (8) 


gestattet. Nach (7) bedeutet dies, daB fiir die AuBere Geschwindigkeitsverteilung U(z) 
eine konvergente Entwicklung der Gestalt 


U (a) Sa? {ug a tt Pater) By 


; (9) 
mit m=,/(2— 6.) und mit wu +0 


méglich sein soll. Dabei kennzeichnet der Koeffizient B, die Geometrie des Zylinders 
an der Vorderkante: Es ist a B, der Keilwinkel der Kontur an der Vorderkante «= 0. 
Da wir uns in der vorliegenden Untersuchung nur fiir den Problemkreis mit spitzer 
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Vorderkante (8, = 0, somit m=) interessieren, sind hier nach (9) alle auBeren 
Geschwindigkeitsverteilungen zugelassen, die sich iiberhaupt in eine konvergente 
Potenzreihe : 

U(a) =u tue tune +t... mit u +90 (10) 
mit sonst beliebigen Koeffizienten uw, (k = 0, 1, 2,...) entwickeln lassen (0 S a S &p). 


Fiir die gesuchte) Lésungsfunktion F(§, 1) wird die Entwicklung 


EE, 0) = 2 Fale (11) 
k=0 

angesetzt, wie dies durch (8) und (5). nahegelegt ist. Kinsetzen in (5) ergibt nach 
Koeffizientenvergleich ein rekursiv zu losendes System von Differentialgleichungen 
fir Fy, F,, F>,... Zunachst zeigt sich, was nicht iiberrascht, daB Fy (7) die bekannte 
Blasiussche Funktion der Plattengrenzschicht bei konstantem Druck ist, jetzt jedoch 
als Funktion der allgemeineren Variablen 7 nach (3), die nur im Blasiusschen Sonder- 
fall U(x) = U, = konst. mit der von ihm benutzten Variablen tibereinstimmt, im 
iibrigen aber — was wieder nicht tiberrascht — in der Grenze « —> 0 in diese tiber- 
geht. (In anderen Worten: In naichster Umgebung der Vorderkante stimmt die Grenz- 
schichtlésung im allgemeinen Fall (10) approximativ mit der Blasiusschen Platten- 
grenzschicht tberein.) 


Die Randbedingungen (6) gehen tiber in 
F,(0) = Fy(0) = 0, Fy'(co) = 1; \ 
FiO) PF (0) Fy (cols Oe ene eete) 
Ein wesentlicher Vorteil der neuen Reihe gegeniiber der Reihenentwicklung vom 
Blasiusschen Typus wird hier erkennbar: Das Glied nullter Ordnung der neuen Reihe 
erfiillt die duBere Randbedingung w —> U(x) bereits exakt fiir alle x (alle €) des zugrunde 
liegenden Intervalls; die Blasiussche Reihe hingegen erfiillt diese Randbedingung 


durch Abgleichen, Glied fiir Glied, fiir 7 —> co mit den Gliedern der Reihe (10). (Fiir 
EKinzelausfiihrungen hierzu vel. 1.) 


(12) 


Zwischen den Koeffizienten 6, der Hauptfunktion und den Koeffizienten w, der 
vorgegebenen auBeren Geschwindigkeitsverteilung bestehen gem4B (7) und (3) einfache 
Relationen, die in! fiir praktische Bediirfnisse angegeben sind. Mit ein paar Multi- 
plikationen und Additionen sind damit die £;, ebenfalls gegeben. 


Setzt man nun fiir die Funktionen F;,(7) Linearkombinationen nach folgender 
Bauart an: 


f= Bi hi 
F, = BY hi + Bele, (13) 
Fy = BY fin + Bi Behe + Bs fz usw., 


so erweisen sich die Funktionen /..(7) als Lésungen universeller Randwertprobleme. 
Man kann sie somit ein fiir allemal bis zu einer ausreichend hohen Ordnung und mit 
fiir alle Anwendungen geniigender Genauigkeit vertafeln. Das ist fiir alle universellen 
Funktionen bis zur 5. Ordnung einschlieBlich geschehen, Neben den von der elektroni- 
schen Rechenanlage MARK IV (Harvard University) gelieferten umfangreichen Tafeln 
mit 7-Werten in 0(0°01)5-99, die Interessenten auf Wunsch zur Verfiigung gestellt 
werden k6nnen, sind Ausziige mit einer fiir praktische Bediirfnisse ausreichenden 
Anzahl von 7-Werten, namlich 0(0°2)6-0 bereitgestellt worden’, 

* H. Gortler: Zahlentafeln universeller Funktionen zur neuen Reihe fiir die Berechnung 


laminarer Grenzschichten. Bericht Nr. 34 der Deutschen Versuchsanstalt fiir Luftfahrt (1957) 
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Man erkennt aus (13): Die Anwendung der Methode in jedem Einzelfall besteht 
nur noch im Ausrechnen der Kombinationen (13) oder von deren Ableitungen unter 
Verwendung der Tafeln. Fiir die lokale Wandschubspannung bendtigt man, um ein 
Beispiel zu nennen, aus den Tafeln nur die Werte f..’’(0), bildet daraus nach (13) 
rasch die F’,'’(0) und hat damit mit den bis k — 5 verftigbaren Tafeln nach (11) eine 
Polynomapproximation 5. Grades fiir die Wandschubspannung. 


Der Wert des Verfahrens hingt von der Konvergenzgiite der Reihe (11) ab. In 
allen Beispielen von Plattenstrémungen mit 4uBerem Druckanstieg, die bisher unter- 
sucht worden sind! *, reichte diese Approximation aus, um mit guter Genauigkeit 
die Grenzschicht in ihrem gesamten Verlauf von der Vorderkante x = 0 bis un- 
mittelbar vor der Ablésungsstelle zu berechnen. 


Ilf. Die Tanischen Grenzschichten und deren Berechnung durch Tanis Reihen 
vom Blasiusschen Typus 


Zur Zeit der Untersuchungen von I. Tani® war die soeben kurz zusammengefaBte 
Methode der neuen Reihe noch nicht bekannt, und dementsprechend waren ins- 
besondere Plattenstromungen bei allgemeinem auBerem Druckverlauf noch nicht 
einer exakten Behandlung zugiinglich. Indem Tani duBere Geschwindigkeitsver- 
teilungen der Form 

Oe) ge Ue MN hg ste OS ite O) (14) 


seiner Arbeit zugrunde legte, war es sein Bestreben, tiber die zuvor allein exakt be- 
handelten Sonderfille n = 0 (Blasius) und » = 1 (Howarth) hinaus zu einer Klasse 
von nichttrivialen Grenzschichtstromungen vorzudringen, deren Verhalten wertvolle 
Aufschliisse bei allgemeineren Uberlegungen bieten kénnte. Fir w,,/u) <0 liegt 
Druckanstieg in Stromungsrichtung vor, und es interessierte, die Abhingigkeit der 
Ablosungsstelle vom Parameter » zuverlassig zu ermitteln. Tani hat die Falle 
nm = 2,4 und 8 numerisch ausgewertet. 


Die von Tani benutzte Lésungsmethode ist eine unmittelbare Ubertragung der 
von Howarth entwickelten Mittel. In Anlehnung an Howarth werden an Stelle von x 
und y die unabhangigen Veranderlichen 


F = 8Uy XU, t= (y/2) Vuo|v x (15) 


eingefiihrt. Der Lésungsansatz einer Potenzreihe vom Blasiusschen Typus, also nach 
Potenzen von x, reduziert sich bei der speziellen Vorgabe (14) auf eine solche nach 
Potenzen von 2”, also von o. Fiir die Stromfunktion y(x, y) der Grenzschicht ange- 


schrieben, lautet die Reihe formal: 
co 


p(x, Y)[Vv © Uy = D/ tet) oF. (16) 


Es sei hier angemerkt, daB fiir den besonders interessierenden Fall des Druckanstieges 
iiber das Grundintervall 0 < x S x, wegen wy/uv) < 0 in (15) die zugehérigen Werte 
von o negativ sind. 

Die Einfiihrung des Ansatzes (16) in die Grenzschichtdifferentialgleichung liefert 
zusammen mit den Randbedingungen ein rekursiv zu losendes System von Randwert- 
problemen fiir die Funktionen ¢,. Dabei ist, wie zu erwarten, f)(t) die bekannte 
Blasiussche Funktion der Plattengrenzschicht bei konstantem Druck. (Wegen der 
etwas verschiedenen Variablenwahl ist zwar t) formal Losung derselben Differential- 
gleichung wie Fy, aber die um den Faktor 2 verschiedenen Variablen bedingen, dab 
z. B. t)'(t) = 2 Fo'(n) wird, wo Striche Ableitungen nach den zugehérigen Variablen 
bezeichnen.) f(t) erfiillt die Randbedingungen 4)(0) = t)'(0) = 0, fo'(co) = 2. 
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In den Differentialgleichungen fiir die Funktionen ¢,(7) fiir & = 1, 2,3,... er- 
scheint der Parameter » der zugrunde liegenden auBeren Geschwindigkeiten (14) 
explizit in den Koeffizienten. Somit ist das System der Funktionen ¢, (bid, 
3,...) fiir jedes spezielle n erneut zu berechnen: ¢, = t,(7; ”) fir 4 = 1,2; /3).4% 
Das zwingt zur Beschrinkung der numerischen Auswertungen auf eine Auswahl 
von n-Werten. 

Die Randbedingungen fiir die ¢, lauten gemeinsam fiir alle n 


(0) 4 (0) Oodle: ead 2e Sane 
£5t( doy aah, 4a OE et 


Man erkennt, daB die auBere Randbedingung wu — U, das heiBt u/uy — 1 + o/8 fir 
t —> co von der Blasius-Reihe (16) fiir die vorgegebene Klasse (14) von auBeren Ge- 
schwindigkeiten bereits mit den zwei ersten Termen der Reihe exakt erfiillt wird, 
wihrend die Terme hodherer Ordnung am duBeren Rande der Grenzschicht keine 
Beitrage liefern. Somit ist beim vorliegenden Problemkreis die neue Reihe hinsichtlich 
der Erfiillung der 4uBeren Randbedingung (exakt bereits mit dem ersten Term der 
Reihe) in keinem wesentlichen Vorteil gegeniiber der Blasius-Reihe von Tani; denn 
bei der praktisch notwendigen Beschriinkung auf abgebrochene Reihen wird man 
immer wesentlich tiber zwei Glieder hinausgehen konnen. 

Fir die Fille n = 2,4, 8 ermittelte Tani die folgenden Funktionen numerisch: 


t.(t:2) fir k <6, 
bAlry4). 2 Tir ee 
ty(t5 8) far: 


In Potenzen von x ausgedriickt, war damit die Reihe (16) fiir n = 2 mit 7 Gliedern 
bis zum Gliede mit x! verfiigbar; entsprechend fiir n = 4 mit 6 Gliedern bis zum 
Gliede mit 2°, fiir n = 8 mit 5 Gliedern bis zum Gliede mit x2. Die Konvergenz- 
giite erwies sich jedoch nicht als ausreichend, um damit die Grenzschicht in ihrem 
gesamten Verlauf bis zur Ablosungsstelle zu ermitteln. Tani muBte daher ein nicht 
streng zu rechtfertigendes Fortsetzungsverfahren anwenden; er schloB sich dabei 
wieder eng an das Vorbild von Howarth an. 

Uber die numerischen Ergebnisse wird im folgenden Abschnitt im Zusammenhang 
des Vergleiches mit den entsprechenden Ergebnissen der Methode der neuen Reihe 
berichtet. 


(17) 


IV. Berechnung der Tanischen Grenzschichten mit Hilfe der neuen Reihe und Vergleich 
der Leistungsfahigkeit beider Verfahren 


Bei Anwendung der Methode der neuen Reihe hat man es bei allen auBeren Ge- 


schwindigkeiten der Gestalt U(x) = »’u,a* mit uy +0, insbesondere auch bei 


allen Tanischen Verteilungen U(x) = uw) + Uy 2” (n = 0, 1, 2, 3,...) formal mit ein 
und derselben Reihe zu tun, nimlich mit der Reihe, welche das System universeller 
Funktionen mit 2, = 0 besitzt. Dieses System ist gegenwirtig (siehe oben) bis zur 
Ordnung 5 einschlieBlich numerisch bereitgestellt. 

In den folgenden Ausfiihrungen bzw. Figuren erscheinen die dimensionslosen 


GroBen 
c=allL, y=yRer/L, wo L=(—U/u,)'/n und Re =u Lr; 


U(x, Y) = U(x, y)/Ug, U(x) = U(x) [9 = 1 — 2”; (18) 
E=£/Re, B, =f, Ret (somit: p, & = B, é). 
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Die neue Reihe ist nach Potenzen von é, somit nach (3) nach Potenzen von 


E= x — o"+1)(m +1) (19) 
geordnet. 

Die spezielle Gestalt von U(x) = 1 — x" hat zur Folge, daB von den Koeffizien- 
ten 6, der Hauptfunktion 8(€) nur B,,, By, Bsn, ... von Null verschieden sind. Dadurch 
erhalten sehr viele universelle Funktionen den Koeffizienten Null in der Losungsreihe. 
Angeschrieben fiir die Geschwindigkeitsprofile u(a, y) erhalt man 


WO i, Bf a BP, Ban fan) be? (20) 


Hieraus ergibt sich unmittelbar die einleitend angekiindigte Tatsache, daB die Klasse 
der Tanischen Grenzschichten fiir die Anwendung der neuen Reihe beim gegen- 
wartigen Stand der Vertafelung universeller Funktionen eigentlich ein ausgesprochen 
ungtnstiges Ubungsfeld bietet. Da fiir die allgemeine Reihe und damit auch fiir die 
spezielle Reihe (20) universelle Funktionen nur bis zur Ordnung 5 vertafelt vorliegen, 
haben wir namlich fiir » = 1 noch 6 Glieder von (20), fiir » = 2 bereits nur 3 Glieder 
von (20), fiir » = 3,4 und 5 sogar nur 2 Glieder numerisch verfiigbar. 

Natiirlich hatte man bei Beschriinkung auf eine Untersuchung allein der Tanischen 
Grenzschichten auf die Fiille der in (20) nicht in Erscheinung tretenden Funktionen 
verzichten und dafiir in (20) fiir spezielle n-Werte die numerische Vertafelung zu weit 
héheren Ordnungen mit angemessenem Aufwand fiihren kénnen. (Es werden ja fiir 


jedes n jeweils neben Fy nur fy, fnns fons fnnn +--+ bendtigt.) Das wiirde dem Vorgehen 
von Tani entsprechen. Da uns jedoch der numerische Zugang zu Problemen mit 
allgemeinem U(x) =u) +u,% + u,u? +... interessierte, ist das bisher nicht 
geschehen. 


Uberraschenderweise zeigt es sich, da wir fiir die Falle n = 2 und » = 4 von 
Tani bereits mit unseren wenigen Reihengliedern die Grenzschicht ebenso weit strom- 
abwarts berechnen kénnen, wie Tani mit wesentlich mehr Gliedern der Blasius-Reihe. 
Die Ergebnisse stimmen im Rahmen der Zeichengenauigkeit vollkommen miteinander 
tiberein. Dies ist zunichst eine sehr erfreuliche gegenseitige Bestatigung der 
numerischen Richtigkeit beider Ergebnisse. Dariiber hinaus wird aber die wesentlich 
bessere Approximationsméglichkeit mittels der neuen Reihe — namlich mit erheblich 
weniger Gliedern als mit der Blasius-Reihe — tiberzeugend demonstriert. 

Da bei den heute mit Hilfe der neuen Reihe numerisch in gréBerer Zahl}:? be- 
rechneten und in weiterer Fiille bequem zuginglichen Stromungen dieser Art eine 
weitere Vertafelung der in (20) auftretenden Funktionen nicht lohnend erscheint, 
mu8ten auch wir ein Fortsetzungsverfahren anwenden, um von den Ergebnissen 
der wenigen Glieder der neuen Reihe (und damit der groBeren Zahl von Gliedern der 
Tanischen Reihe) weiter stromabwarts bis zur Ablésungsstelle zu gelangen. Wir 
schlossen uns nicht dem Howarth-Tanischen Fortsetzungsverfahren an, sondern be- 
nutzten die von uns entwickelte Differenzenmethode der schrittweisen Fort- 
setzung® 1°, Diese ist zugleich neben sonstigen Proben (sog. ,,erste Grenzschicht- 
bindung‘‘ etwa) eine willkommene Kontrolle der Brauchbarkeit der von der Reihe 
her gelieferten Approximation. Indem man wesentlich weiter stromaufwarts zuriick- 
geht mit dem Hinsatz der Differenzenrechnung, als dies notwendig ist, kann man strom- 
abwirts verfolgen, wie weit die vom Reihenverfahren gelieferte Approximation 


9 H. Gortler: Ein Differenzenverfahren zur Berechnung laminarer Grenzschichten. Ing.- 
Arch. 16, 173—187 (1948). 

10 H. Witting: Verbesserung des Differenzenverfahrens von H. Gértler zur Berechnung 
laminarer Grenzschichten. Diss. Freiburg i. Br. (1953); Z. angew. Math. Physik 4, 376—397 
(1953). 
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brauchbar bleibt und von wo an stromabwirts ins Gewicht fallende Abweichungen 
einzusetzen beginnen. Unsere weniggliedrigen Ergebnisse sind somit nicht nur in 
Ubereinstimmung mit den Ergebnissen der mehrgliedrigen Tanischen Reihe, soweit 
diese nach Tani eine brauchbare Approximation bieten, sie sind vielmehr auch mit 
Hilfe der Differenzenmethode unabhingig von den Tanischen Ergebnissen sorefaltig 
kontrolliert und zuverlissig fortgesetzt worden. 

Neben den Tanischen Fallen » = 2 und n = 4 und dem Howarthschen Fall n = 1 
(iiber den wir ausfithrlicher schon in ! berichtet haben) haben wir auch die Fallen = 3 


Q7 02 03 t ) G2 Q4 06 z 


Abb. 1. Verlauf der lokalen Wandschubspannung 


Abb. 2. Verlauf der lokalen Wandschubspan- 


fiir den Fall n =|2 (U = 1— 2). R und Fortsetzung 
R, (ausgezogene Kurve) Ergebnis der neuen Reihe 
(R) und des Differenzen-Fortsetzungsverfahrens 
(R,). R, Ergebnis der Approximation 4. Ordnung 
in € (3 Glieder der neuen Reihe) dort, wo sie keine 
ausreichende Naherung mehr bietet. R, und Ry ent- 
sprechende Ergebnisse der Approximationen 2. Ord- 
nung (2 Glieder der neuen Reihe) und 0. Ordnung 
(1 Glied), 7, Tanische Fortsetzung dort, wo die 
Approximation durch die ersten 7 Glieder seiner 
Reihe (die bis dahin mit R tbereinstimmt) un- 
brauchbar wird 


nung fiir den Fall n = 4 (U = 1 — 2). Be- 

schriftung entspricht jener von Abb. 1. Hier 

ist jedoch R, nur eine 2gliedrige Approxi- 

mation durch die neue Reihe. Sie leistet das 

gleiche wie die ersten 6 Glieder der Tanischen 
Reihe 


und » = 5 ausgewertet. Wir berich- 
ten nun iiber Einzelheiten der nume- 
rischen Ergebnisse. 

In Abb. 1 ist zunadchst fiir den Fall 


nm = 2, somit fir U = 1 — 2?, der Ver- 
lauf der ortlichen Wandschubspannung dargestellt, sowohl in der nach dem gegenwirtigen 
Stand der Vertafelungen erreichbaren 3gliedrigen Approximation 4. Ordnung in é 
(Kurve R,), als auch in der vor der Ablésungsstelle verbesserten Berechnung mit 
Hilfe des Differenzen-Fortsetzungsverfahrens (ausgezogene Kurve R,). Zur Demon- 
stration der Konvergenzgiite sind ferner die Approximationen durch das 1. Glied 
der neuen Reihe allein (Approximation nullter Ordnung in £, Kurve R,) und durch 
die zwei ersten Reihenglieder (Approximation 2. Ordnung in é, Kurve R,) eingetragen. 
(Die Kurve R, gibt zugleich den Verlauf der lokalen Wandschubspannung fiir die 
Blasiussche Platte, das hei®t bei konstantem auferem Druck, wieder.) SchlieBlich 
zeigt in Abb. 1 die Kurve 7’, das Ergebnis des Tanischen Fortsetzungsverfahrens 
dort wieder, wo das Tanische Reihenergebnis keine brauchbare Approximation mehr 
lieferte. Das Tanische Fortsetzungsverfahren stiitzt sich auf plausible, aber unbewiesene 
Annahmen iiber die vernachlassigten Glieder seiner Reihe mit dem Anspruch, eine 
brauchbare Approximation zu liefern. Man erkennt: Dort, wo unsere Approximation R, 
mit drei Gliedern der neuen Reihe unbrauchbar wird, das hei&Bt von R; abzuweichen 
beginnt, muBte auch Tani mit seiner Fortsetzungsmethode einsetzen, um das bis dorthin 
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Abb. 3. Verlauf der lokalen Wandschubspannung fir alle ausgewerteten Falle (n = 1, 2, 3, 4 


und 5). Ergebnis der neuen Reihe und der Fortsetzung (R;) mit Hilfe des Differenzen- 
verfahrens. ------- Ergebnis der neuen Reihe (R) mit den ersten 6, 4, 2, 2 und 2 Gliedern im 
Bereich zu langsamer Konvergenz. -.-.-.- Ergebnis der Howarthschen (H, fur n = 1) und der 


Tanischen (7', fir n = 2 und 4) Fortsetzung dort, wo ihre Reihe mit 9 (n = 1), 7 (n = 2), 
6 (n = 4) Gliedern keine ausreichende Approximation mehr liefert 


he Se Se eS Ee 
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Abb. 4. Verlauf der Verdrangungsdicke 6* fir n = 1, 2, 3, 4 und 5. Ergebnis der neuen 
Reihe und der Fortsetzung mit Hilfe des Differenzenverfahrens. -------- Ergebnis der Approxi- 
mation mit 6, 3, 2, 2, 2 Gliedern der neuen Reihe im Bereich zu langsamer Konvergenz 


mit unserem Ergebnis R, tibereinstimmende Ergebnis der ersten 7 Glieder seiner 
Reihe zu korrigieren. (Leider gibt Tani nur seine endgiiltigen Resultate an; man kann 
seinen Angaben nicht entnehmen, wie seine 7gliedrige Approximation nach Unbrauch- 
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Abb. 5. Verlauf der Impulsverlustdicke # fir n = 1, 2, 3, 4 und 5. ———— und ------- wie in 
Abb. 4 
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Abb. 6. Geschwindigkeitsprofile an einer Auswahl von festen Wandstellen % fiir n = 2, 3, 4 und 5 


barwerden verlauft.) Somit leisten hier die drei ersten Glieder der neuen Reihe das- 
selbe wie die sieben ersten Glieder der Tanischen Reihe vom Blasiusschen Typus. 

Im wbrigen erkennt man aus Abb. 1, da8 unsere mit erheblicher Sorgfalt mit 
der Differenzenmethode ermittelte Fortsetzung R, eine etwas spitere Ablosung 
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(x ="0°290) hefert als die Tanische approximative Fortsetzung 7’, (x = 0°271). Doch 
ist die relative Abweichung (zirka 6°5°%%) wohl nur von mathematischem Interesse. 
(In Nachbarschaft der Ablésungsstelle sind die durch die Vernachlassigungen der 
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Abb. 7. Geschwindigkeitsverteilung in konstanten Wandabstaénden lings der Wand fir die 
Falle n = 2 und 4. ———— Ergebnisse der neuen Reihe und der Fortsetzung mit der Differenzen- 
methode. --~----- Tanis Ergebnisse dort, wo sie von unseren abweichen 


Grenzschichttheorie bedingten Fehler vermutlich im allgemeinen von der gleichen 
GréBenordnung einiger Prozente. Selbst wenn eine genauere experimentelle Er- 
mittlung der Ablosungsstelle vorlige, kénnte sie daher keine Entscheidung fiir das 
eine oder das andere der oben genannten Zahlenergebnisse liefern.) 
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Abb. 2 zeigt die entsprechenden Ergebnisse fiir das zweite Tanische Beispiel, 
n = 4, somit U = 1— x. Ganz entsprechend liefern die hier verfiigbaren ersten 
2 Glieder der neuen Reihe (Approximation 4. Ordnung in é, Kurve &,) Abweichungen 
von der Differenzenfortsetzung R, von jener Stelle an stromabwirts, wo auch Tani 
mit seiner Fortsetzung 7’; einsetzen muB, um das unbrauchbar werdende Ergebnis 
der ersten 6 Glieder seiner Reihe zu ersetzen. Gegeniiber dem Wert # = 0°462 fiir 
die Ablédsestelle nach Tani erhalten wir 
wieder eine etwas spatere Ablésung, bei 
x = 0°485 (prozentuale Abweichung zirka 
47%), 

In Abb. 3 haben wir den Verlauf der 
lokalen Wandschubspannung fiir alle von 
uns ausgewerteten Falle (n = 1, 2, 3, 4 
und 5) zusammen aufgetragen. (Fir die 
neue Reihe standen dabei 6, 3, 2, 2 und 
2 Glieder zur Verfiigung.) Dazu sind ein- 
getragen unsere Fortsetzungsergebnisse R, 
nach dem Differenzenverfahren und ferner 
fir n = 1 das Fortsetzungsergebnis nach 
Howarth, fiir » = 2 und m = 4 jene nach 
Tani. (Fir » = 1 standen Howarth 9 Glie- 
der, fiir nm = 2 und 4 standen Tani 7 bzw. 
6 Glieder der Reihe vom Blasiusschen Typus 
zur Verfiigung.) Als Ablosungsstellen fiir 
die zusatzlich berechneten Falle n = 1, 3 
Eo ee ee und 5 finden wir: Fir » = 1 Ablosung bei 

: x = 0°125 (gegentiber x = 0°120 bei Ho- 
warth, Abweichung 4:0%); fiir n = 3 Ab- 
losung bei x = 0°409; fiir n = 5 Ablosung 

sn eee eR ge EE oe 
Abb. 8. Die von der Vorderkante bis zur Der Verlauf der Verdrangungsdicken 
Ablésungsstelle iiberstrémte dimensionslose fiir die behandelten Falle ist in Abb. 4, 


Wandbogenlinge #4 in Abhingigkeit vom der Verlauf der Impulsverlustdicken in 
Exponenten n. Gestrichelt die praktisch nicht Abbe darcoscelle 
wesentlich abweichenden Ergebnisse von : 8 ; 


Howarth (n = 1) und Tani (n = 2, 4), ferner Abb. 6 gibt eine Auswahl der berech- 
das Ergebnis von Tani firm = 8. Es ist (vgl. neten Geschwindigkeitsprofile an verschie- 
(18)) #4 = %4 (u,/— 4,) In denen Wandstellen x wieder, und zwar 


fir n= 2, 3, 4 und 5 (fiir n=1 vel. 2). 

Um noch einen ins einzelne gehenden und bekanntlich empfindlichen Vergleich 
der Tanischen Ergebnisse mit den unsrigen zu bieten, ist in Abb. 7 fiir » = 2 und 
fir n = 4 die Geschwindigkeitsverteilung in der Grenzschicht in konstanten Wand- 
abstaénden entlang der Wand dargestellt. Dort, wo die Tanischen Ergebnisse von 
unseren abweichen, sind sie gestrichelt eingetragen. 

Sowohl in Abb. 6 wie in Abb. 7 sind die Ergebnisse der Fortsetzungsrechnungen 
dort zugrunde gelegt, wo die Reihenergebnisse unbrauchbar werden. 

Um schlieBlich die Abhangigkeit des dimensionslosen Wertes x der Abldésestelle 
von dem Exponenten n zu demonstrieren, sind die von der Vorderkante an bis zur 
Ablésung tiberstrémten Wandbogenlangen fiir n = 1, 2, 3, 4 und 5 aufgetragen. 
Es ist (vgl. (18)) % = %(— Up/Uy)"/n. 

(Eingegangen am 10. Januar 1957) 
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(Mitteilung aus dem Institut fiir Kinfithrung in die Stark- und Schwachstromtechnik) 


Uber den Kraftangriff des allgemeinen, nichtstationdren 
elektromagnetischen Feldes an der Materie 


Von H. Hofmann, Wien 


Zusammenfassung. Die in friiheren Arbeiten durchgefiihrten Untersuchungen iiber den 
Kraftangriff des stationiren elektromagnetischen Feldes an der Materie werden auf das allgemeine, 
nichtstationire Feld erweitert. Da die Maxwellschen Gleichungen prinzipiell zwei Anschauungen 
uber den Magnetismus der Materie, naimlich die Deutung durch Elementarstréme oder durch 
magnetische Mengen, zulassen, werden beide Auffassungen behandelt und die diesbeziiglichen 
Formeln angegeben. Es ergeben sich entsprechend dieser beiden Theorien auch verschiedene 
Ausdriicke ftir Energiestrom- und Impulsdichte des Feldes. 


I. Einleitung 


Die Behandlung des Kraftangriffes im station’ren elektromagnetischen Feld 
int} 2,3 brachte zunichst folgende Resultate: 

a) Es mu8 prinzipiell unterschieden werden zwischen der vom Feld auf einen 
Korper direkt ausgetibten Kraft — man spricht in diesem Fall auch von kérperfestem 
Kraftangriff — und der wegen der Anwesenheit des betreffenden Feldes auf den 
Korper insgesamt einwirkenden bewegenden, ponderomotorischen Kraft. Die direkte, 
kérperfeste Kraftdichte ist durch ein oder mehrere Produkte der Form ,,Mengendichte 
mal Feldstirke‘” gegeben. Die ponderomotorische Kraft beinhaltet sowohl den 
direkten, korperfesten Kraftangriff als auch die von der umgebenden Materie wegen 
Anwesenheit des Feldes auf die Oberflache des betreffenden Korpers ausgeiibte Kraft 
und 1a8t sich aus diesen beiden Teilkraften entsprechend Gl. (2) oder auch zum Beispiel 
aus der bei einer virtuellen Verschiebung umgesetzten Energie, also unter Verwendung 
des Energieprinzips, berechnen. 

Die k6rperfeste Kraft ist also allgemein 


% =| fav, (1) 
Vv 
({ direkte, korperfeste Kraftdichte), wahrend eben fiir die ponderomotorische Kraft 
Boma =|fAV + opddA=F+Hpad (2) 
4 A 


angesetzt werden kann. Da sich die ponderomotorische Kraft im elektromagnetischen 
Feld unter Einhaltung gewisser Bedingungen (s.1, Abschnitte I, IJ und VIII, 1) in 
Integrale iiber das Volumen oder die Oberflache des betreffenden Korpers allein 
umwandeln lat, erhalt man auch die anderen Formen 

Opond or | foona dV = ° Ypond dA. (3) 
Die in Gl. (3) auftretenden ponderomotorischen Volumskrafte {,,,q bzw. pondero- 
motorischen Oberflichenkrafte (Spannungen) },onq sind jedoch rein formaler Natur, 
da sie, wie der Vergleich der Gln. (2) und (3) zeigt, physikalisch deutbare Volums- 
und Oberflichenkrifte in sich vereinen. Sie haben vom Standpunkt einer molekularen 
Theorie der Materie aus gesehen nur Bedeutung als Rechengro8en, die bei Integration — 
in dem einen Fall iiber das Volumen, im anderen tiber die Oberflache des betrachteten 


1H. Hofmann: Uber den Kraftangriff des stationaéren elektromagnetischen Feldes an der 
Materie. Osterr. Ing.-Arch. 10, H. 4 (1956). : f 
2 H. Hofmann: Uber den Kraftangriff des Magnetfeldes an Elementarstrémen. Osterr. 
Ing.-Arch. 11, H. 1 (1957). ; a 
° H. Hofmann: Die Behandlung von Sprungflichen beim Kraftangriff des stationaren 
elektromagnetischen Feldes an der Materie. Osterr. Ing.-Arch. 11, H. 1 (1957). 
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Kérpers — allein bereits die bewegende Gesamtkraft liefern. Demgema6 ist zum 
Beispiel im elektromagnetischen Feld streng zwischen den aus dem Lorentz-Ansatz 
fiir das Makrofeld abgeleiteten direkten, korperfesten und den Maxwellschen pondero- 
motorischen Kraften zu unterscheiden (s.!, Abschn. IT). 

b) Setzt man als direkte, kérperfeste Kraftdichte unmittelbar nur das Produkt 
aus der makroskopisch resultierend in Erscheinung tretenden Mengendichte und der 
in Frage kommenden Feldstarke an (im statischen elektrischen Feld also , treie 
Ladungsdichte mal elektrische Feldstirke), so gelangt man bei K6rpern, die sich 
nicht im Vakuum befinden sondern von anderer Materie umgeben sind, durch Inte- 
gration einer solchen Kraftdichte allein noch nicht auf die direkte, kérperfeste Kraft 
auf einen abgeschlossenen Kérper. Man wendet in diesem Fall ,,mathematische” 
Bereichsabgrenzung an, die bei polarisierter Materie prinzipiell nicht ausschlieBlich 
ganze Materieteilchen erfaBt (s.4 und!, Abschn. III). An der Bereichsgrenze decken 
sich namlich die (makroskopisch flachenhaften) Polarisationsladungen oder Klementar- 
stréme benachbarter Materiebereiche und treten nicht in Erscheinung (ausgenommen 
bei Sprungflichen, von denen hier zunichst abgesehen wird). Die beiden gleich- 
groBen, sich deckenden Schichten verschiedenen Vorzeichens zahlen jedoch physikalisch 
gesehen zu verschiedenen Materiebereichen; die eine der beiden Schichten muf somit 
bei dem ins Auge gefaBten Korper zusitzlich zu der makroskopisch feststellbaren 
Mengendichte mitberiicksichtigt werden. Die Beziehungen, die diese andere ,,physi- 
kalische“’ Art der Bereichsabgrenzung zur Voraussetzung haben — sie erfaBt also 
prinzipiell nur ganze Materieteilchen (Dipole, Elementarstréme; s. wieder* und}, 
Abschn. III) — fiihren unmittelbar auf die direkte, korperfeste Kraft auf den be- 
treffenden Materiebereich. 

Aus dem Lorentz-Ansatz fiir die Kraftdichte oder aus dem entsprechenden Ansatz 
im Bereich der ,,Theorie der magnetischen Mengen“ gewinnt man nun durch direkten 
Ubergang auf die GréBen des Makrofeldes, also auf die gemittelten MikrogréBen, 
die Ausdriicke fiir ,,mathematische“ Bereichsabgrenzung, da die Mittelwertsbildung 
der MikrogroBen stets ,,.mathematische“ Abgrenzung zur Voraussetzung hat. Wenn 
in dieser Arbeit wieder die Kraftansatze fiir beide Arten der Begrenzung angegeben 
werden, so wird damit folgende Absicht verfolgt: Der direkte Anschlu8 einerseits 
an die Gleichungen des Mikrobereiches und anderseits an Literaturstellen, die sich 
zum Beispiel mit dem Impulssatz der Elektrodynamik unter Verwendung der GréBen 
des Makrofeldes befassen (s. etwa °, Abschn. 6), wird durch die mehr formalen Ansitze 
der ,,mathematischen® Abgrenzung gewonnen. Fiir die praktische Anwendung der 
direkten Kraftdichten hingegen kommen im Grunde genommen wieder nur die 
Ausdriicke fiir die ,,physikalische‘‘ Bereichsbegrenzung in Frage. Lediglich bei der 
Berechnung der Kraft auf Korper, die sich véllig im Vakuum befinden, spielt begreif- 
licherweise mangels angrenzender Materie der Unterschied der beiden Abgrenzungs- 
arten beim Integrationsergebnis keine Rolle, so daB auch die Ausdriicke fiir ,,mathe- 
matische‘ Begrenzung unmittelbar angewendet werden diirften. In diesem Falle 
miissen tibrigens auch die direkte und die ponderomotorische Kraft gleich sein (vgl. }, 
Abschn. VIII). 


II. Direkte, kérperfeste Volumskrafte im nichtstationairen elektromagnetischen Feld der 
» Theorie der Elementarstréme“. 


Im nichtstationéren Fall ist eine getrennte Betrachtung des elektrischen und 
des magnetischen Feldes nicht mehr zweckmafig, da die Bedeutung der jetzt auf- 
* H. Hofmann: Zur Frage der Bereichsabgrenzung elektrisch und magnetisch polarisierter 


Materie. Elektrotechn. u. Maschinenbau 73, H. 10 (1956). 
° F. Hund: Materie als Feld. Berlin-Géttingen-Heidelberg: Springer-Verlag. 1954. 
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tretenden zeitabhingigen Glieder der beiden Felder zum Teil nur gemeinsam erfaBt 
werden kann. Legt man den Untersuchungen die ,,Theorie der Elementarstrome“, 
also die Prinzipien der Elektronentheorie, zugrunde, so mu8 man vom Lorentz-Ansatz 
fiir die Kraftdichte im Mikrobereich Ae 


fo — oer (EB) stag fee, Ae. (4) 


[Index (E) Elektronentheorie bzw. ,,Theorie ao Elementarstrome“, 0) el. Ladungs- 
dichte im Mikrobereich, e™, §(® elektrische und magnetische Feldstiirke im Mikro- 
bereich nach der Elektronentheorie.] Ferner mu8 man beachten, daB man von 
den Grundgleichungen des Mikrobereiches auf umgeformte Maxwellsche Gleichun- 
gen fiir das elektromagnetische Makrofeld gelangt (s.*, Gl. I A bis I D, bzw. ®, Abschn. 4, 
Gl. 5, oder’, Gl. Ib bis IV b), wobei die Vektoren © und $ die Bedeutung der beiden 
Feldstarken des elektrischen und des magnetischen Teilfeldes erhalten. Sieht man 
die Materie als ideales Kontinuum an (kontinuierliche Verteilung wahrer Ladungen 
und von Dipolen bis zum (makroskopisch) unendlich kleinen Volumselement), so 
kann man auch bei der Kraftdichte Gl. (4) auf die GréBen des Makrofeldes iibergehen. 
Wegen der beiden Moglichkeiten der Bereichsabgrenzung werden wieder zwei ver- 
schiedene Ausdriicke fiir die Kraftdichte erhalten. 


1. ,Mathematische“ Bereichsabgrenzung 


Bei ,,mathematischer“* Bereichsabgrenzung kann man die fiir den Mikrobereich 
giltige Form der Gleichungen unmittelbar verwenden, das heiBt die entsprechenden 
GréBen des Makrobereiches direkt einfiihren: Gl. (4) fiihrt dann wegen 


di = of = oP + of = div D — div B = a6, 6) 
oy =g +e .rotM +e ee (6) 
e® = & (7) 
und ae 
5% = 8 (8) 


(en eGll5-7,'1) und 12). auf 
=o -E+— |g +cprotm+—F, 9}. (9) 


[Index (m) ,,mathematische* eee Dabei ist g die Leitungsstromdichte, 
der Ausdruck 

c,rot M = g (10) 
bedeutet die Dichte der Elementarstréme (Magnetisierungsstrome) und é%3/ot die 
Dichte der (elektrischen) Polarisationsstréme. Diese letzte Grof8e tritt zusatzlich 
zu den beiden ersteren, auch im stationadren Feld [s.1, Gl. (75)] wirksamen Stromdichten 
auf, da sich bei zeitlicher Anderung des Polarisationszustandes Elektrizitatsmengen, 
namlich Polarisationsladungen, bewegen. Mit Hilfe der urspriinglichen, besser aber 
noch der umgeformten Maxwellschen Gleichungen (s.’7, Gl. 1b bis IV b) erhalt man 
dann aus Gl. (9) 


G 
(® = Cdiv + [rob B 24,9] = EVE) +|VB1—~- 4.8] = 
= (7,66) — (EV) E + (BV) B—V(BB,) — + | 9]. (11) 


Colas Kneissler: Die Maxwellsche Theorie in verainderter Formulierung. Wien: Springer- 


Verlag. 1949. 
= a Hofmann: Uber die Deutung der Maxwellschen Gleichungen mit Hilfe elektrischer 


und magnetischer Mengen. Acta Physica Austriaca, 11, H. 2 (1957). 
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(Der Index c bedeutet, daB die GréBe beziiglich des vorherstehenden Differentiations- 
Legere konstant zu halten ist.) Wegen der IV. Maxwellschen Gleichung 


div 8 = VV 8) =0 (12) 

und der Umformungen 
V (BP = 2V (BB,), (13) 
(VI)=V (14) 


(I Einheitstensor) und 
(6 (V G]] =V (GG,) — (EV) E =>V (GF - GV)E, (15) 
ergibt sich aus Gl. (11) 


@ = (V7, 66 — $61) + CVE] + B7)8 +308) —FV @Fl-2[F3| - 
=(r,66—5 @1)+ (7, BB-F@rl+eEven- zs} 6) 


was unter Beriicksichtigung der II. Maxwellschen Gleichung schlieBlich auf 
1 1 We 
‘@ =(V,66 —5 @#1) + (v7, BB — 5 Bl) —— F IE] (17) 


fiihrt. Fiihrt man die in! definierten Tensoren des elektrischen Feldes und des 
magnetischen Feldes der ,, Theorie der Elementarstrome“ fiir den Fall ,, mathematischer“ 
Hoveichaaitereneung ein [1, Gl. (24) und (79)] 


TS =E6-5 @FI=TE, ie) 
1 ~ 
TiS (Ole (19) 
so ergibt sich weiterhin unter Bentitzung von 
Th” = TS +TS =€€ + BB => (6) + 1 = Th” ey 
die Beziehung 
a : ao™ 
i dv ie zs ss (ESpadiv-Tia = — (21) 
(vgl.*, Kap. IV), oder 
16 uae 
div 7 ee E) Sate | oe ee [EB] = tet é (22) 


Die Divergenz des Tensors T{*)™ liefert also die Summe aus der an den Ladungen 
des ,,mathematisch‘‘ begrenzten Bereiches angreifenden direkten Kraft des elektro- 
magnetischen Feldes und der zeitlichen Ableitung einer Gré8e, die als Impulsdichte G6 
des elektromagnetischen Makrofeldes im Sinne der Elektronentheorie gedeutet werden 
mu : 


1 
G® = ae [EB]. (23) 
Die Impulsdichte 6 und die mit ihr in Zusammenhang stehende Energiestromdichte 
CO == 67 OM =. 6, (G3) (24) 


sind dem Vektorprodukt der beiden Feldstirken des elektromagnetischen Feldes 
proportional, wobei im Rahmen dieser Theorie sinngemif das Vektorprodukt der 
GroBen © und B erhalten wurde. 

Es 148t sich leicht nachweisen, da die Anwendung des Impulssatzes der Elektro- 
dynamik, der direkt aus den Feldgleichungen abgeleitet wird, unmittelbar auf dasselbe 
Resultat fiir ,,mathematische“ Abgrenzung fiihrt, was zu erwarten ist, da ja die 
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Feldgleichungen des Makrobereiches nur die makroskopisch in Erscheinung tretenden 
Mittelwerte der Ladungs- und Stromdichten enthalten und somit ,mathematische“ 
Begrenzung der Volumselemente bei der Mittelwertsbildung der GréBen des Mikro- 
bereiches zur Voraussetzung haben. Man geht am besten gleich von den im Sinne 
der Elektronentheorie umgeformten Maxwellschen Gleichungen (?, Gl. Ib bis IV b) 
aus und folgert aus ihnen (s. auch *®, Abschn. 6): 


[rot 8, B] + [rot , €] = [V B] B] + [VE] E] = 
= (BV) 8 —7(B%.) + (EV)E—V (6G) — + 


Co 


£68] +2 g +eorot m +, 8]. 


Summiert man hiezu (25) 
(VG) = 9’ - G, (26) 


so erhilt man mit den Umformungen Gl. (13) und (14), den Definitionen (18), (19) 
und (20) und der IV. Maxwellschen Gleichung (12) wieder 


(7, BB — 5 (BPI) +(V, CE —F ©1) = (V Thy”) = div TS” = 
' 1 l 
=o - E+ |g +cprot ® +>, B] + — | (EB) =(B +2 Z1EB) (27) 


in Ubereinstimmung mit Gl. (22). Fiir theoretische Betrachtungen also, die die 
Klarung allgemeiner Prinzipien aber keine weiteren praktischen Auswertungen zum 
Ziele haben, wird man somit mit Vorteil die Ausdriicke fiir ,,mathematische‘‘ Bereichs- 
abgrenzung heranziehen kénnen. 


2. ,Physikalische“ Bereichsabgrenzung 


Bei der Ableitung der Kraft auf einen ,,physikalisch‘‘ begrenzten Bereich gibt es 
wieder die beiden Wege: Entweder es werden zur Erganzung auf ganze Dipole zu der 
Kraft auf den ,,mathematisch“ begrenzten Bereich, die z. B. mit Hilfe der Kraft- 
dichte Gl. (9) oder Gl. (21) angeschrieben werden kann, die Krafte auf die in der 
Bereichsgrenze latent vorhandenen flachenhaften Dichten der elektrischen Polari- 
sationsladungen und der Elementarstro6me hinzugezahlt: 


1 
BEY = \{o'?- E+ = |g +eorotm +P, Bll av +-h{o— -€+— jm Bilas = 
V A 
Co 


= | {div TS” — = [EB AV + O{(By E+ [Mn] Bh aA. (28) 
¥ Vig 


(Index (ph) ,,physikalische‘‘ Abgrenzung.) In Gl. (28) fanden dabei die Beziehungen 
fiir die Flachendichte der elektrischen Polarisationsladungen [1, Gl. (16)] 


op = (Pn) (29) 
und die Flachendichte der Magnetisierungs- (Elementar-) Strome [*, Gl. (17)] 
j™ = eg [Mn] (30) 


Anwendung. Oder aber man geht direkt von der Kraftdichte fiir ,,physikalisch‘‘ 
begrenzte Materie aus, wobei zu den Ausdriicken des stationaéren, von elektrischen 
Ladungen und Elementarstrémen erzeugten elektromagnetischen Feldes [', Gl. (29), 
und 2, Gl. (11)] im allgemeinen, nichtstationiren Feld die ausschlieBlich raéumlich 
verteilte Stromdichte des Polarisationsstromes 0%3/ot hinzutritt: 

@, — of .G + (B grad) € + — [gB] + grad (MB) +—-[G 8] BD 


Q* 
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Beide Wege fiihren natiirlich auf dieselbe Kraftdichte, die jetzt aus der urspriing- 
licheren Gl. (31) entwickelt werden soll: Mit Hilfe der Maxwellschen Gleichungen 
folgt aus ihr 


fiom = EVD) + (BV)E + jrot § — 


Diese Beziehung wird weiter umgeformt, wobei die aus der II. Maxwellschen Gleichung 
gewonnene Umrechnung 


B] +7 (MB) + 2-| GB] +508). (32) 


mle 
oe ee 


—|¢- 4] = rot ] = (CV EN =" CE) — (EV) E= ZI GP—(CV)E (33) 
und die Definitionsgleichungen 
d-B+E (34) 
und 
B= WH (35) 


Verwendung finden: 
o€ 


EVD) + (D—-GV)E + BI) H —V (HB) +7 MB) — 7 [FB] + 97 B) = 
— (V,D6) — (EV) E+ (VV, BH) — 78 COC ele 
=(V,D6)— 37 © —— [EZ] +, BH) —FV (BP +7 MB) —— || Bl = 
Deer arty iether) (36) 


Mit Hilfe der Beziehung (14) und der in! [GIl. (32) und (86)] definierten Tensoren fiir 
» physikalische“ Abgrenzung, 


T2, =ED—> (GPI = CD—4 (DG) ees (BE) I (37) 
S Reece A EA TERA hI e = §B— + (BH) 1+ > (MB) T (38) 


die fiir das elektromagnetische Gesamtfeld der elektrischen Ladungen und Elementar- 
str6me zu dem Tensor 
Baye Se) py al eet © Es GBF {( (DE) — (VE) + (VH)—(MB)} I (39) 


zusammengefaBt werden konnen, folgen wieder die den Gl. (21) und (22) fiir ,,mathe- 
matische‘ Abgrenzung analogen Beziehungen . 


ERE Pe oft Ri a oR Bie cia Sie hy pes 
phy = div Tony” — GG [EB] = div Ton” — a (40) 
und 
div ie E)o = f@, + = (6 ¥]= en fe ae a : (41) 
(ph) (phy + (ph) 


III. Direkte, kérperfeste Volumskrafte im nichtstationiren elektromagnetischen Feld 
der ,,Theorie magnetischer Mengen‘ 


Die Maxwellschen Gleichungen lassen sich nicht nur im Sinne der Hlektronenthewee 
deuten, sondern kénnen auch durch Mittelwertsbildung aus den fiir den Mikro- 
bereich giiltigen Grundgleichungen einer ,,Theorie (elektrischer und) magnetischer 
Mengen“ gefolgert werden (s.7, auch1, Abschn. V). Wahrend nach der Elektronen- 
theorie neben der Leitungsstromdichte g die Dichte der Elementarstréme g® die 
Bedeutung der im Magnetfeld in Frage kommenden Mengendichte besitzt und als 
Feldstairke des magnetischen Makrofeldes der Vektor 8 anzusehen ist, tritt im Bereiche 
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der ,,Theorie magnetischer Mengen‘ neben der Leitungsstromdichte g die Dichte 
der ,,freien“’ magnetischen Mengen 


o/men) — div § te div Ih = omen) (42) 


als Mengendichte auf, als magnetische Feldstiirke ist hier der Vektor § anzusprechen. 
Die sinnwidrige Verkniipfung von Groen der beiden Theorien kann natiirlich keines- 
falls auf physikalisch richtige Begriffe fiihren, so muB z. B. eine von Sommerfeld 
und Bopp in® vorgeschlagene ,,Kraftdichte“ des Magnetfeldes, die durch Produkt- 
bildung der Dichte der freien Magnetmengen 0’ und der Feldstarke der ,,Theorie 
der Elementarstrome“ 8 entsteht, abgelehnt werden. Die GréBe B oder, wie an anderer 
Stelle® (§ 28, B) behauptet wird, im Falle ec mit der Geschwindigkeit » bewegten 


magnetischen Ladung der Ausdruck 8 — — * Tv €], kann nicht zur Berechnung einer 


Kraft auf einen magnetischen Einheitspol there cozoeen werden, da 8 ja nur dann 
die Bedeutung der Feldstiirke des Magnetfeldes erhailt, wenn die Maxwellschen 
Gleichungen im Sinne der Elektronentheorie ausgelegt werden, in deren Rahmen es 
prinzipiell keine Magnetmengen gibt. Setzt man jedoch die Existenz von magnetischen 
Mengen voraus, die, solange die IV. Maxwellsche Gleichung Giltigkeit besitzt, ent- 
sprechend Gl. (42) nur in Form von Polarisationsmengen auftreten, so ist, wie in? 
und auch schon in ! (Abschn. V) gezeigt wurde, als magnetische Feldstarke der Vektor § 
zu verwenden [s. die folgende Gl. (50)]. Die physikalisch sinngemife Anwendung 
der Vektoren § und & stellt auch die Voraussetzung dafiir dar, da man die direkten 
korperfesten Kraftdichten im stationiren elektromagnetischen Feld als Divergenzen 
einander analoger Tensoren darzustellen vermag. Ebenso nur unter dieser Voraus- 
setzung kann man die Ubereinstimmung der aus dem Energieprinzip gewonnenen 
Maxwellschen bewegenden, ponderomotorischen Krafte mit den bewegenden Kraften 
nachweisen, die man mit Hilfe der aus der Elektronentheorie (Lorentz-Ansatz) oder 
der ,,Theorie der magnetischen Mengen“ abgeleiteten direkten Feldkrafte erhalt 
(s. 1, Abschn. VIII). 


Zu den Grundgleichungen der ,,Theorie der elektrischen und magnetischen Mengen‘‘ 
(7, Abschn. 2) tritt also noch der Ansatz fiir die Kraftdichte im Mikrobereich 


ND) isn oe +> [ocr A 7] + ene) = ee v, a} (43) 


[Index (M) ,,Theorie der elektrischen und magnetischen Mengen“, 0), 0") el. 
bzw. magn. Ladungsdichte im Mikrobereich, e™, ) el. bzw. magn. Feldstarke im 
Mikrobereich nach der ,,Theorie elektrischer und magnetischer Mengen‘‘]. Die beiden 
letzten Summanden in Gl. (43) stellen das magnetische Analogon zum Kraftansatz 
auf elektrische Ladungen dar, dabei wird das negative Vorzeichen im letzten Summan- 
den (der ,,magnetischen Lorentz-Kraft‘) von grundlegenden Prinzipien, wie etwa 
dem Rélativitateprinzip, gefordert; es hangt letzten Endes mit dem negativen Vor- 
zeichen in der II. Maxwellschen Gleichung zusammen. (Uber die Invarianz der 
Grundgleichungen der ,,Theorie der elektrischen und magnetischen Mengen“ gegen- 
iiber einer Lorentz-Transformation wird noch an anderer Stelle berichtet werden.) 


8 Sommerfeld und Bopp: Zum Problem der Maxwellschen Spannungen. Ann. Physik 8, 
H. 1/2 (1950). 


9 Sommerfeld: Vorlesungen iiber theoretische Physik, Bd. III, Elektrodynamik. Leipzig: 
Akad. Verlagsgesellschaft Geest und Portig. 1949. Die Beziehung wird von Sommerfeld wegen 


1 
der Verwendung anderer GréBendefinitionen in der Form 8 — ai [v €] angegeben. 
0 
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1. ,Mathematische“ Bereichsabgrenzung 


Wie schon mehrfach erwihnt, gestattet die ,,mathematische“ Bereichsabgrenzung 
einen direkten Ubergang von den Gleichungen des Mikrobereiches zu denen des Makro- 
bereiches, wenn man die Materie als ideales Kontinuum auffa8t. Es konnen in diesem 
Fall auch beim Ansatz fiir die Kraftdichte die MakrogréBen, also die gemittelten 
GréBen des Mikrobereiches, eingefiihrt werden, wobei im Bereich der _,,Theorie der 
elektrischen und magnetischen Mengen“ die folgenden Beziehungen zu gelten haben: 


pe) = op’) = ped + ge) — div D — div f = divG, (44) 
gimen) — 9/(mgn) — o(men) — _ div NQ = div §, (45) 
oar Wa 4) 
eu =g +S, (46) 
etmen) y = ae (47) 
em — G, (48) 
= §, (49) 
(s. 7, Gl. 16, 20, 22, 25, 27 und 28). Man erhalt dann aus Gl. (43) 
, ; é 
Mao E+ — [ote S + — lai e e. (50) 
0 


Zu den im stationiaren Feld giiltigen Ausdriicken [!, Gl. (18) und (61)] kommt also 
jetzt im allgemeinen Fall noch der Kraftangriff auf die Dichten des elektrischen und 
des magnetischen Polarisationsstromes o8/ot und ot/ot. Unter Verwendung der 
Maxwellschen Gleichungen und der Beziehungen (34), (35), (14), (15), (44) und (45) 
folgt weiter aus Gl. (50) 

oD 


iv EV 6) +H) +|5I-— 4 9) +-[F 9]-s 
= V, 66) — (EV) E + V7, HH) — (GY) H + (HV) § — FV (8 -2]F 9]— 

a §|— (V7, €€—5 1) + CVE] +(7, 99 — | (HP 1) — 

—~ 2% s|+ ze F=(". ce -$@ 1) +(7,99— Foy) - 

— > @ (E91. (51) 


Beniitzt man wieder die in! [Gl. (24) und (66)] definierten Tensoren fiir ,,mathe- 
matische‘‘ Bereichsabgrenzung 


om 
at 


g|= 


Te} =€E—+ (GP 1= Te, (52) 
T= $$ —F (6PI=T%, (53) 


bzw. deren Summe 
TEs? = TS + TR = CE + 9H — 51? + (H}1= TA”, (54) 
so folgt aus Gl. (51) 
e. g : e. 
(din ay ee cs ene [(EG)= div Ts 


EGO 
ee et (55) 
oder 


- (el, M) __ ¢(M a@ 
div Tm) im) + 5% | ($1 = (08 + ays (56) 
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Fir Impuls- und Energiestromdichte erhailt man also sinngemé8 im Rahmen der 
Theorie elektrischer und magnetischer Mengen“ die aus den beiden Feldstarken & 
und § dieser Theorie aufgebauten Ausdriicke 


G% = — [EG] (57) 


SOO = 04? GOO = 0, [CH]. (58) 
Die Gl. (55) und (56) besitzen wieder dieselbe Gestalt wie die analogen Beziehungen (21) 
und (22) bzw. (40) und (41) der ,,Theorie der Elementarstréme‘‘. 

Ks ist zu erwarten, da man auch im Rahmen der ,,Theorie elektrischer und 
magnetischer Mengen“ einen Impulssatz der Elektrodynamik ableiten kann, der 
auf das Resultat (55) bzw. (56) fiir ,mathematische“ Abgrenzung fiihrt. Die Feld- 
gleichungen, aus denen der Impulssatz direkt gewonnen wird, setzen ja wieder bei 
ihrer Ableitung durch Mittelwertsbildung der Gleichungen des Mikrobereiches ,,mathe- 
matische‘‘ Begrenzung der Volumselemente voraus. Der Impulssatz folgt jetzt aus 
den Maxwellschen Gleichungen, die im Sinne der ,,Theorie elektrischer und magneti- 
scher Mengen‘‘ umgeformt sind (?, Gl. Ic bis IV c). Man erhalt aus ihnen: 


[rot 6, €] + [rot §, 1 = — =| | - |e] +< to 91+—[ 4 9] + 


Ot Ot 
+z|4 |= (7) 6-37 GP +9”) H- FV (GP=+|F 9+ 


+= [981-< |G 6] +<-% (E91. (59) 


und 


Summiert man hiezu 
EVE) +HVH)=o.C + omg, (60) 
so ergibt sich mit. Gl. (14) 
(V, 6 — > (G1) +(V, 8 — F191) = VTS) + VTS) = TS”) = 


Se eC ee 


- 7H ~ [ES] =18 +4 =) ~ (66). (61) 


Der direkt aus den Feldgleichungen des Makrobereiches abgeleitete Impulssatz fiihrt 
also im Rahmen der ,,Theorie elektrischer und magnetischer Mengen“ ebenfalls auf 
die entsprechende, aus dem Kraftansatz Gl. (43) gewonnene Beziehung Gl. (56). 


2. Physikalische“ Bereichsabgrenzung 


Summiert man zu den Kraftdichten des stationiren elektromagnetischen Feldes 
fiir ,,physikalische‘‘ Abgrenzung [!, Gl. (29) und (69)] die im allgemeinen Feld zu 
beriicksichtigenden Kraftdichten auf die elektrischen und magnetischen Polari- 
sationsstréme, so erhadlt man 


1 1 [ om 
(29, = eva € + (Berad) E + [FP 9] +2 [9 $1 + (Merad) § — 3-| FG] = 


=€VD)+ (DV) E—(EN)E- 2] F cal 


+ (7 $191-3.|G S|=.20-F eo, aacneai 


—}rir Lf F s|-HfBg-("-26—4 01) 


+ (7,85 —(9"1) -—<. 5 (ES). (62) 
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Die Definitionen der Tensoren fir ,,physikalische‘ Bereichsabgrenzung des elektri- 
a Feldes und des magnetischen Feldes der ,,Theorie magnetischer Mengen“ 


[', Gl. (32) und (72)] 


T?,=ED—+4 GP1=ED-—4FMGH1+ 5 (BH! (63) 
Tim = 8B Liga ies eee ee (64) 
und der daraus gebildete Summentensor des Gesamtfeldes 
Ten = Oa DO 1 62+ (HP}1 (65) 
erméglichen wieder, zusammen mit Gl. (57), die schon bekannten Schreibweisen 
F) ei oe. ag” 
(08, = div Tei 2 2 eg) — aiv THY — (66) 
und 
é ao 
div Ti)” = fom + a (ES) = fom + —Q- (67) 


IV. Zusammenfassung der Ergebnisse 


Mit der Darstellung der korperfesten Krafte des allgemeinen elektromagnetischen 
Feldes ist ein gewisser AbschluB bei der Untersuchung des Kraftangriffes dieses Feldes 
an der Materie erreicht. Es zeigte sich, daB die direkte, korperfeste Kraftdichte die 


Form 
fe = div T,— 8 (68) 
besitzt [Gl]. (21), (40), (55), (66)], wobei je nach te verwendeten Theorie, und inner- 
halb einer solchen wieder entsprechend den Begrenzungsarten der Materie, ver- 
schiedene Tensoren T; [GIl. (20), (39), (54) und (65)] und Impulsdichten ©, [Gl. (23) 
und (57)] zur Anwendung kommen. Die einzelnen Tensoren und Impulsdichten 
besitzen einander vollig analoge Gestalt, wobei beziiglich der Tensoren insbesondere 
auf 1(Abschn. VII) hingewiesen wird; die Impulsdichten sind dem vektoriellen Produkt 
der jeweiligen elektrischen und magnetischen Feldstirke der verwendeten Theorie 
proportional. 

Leitet man den Impulssatz im Rahmen der klassischen Maxwell-Theorie ab, so 
ergibt sich ebenfalls ein Ausdruck der Form Gl. (68), [s.1, Gl. (1)], der im Rahmen 
dieser Theorie als ponderomotorische juntos gedeutet wird: 


lnoadc= a div T— = Cy * 1D9] (69) 


(T = klassischer Maxwell-Tensor). Im stationaren Feld fallt der zeitabhangige zweite 
Summand in den Gln. (68) und (69) weg, die direkte und die ponderomotorische Kraft- 
dichte lassen sich allein durch die Divergenz eines Tensors ausdriicken. In diesem 
Fall 1aBt sich der Nachweis erbringen, da8 man von den direkten Kraften ausgehend 
zu denselben bewegenden Kraften gelangt, die durch die Divergenz des klassischen 
Maxwellschen Tensors gegeben sind, wenn man dabei auf die Bedingungen Bedacht 
nimmt, unter denen die Maxwellschen ponderomotorischen Krafte angewendet werden 
diirfen (s.1, Abschn. VIII). 

Ks erscheint aber aussichtslos eine derartige Ubereinstimmung auch im allge- 
meinen Feld nachweisen und damit zeigen zu kénnen, dali die Maxwellsche Hypothese, 
die ponderomotorischen Krafte Gl. (69) behalten auch im nichtstationiren Feld ihre 
Giltigkeit, den Tatsachen entspricht. Vom Standpunkt der beiden molekularen 
Theorien der Materie, der Elektronentheorie (,,Theorie der Elementarstréme“‘) und 
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der ,,Theorie elektrischer und magnetischer Mengen“, aus gesehen enthilt die zeitliche 
Anderung der Maxwellschen Impulsdichte, also der zweite Summand in Gl. (69), 
zum Teil noch Ausdriicke, die einen Kraftangriff an der Materie liefern, namlich 
die Kraftdichte auf die elektrischen und allenfalls auch auf die magnetischen Polari- 
sationsstréme. Die tatsichliche bewegende, ponderomotorische Kraft im allgemeinen 
Feld kann, auch im Rahmen der hier behandelten Theorien der kontinuierlich aus 
Ladungen und Doppelschichten (Elementarstrome) aufgebauten Materie, schon aus 
folgenden Griinden nicht die einfache Gestalt Gl. (69) besitzen: Die ponderomotorische 
Kraft besteht bekanntlich aus der im Inneren des Korpers durch das Feld hervor- 
gerufenen direkten, kérperfesten und der im jeweiligen Augenblick an der Oberflaiche 
durch die umgebende Materie angreifenden Kraft, es gilt im Falle eines starren Kérpers 
auch im nichtstationiiren Feld die Gl. (2) 


Sona = | FAV + Opa. (70) 
Vv A 

Der direkte Kraftangriff an jedem Volumselement ruft im zeitlich veranderlichen 
Feld, selbst wenn der betrachtete Kérper als starr angenommen werden kann, zumin- 
dest in der umgebenden Materie nichtstationare elastische Spannungen und Defor- 
mationen hervor, die sich mit Schallgeschwindigkeit fortpflanzen. (Im stationiren 
Feld geniigte begreiflicherweise zur Erfassung der Kriafte der umgebenden Materie 
die Gleichgewichtsbedingung zwischen elektromagnetischen und mechanischen Volums- 
kraften.) Zur Berechnung der bewegenden Kraft miiBten also gemaB Gl. (70) auBer 
der direkten Kraft auf den betreffenden starren Kérper auch die im betrachteten 
Moment auf seine Oberfliche wirkenden Druckkrafte bekannt sein. Diese miiBten 
ihrerseits aus der direkten Kraftdichte auf die umgebende Materie unter Beriick- 
sichtigung aller Rand- und Anfangswerte und unter Verwendung der in Frage kommen- 
den mechanischen Materialkonstanten berechnet werden. Es erscheint daher schon 
auf Grund der Tatsache, daB die ponderomotorischen Maxwell-Krafte Gl. (69) keinerlei 
Bezug auf mechanische Materialkonstanten enthalten, unmoglich, daB sie die richtigen 
bewegenden Krafte auf einen Kérper liefern, wenn sich dieser nicht im Vakuum be- 
findet. An dieser Stelle sei besonders darauf hingewiesen, daB die direkten Volums- 
krafte zwar allgemeine Giiltigkeit im Rahmen der betreffenden Theorie besitzen, 
im nichtstationjren Feld jedoch ihre Integration im ersten Summanden der Gl. (70) 
nur bei einem starren K6rper sinnvoll sein kann. 

Es bleibt also noch der einfachere Fall des vollig vom Vakuum umgebenen starren 
Korpers zu untersuchen, Direkte (kérperfeste) und ponderomotorische Kraft miissen 
jetzt mangels an Kraften angrenzender Materie gleich sein; auch darf kein Unter- 
schied in den Rechenresultaten bestehen, die mit Hilfe der Ausdriicke fiir ,,mathe- 
matische“ und ,,physikalische‘‘ Abgrenzung gewonnen werden. 

Zunichst soll die Kraft auf den starren Korper mit Hilfe des GauBschen Satzes 
umgeformt werden: 

Si. = | div T.av — | SE dV = (Te ae) — | dvs (71) 

V 4 A Vv 

Das erste, umgeformte Integral gibt jetzt mit Hilfe jedes beliebigen Tensors, also 
auch des klassischen Maxwell-Tensors, dasselbe Resultat, da an der Oberfliche des 
Korpers Dy = E), By = Ho, Bo = 0 und M, = O gilt. Der zweite Summand in Gl. (71) 
liefert jedoch im Rahmen jeder Theorie einen anderen Wert, da die Ausdriicke fur 
die Impulsdichten verschieden sind. Diese sind zwar von der Begrenzungsart unab- 
hangig, so daB im Rahmen der ,,Theorie der Elementarstro6me“ oder der. ,,Theorie 
elektrischer und magnetischer Mengen“ beide Begrenzungsarten erwartungsgemaB 
auf dasselbe Resultat fiihren ; die direkte und in diesem Fall mit der ponderomotorischen 
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identische Kraft des nichtstationiren Feldes dieser beiden Theorien ist aber ver- 

schieden und stimmt schlieBlich auch mit der der Maxwellschen Theorie nicht tiberein. 
Die neuere Literatur lehnt in der iiberwiegenden Mehrheit die Maxwellschen Krafte 

des nichtstationaren Feldes, insbesondere wegen der Impulsdichte dieser Theorie 


6 = = [DB], (72) 


ab (vgl. 9, § 35). Es wird praktisch nur die auf dem Lorentz-Ansatz aufbauende Kraft- 
dichte der ,,Theorie der Elementarstrome“ (und zwar fiir ,,mathematische* Bereichs- 
abgrenzung) anlaBlich mehr formeller theoretischer Betrachtungen (Impulssatz, 
vierdimensionale Schreibweise u. 4.) verwendet. Dieser Kraftdichte ist aber vorlaufig 
der Ausdruck der ,,Theorie elektrischer und magnetischer Mengen“ als ebenfalls be- 
rechtigt gegeniiberzustellen, da einerseits die Feldgleichungen dieser beiden Theorien 
im Makrobereich dieselben sind und mit den Maxwellschen Gleichungen tberein- 
stimmen; anderseits liefern beide Theorien im stationéren Feld in den tiberpriifbaren 
Fallen gleiche bewegende Krafte, die ebenfalls mit denen der Maxwell-Theorie in Uber- 
einstimmung sind. Die Differenzen der drei Theorien hinsichtlich des Kraftangriffes 
im nichtstationiren Feld bediirfen dagegen noch der Klirung, die allerdings auf 
direktem experimentellem Wege schwer durchfiihrbar erscheint. 


(Hingegangen am 5. Februar 1957) 


Beitrag zur nichtlinearen Torsion 
Von F. Jindra, Stuttgart 
Mit 4 Textabbildungen 


Zusammenfassung. Mit einem nichtlinearen Elastizitaétsgesetz fur kleine Verzerrungen werden 
die Losungsansaétze fiir die reine Torsion prismatischer Stabe angegeben. Unter Beschrankung 
auf einfache Formen der Schubfunktion wird die Weiterbehandlung der aufgestellten Differentiai- 
gleichung erméglicht. Die Losungen betreffen die Querschnittsformen: Kreis, Ellipse, gleichseitiges 
Dreieck und Quadrat. Die Zahlenrechnungen ergeben, da sich die zusatzlichen Spannungen 
infolge der Nichtlinearitat bei diesen Querschnitten etwa in gleichen Grenzen halten. 


I. Einleitung 


Die Aufgabe der Elastomechanik ist die Ermittlung der inneren Krafte und der 
Formanderungen eines in bestimmter Weise belasteten und gelagerten elastischen 
Korpers. Zur Loésung verfiigen wir tiber die Gleichgewichtsbedingungen, den geo- 
metrischen Zusammenhang zwischen den VerzerrungsgréBen und den Ableitungen 
der Verschiebungen und das Elastizitatsgesetz. Die klassische Elastizitatstheorie 
beruht auf dem Hookeschen Gesetz, das einen linearen Zusammenhang zwischen 
Spannungen und Verzerrungen voraussetzt. Dabei wird die Giiltigkeit des Gesetzes 
auf unendlich kleine Verzerrungen beschrankt. Manche zihe metallische Werkstoffe, 
wie Kupfer, Bronze, ziher Stahl, zeigen aber schon bei diesen kleinen Verzerrungen 
ein nichtlineares Verhalten. Hier beniitzen wir deshalb ein nichtlineares Elastizitits- 
gesetz fiir kleine Verzerrungen, das dem elastischen Verhalten solcher Werkstoffe 
gerecht wird. 

In rechtwinkligen Koordinaten x, y, z sind die von den Schubspannungen hervor- 
gerufenen Winkelénderungen nach dem allgemeinen nichtlinearen Elastizitatsgesetz 
fiir kleine Verzerrungen! gegeben durch die Ausdriicke: 


1 1 1 
Wey = E Ito”) Tx ys Py 2 = “G g(t) Ty 2» Vea = @ g(t”) Tz a (1) 


1H. Kauderer: Ing.-Arch. 17, 450 (1949). — C. B. Biezeno und R. Grammel: Technisch 
Dynamik, 2. Aufl., Bd. I, 8.35. Berlin. 1953. eine te 
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Dabei bedeuten 1,,,, t,,, T die Schubspannungen, Yow Yux Yee die Winkelanderungen, 
G den Schubmodul, ¢,2 die dimensionslose SpannungsgrdBe 
2 61 
ty e Ce 2 (Fe Oo)” (oy Oo)” fe (o, 0o)*] can ea =| noe (2) 


Weiter bedeuten hier o,, o,, o, die Normalspannungen, ¢,, ¢,, ¢, die Dehnungen und o, 


i 1 
den Mittelwert o, = 3 (9, +o, + ,). Die Schubfunktion g(é,2) wird als Potenzreihe 
der SpannungsgréBe (2) angesetzt: 


(to?) = 1 + gate? + Gate’ + Goto +... (3) 
Die Gl. (1) lassen sich umkehren. Man erhilt 
Tay = G v(o") Puy Ty, = G v(Wor) YPuz Ta =G ¥ (Wo) Dz x» (4) 
wo die DehnungsgréBe y,? durch 
2 
Yo = 3 12 [ex — eo)? + (&y — &)? + (& — &)*] + Pos? + Yur? + Yrart (5) 


mit 


0 = (es + &y + &) 


gegeben ist. Die Scherungsfunktion y(y,2) nimmt als Potenzreihe von y,? die Ge- 
stalt an: 


V(Yo) =1t+ yey + yao + repo t+.--- (6) 
Nach den Regeln tiber die Umkehrung von Potenzreihen gilt fiir die Koeffizienten 
6 
Nes Gu Peeler Gar hs 302 Oa 
Ye = — 129.7 + 89294 — Go (7) 


Ya — 559.4 — 559,79, + 59,7 + 109.95 — gg, USW. 


Die Schubfunktion (3) bzw. die Scherungsfunktion (6) bestimmen also mit G zu- 
sammen das elastische Verhalten des Werkstoffes. Ihr Verlauf ist experimentell 
durch geeignete Versuche zu ermitteln. 

Die vorliegende Untersuchung beschaftigt sich mit den Anwendungen des nicht- 
linearen Elastizitaétsgesetzes fiir kleine Verzerrungen auf das Problem der reinen 
Torsion prismatischer Stabe. Wir leiten zuerst allgemein die Differentialgleichung 
fiir die Verwolbung und die Spannungsfunktion her. Wir betrachten kurz die exakte 
Losung fir den Kreisquerschnitt. Das elastische Verhalten von zihem FluBstahl 
bei einem Torsionsversuch wird als Beispiel angegeben. Unter Beschrankung auf 
bestimmte einfache Formen der Schubfunktion (3) behandeln wir den elliptischen 
und den dreieckférmigen Querschnitt. Zum SchluB berechnen wir naherungsweise 
die maximalen Schubspannungen im quadratischen Querschnitt. Die numerische 
Auswertung der Ergebnisse ergibt Anderungen der Spannungsverteilung gegeniiber 
der linearen Theorie und insbesondere einen Abbau der Spannungsspitzen durch 
Abweichungen des Elastizitatsgesetzes von der Linearitat. 


II. Die reine Torsion prismatischer Stabe 
Wir behandeln hier die Torsion eines zur z-Achse parallelen prismatischen Stabes 

von beliebigem Querschnitt unter einem Verdrillungsmoment M. Fir das reine 
Torsionsproblem? sind bekanntlich nur die Schubspannungen 1,, und 1,, beizu- 
behalten und alle anderen Spannungen gleich Null zu setzen, also 

Cg 10) Gy =F 4 == 0, 

Ty 2 = Ty(%, Y); Tae = Ta(L; Y). } 
2 Vgl. z. B. S. Timoshenko: Theory of Elasticity, 1. Aufl., 8.229. New York. 1934. 


(8) 
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Dann bleibt als einzige Gleichgewichtsbedingung tbrig 
Chik Ty 


wobei die Volumkrafte vernachlassigt werden. Die Randbedingung ist die Forderung 
tangentialen Verlaufs der Schubspannungen. Wenn dz, dy die Projektionen des 
Bogenelements ds der Randkurve des Stabquerschnitts sind, so mu8 langs des Randes 
gelten “eal dy ty de 


dy 
gee oder T, 7 eae i a 5s 0. (10) 


Die Verschiebungen nach den Achsenrichtungen seien u, v, w. Wegen (8) werden 
die Verzerrungen 


ou ov aw ou Ov 
ame at ee 


0 0 
: a ow ou ow i) 
Vis og Gy Pra ar aes 
woraus die Vertraglichkeitsbedingung des Verzerrungszustandes 
Wee Wve _ Konst, (12) 


oy 0x 
folgt. Fiir die Aufstellung der Gleichungen braucht man noch den Zusammenhang 
zwischen den Spannungs- und den Verzerrungsgr6Ben. Das nichtlineare Elastizitats- 
gesetz liefert nach (1) bzw. (4) die Beziehungen 


1 1 ‘ 
Yuz = @ Ite) Ty) Yaw = @ G(to’) Tx (13) 
mit 
2 
to’ = -3 qm (t.” + 7”) (13a) 
oder 
Fy = EB y(Ye) Yun Te =F VME) Pre (14) 
mit 
2 
Yor = 3 (v2 ear) (14a) 


Fir das Integrationsproblem ergeben sich folgende zwei Fassungen: 

a) Die Differentialgleichung ftir die Verw6lbung. Die Verdrehung um 
die z-Achse fiir die Langeneinheit des Stabes sei w. Wir machen den iiblichen Ansatz 
der linearen Theorie 

U=— OY 2,  V= 082, “w= w Oa, ¥), (15) 
wo 9(z, y) die Verwélbung des Querschnitts ist, die fiir alle Querschnitte gleich ausfallt. 

Damit wird nach (11) und (14) 


t2 =G yy?) (Zo — 9), 


ty = Ay (yi) (32 + 2) wh 


2 ap \2 / ap \2 a @ 

ye = 3 0 (32) et) re ay 2y a, +a 4 |. (17) 

Fihren wir diese Werte in (9) ein, so entsteht die folgende nichtlineare Differential- 
gleichung zweiter Ordnung fiir die Verwélbung 9(, y): 


Op ay 4, Ayo") |/ Oy 2 &y op é o 
rive) (aor + ar) +5 ay [Cae 8) aoe +2(38 "(ae +2) Grey + 


-k ee aie x) ia a=): (18). 
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Die Randbedingung (10) nimmt fiir g die Form an: 


a) 0 1 
=e dy dx => d(x? + y?). (19) 


Die GréBe des Torsionsmomentes M erhalten wir durch Summation der Momente 
der Schubspannungen um z iiber dem Querschnitt: 


M=\\ (t,x — t, y) dx dy a 
a rh op op 
= Go| ylpe)(2 50 — y 3 + ot + de dy. (20) 
b) Die Differentialgleichung fiir die Spannungsfunktion. Zu einer 
direkten Ermittlung der Spannungsverteilung gelangt man durch den Ansatz einer 


Spannungsfunktion D(x, y) 
o@ o®D 


ae han eee hr e) 
wodurch die Gleichgewichtsbedingung (9) identisch erfiillt ist. Dann wird nach (13) 
1 o®@ 1 o® 
Yyz = — Yeh g(to”) “BvD e Vea = GP g(to”) “Oy (22) 
mit 
2 a®@ \* o®@ \* 
om see (2) (8 es 


Setzt man (22) mit (23) in (12) ein, so folgt unter Beachtung von (11) und (15) die 
nichtlineare Differentialgleichung zweiter Ordnung fiir die Spannungsfunktion O(a, y) 


2D ae 4 1 dg(t,2) oi aD 
g(t?) ( On? a SS 3 G2 A(ty?) ( Ox Ox? 


ee ace * 
Die Randbedingung ist fiir einfach zusammenhiingende Querschnitte 
D=0 (25) 
am Rand. Das Torsionsmoment errechnet sich dann zu 
M =\\ (t,e — ty y) dady = 2\ | ®(x, y) dx dy. (26) 


Ill. Der Kreisquerschnitt 


Zuerst behandeln wir kurz den Stab mit kreisfO6rmigem Querschnitt?, fiir den sich 
an Hand der Gl. (16) bis (20) eine genaue Losung einfach angeben laBt. Hierzu setzt 
man 


‘ Y(x, y) = 0, 
das hei®Bt der Querschnitt bleibt eben. Gema8 (16) folgen die Spannungen 


Tz = —Gy(y?)wy, 


(27) 
Ty =Gy(por) w x 
mit é 
yo = 3 oF (a? + 9) 
oder die resultierende Schubspannung 
ary = Gy (5 aw? | wr, (28) 


3 Vgl. auch H. Kauderer, a.a. O., 8. 461. 


138 F. Jindra: 


wo r den Nullpunktabstand bedeutet. Die Randlinie erweist sich nach (19) als Kreis. 
Mit Hilfe von (20) findet man das Torsionsmoment 


- 
; 2 a 
M=22Go | y (5 ot) ar 258 Go| Yi +9 V2 w? RS + 9 408 Be mage ih (29) 
0 
Mgfomt wenn R der Halbmesser der Welle ist. Die Um- 
2000 kehrung der Beziehung (28) wird 


ylr) = or = aoa a) * (30) 


Als Beispiel ist das Ergebnis eines Torsions- 
versuches an Hohlstében aus Flu8stahl in einem 
Schubspannungs-Verdrehungsdiagramm in Abb. 1 
wiedergegeben‘. Der lineare, idealisierte Verlauf der 

Pmt aks Schubspannung ist als Gerade gestrichelt eingetragen. 

0 G25 G50 075 0 1é0%tm’ Nan findet daraus G = 0°87: 10%kg/cm?. Wertet 

Abbe.) Schubspanacnga ver man die Versuchsergebnisse nach (30) aus, so kann 

drehungsdiagramm fiir Flu8stahl innerhalb der Versuchsgenauigkeit die Schubfunk- 

tion g(t?) in guter Naherung durch die ersten 

zwei Glieder der Reihenentwicklung (3) dargestellt werden, also g(t?) = 1+ got, 
mit g, = 0°065- 10°. 


7600 


1200 \— 


G00 


400 [ 


LV. Der elliptische Querschnitt 


Die Auflésung der nichtlinearen Randwertaufgaben (24) mit (25) fir die 
Spannungsfunktion macht bei willkiirlich vorgegebenen Schubfunktionen groBe 
mathematische Schwierigkeiten. Im folgenden wollen wir uns deshalb auf bestimmte 
einfache Formen der Schubfunktion beschranken, was fiir praktische Falle ausreicht. 


Fir den Sonderfall, daB die Funktion g(t,2) in der Form 
gto) =1 4+ gete (31) 
darstellbar ist, vereinfacht sich die Differentialgleichung (24) zu 
Dyy + Py, +e[(3 Dy? + O,?) ®,, + 4B, D, O,, + 


+ (@,2 + 3 @,2) @,,] = — 2@a, (32) 
wobei zur Abkirzung 
_ 292 
ee (33) 


eingefiihrt wird. Die Stoffzahl (33) kennzeichnet hier allein die Abweichung vom: 
Hookeschen Gesetz. 


Zur genaherten Integration von (32) wenden wir die Methode der Stérungsrechnung 
an. Kine notwendige Voraussetzung fiir die Anwendbarkeit des Verfahrens ist, daB 
der gewahlte Stérungsparameter auf hinreichend kleine Werte beschrinkt bleibt, 
um die Konvergenz des Verfahrens zu verbiirgen. Man setzt hier die Lésung des 
Problems als Potenzreihe im St6rungsparameter ¢ an 


D(x, y) = © (x, y) + e OY (x, y) + 2& S@ (a, y) +... (34) 


Indem man (34) in (32) eintragt, ordnet und die Koeffizienten aufeinanderfolgenden 
Potenzen von ¢ Null setzt, entsteht eine Folge von linearen, inhomogenen Differential- 
gleichungen, die man der Reihe nach integrieren kann. 


* Nach M. Ros und A. Eichinger: Versuche zur Klarung der Frage der Bruchgefahr, S. 26. 
Zurich. 1926. 
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Die von « freien Glieder liefern die Differentialgleichung der ungestérten Aufgabe 
| AG — — 2G a, (35) 
2 2 
wo A= a ++ aA den Laplaceschen Operator bedeutet. Dabei ist die nullte 


Funktion ®© (x, y) die bekannte Lésung der linearen Theorie fiir ¢—0. Der 
Koeffizient von « ergibt als Differentialgleichung fiir O(a, y): 


ADY + (3 BY” + BP") DP? + 4BO GO GO + (GL¥ 4 3G") GO — 0, (36) 
Fir die Funktion ®® (x, y) zweiter Ordnung folgt 
AG + (3 BL” + G%) B®, 4 40 G GH + (GO 4 3. G0") OY 
+ 2(3 BY B+ 2G GO 4 GO GO) GY + 
+ 2 (DY B® + 269 6 4 360 GO) GY — 0, (37) 
Fir samtliche Funktionen 6, 6%, D2 usw. gilt die homogene Randbedingung (25). 


Die Spannungskomponenten stellen sich als Potenzreihen in « dar; denn durch 
Differentiation von (34) erhalt man gemif (21) die Beziehungen 


Te a ee Pier Pee a, 
5 (38) 
ty =—3, = es eel a thd ee 
wo 
fot ap) aja. eo” 
en we DY 


usw. die partiellen Ableitungen der Funktionen bedeuten. In gleicher Form 1laiBt 
sich auch das Torsionsmoment (26) schreiben: 
M=MO+¢e¢MYV+2M2+..., (39) 
durch die Festsetzung 
M=2\|GOdady, MY =2\\ O%daxdy usw. 
Ist 


die Gleichung der Randellipse, so wird in der linearen Theorie der Spannungszustand 
im gedrillten Stab mit elliptischem Querschnitt durch die Spannungsfunktion nullter 
Ordnung® 


BD = 6 (a? + 9? y? — a?) (40) 
beschrieben, wobei zur Abkiirzung 
Go a 
C= — nary und Qs i 


gesetzt werden. 
Daraus folgen gemif (21) die Spannungen nullter Ordnung 


tC =20py, t= —2¢x (41) 
und nach (26) das Torsionsmoment nullter Ordnung 
Me — nea, (42) 
Fiir die Niherungsfunktion erster Ordnung erhalt man nach (36) die Gleichung 
AD = — 808 [(3 + g) 2® + (1 + 3Q%) ot y?] (43) 
8 Vgl. 8. Timoshenko, a.a. O., S. 234. 
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mit der Randbedingung 
@®1) — 0 (am Rand). 
Man macht den die Randbedingung befriedigenden Lésungsansatz 
OY” = — 83 (22+ oy? — a?) (A, v2 + A, 0? y? + Ao’). (44) 


Durch Einsetzen findet man die Festwerte 


1 
YA gle ae aie Olen ee 


1 
Ay= 57 (—3 +5 e +19 ott 3 0%), 


er 
a" nm, (1+ 9?) 


3 +16 62 + 1004 + 16 0° + 3 08) 
mit 


n, = 12(1 + 6 6? + 0°). 


Die Verbesserungsglieder erster Ordnung der Spannungen errechnen sich dann 
nach (38) zu 


= euheae [12 (1 + 0%) ga® +(—3 +50? +19 of + 3 08)? y? + | 
1 
Heal - A+ MAF3 ea, 


(3 +19 &§ +5 ot — 38) a + 12 (1 + %) ot ye — 


(1) 8 c? x 
© — 3 
y 3 1+ 607+ 0 


—zewl-AB+AIL+3e) ee 


und fiir die Verbesserung erster Ordnung des Torsionsmomentes ergibt sich nach (39) 


LL 
1 VeWe 


a=a 
z 2 4 
MO = 8 | OY (x, y) da dy = ae eae as. (46) 
Z=0 y=0 


Nunmehr kann die Differentialgleichung fiir die Naherungsfunktion ®@) zweiter 
Ordnung nach (37) angeschrieben werden. Sie lautet: 


32 Ge 
Bia es IY 
(45 + 291 9? + 125 ot — 23 o — 6 o*) x + 108 (1 + 0%)8 gt a? y? + 
+ (— 6 — 23.07 4 125'9% 4 291 @% 45 o?)0* ye — 


1— @? } 
— Tot (3 +e) (8 + 26 G* + 6 ot) pa? a? + 


A@G2) — 


+755 (1+3 0) (6 +2602 +8 o4) ota? y). (47) 


Fir die gesuchte Funktion wahlen wir den am Rand verschwindenden Ansatz 
D2 — 32.05 (a® + 9% y? — a) (Bat + Byg?a®y? + By ot yt + 
+ Ba 2? + B, ea? y? + Bya*). (48) 


Nach einiger Zwischenrechnung findet man 


il 
B, ==, (45 + 921 o? + 4226 04 + 1661 9% — 1024 9% — 114 gl + 45 91), 


B,= — =~ (45 + 651 9? — 1129 94 — 10654 0% — 1129 98 + 651 9" +. 45 12), 
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1 
B, = =~ (45 — 114 9? — 1024 0 + 1661 9° + 4226 9° + 921 0! + 45 912), 


1 
Bs = Fos (45 + 1176 0% + 9483 94 + 27742 0% + 24785 0% + 


+ 22330 9! + 10129 g!2 — 2198 o14 — 1242 916 — 90 018), 
1 


By = Fra gt (— 90 — 1242 o? — 2198 o4 + 10129 go + 22330 0% + 
+ 24785 9! + 27742 ol + 9483 o!4 + 1176 018 + 45 018), 
1] 
Bo = Fre ope (45 + 1086 0? + 8241 gf + 25544 9% + 
+ 34914 08 + 44660 0! + 34914 9l2 + 25544 ol! + 8241 9! + 1086 018 + 45 929) 
mit 


a= 90 (Lice 6.9? 4-6%)(1 +14 9? + 0%). 

Hiermit leiten sich die Verbesserungen zweiter Ordnung der Schubspannungen ab: 
tg = 640 ey (Bs + By) a + 2(B,+ Be vey + 

+ 3 Bs ot y* + (B, + B, — B,) a? x + 2(B, — B) ea? y* + (By — B,) a‘), 
7) — — 642 [3 Boot +2(B,+ Beaty t 

+ (By + Bs) ot y* + 2(B, — B;) a? 2? + (B, + B, — By) ea? y + (By — B,)a*]. 
Mit den gefundenen Werten der Konstanten nimmt die Verbesserung zweiter Ordnung 
des Torsionsmomentes die Form an: 


1 45 + 318 g? + 467 of + 644 0° + 467 9? + 31894 450% 
3 oe (1 + oe)? (1 + 6g? + oF) ORE: 


(49) 


M® = 


(50) 


Als Zahlenbeispiel betrachten wir zihen Stahl von Ziffer 3, fiir welchen man 
6 = 0:0573 - 10-* cm*/kg* errechnet. Es sei: a = 50m, 0 = 2, ow = 0°46: 10% em—. 

Mit diesen Werten ergibt die nach dem dritten Glied abgebrochene Naherungs- 
lésung der Spannungen gemaB (38), (41), (45) und (49) die maximale Schubspannung 
am Endpunkt der kleinen Halbachse ty, = — 1452 kg/cm?. Die lineare Theorie 
liefert hier 70), = — 1600 kg/cm?. Die maximale Schubspannung verringert sich 
demnach gegeniiber dem Wert der linearen Theorie bei gleichem Verdrehungswinkel 
um 9°3%. Wir ermitteln noch das Tor- ty 
sionsmoment. Wir erhalten MO = 
= 785 mkg, wahrend fiir die Naherungs- 
lésung nach (39), (42), (46) und (50) 
M = 738 mkg folgt, also eine Abnahme 
um 6'1% gegentiber dem Wert M©) der 
linearen Theorie. In Abb. 2 ist der Ver- 
lauf der Schubspannungen langs der Halb- 
achsen maSstiblich aufgetragen. Gestri- 
chelt eingezeichnet ist das Ergebnis der 
linearen Theorie. Man erkennt einen Eee: 
Abbau_ der Spannungsspitzen an den Schubspannungen in elliptischem Querschnitt 
Randern infolge der Nichtlinearitat. 

Nun wollen wir uns noch dem Sonderfall zuwenden, da die Reihenentwicklung (3) 
wenigstens nach dem dritten Glied abgebrochen werden kann, also die Schubfunktion 


in der Form 


g(to?) = 1 + Gato" + Gato (51) 
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darstellbar ist. Aus (24) wird dann 


A® + e[(3 G2 + G,2) O,, + 4D, D, B,, + (2 + 3 G,) Oy] + 
+é [(5G,4 + 66207 + By) Oz_ + 8B, Dy (O? + P,?) Pay + 


+ (®,! + 66,207 + 5@,*) B,,] = — 24a, (52) 
wenn folgende Stoffzahlen eingeftihrt werden: 
e=-2% wd ¢=5%. (53) 


Um die geniherte Lésung der Gl. (52) zu erhalten, erweitern wir den Ansatz (34) 
der Stérungsrechnung folgendermaBen: 
P(x, y) = B (x, y) + e DY (x, y) + PD (x,y) +...€ OY (x,y) +... (54) 


Geht man mit (54) in (52) ein und entwickelt nach Potenzen von « und ¢, so folgen 
wieder durch Vergleich gleichhoher Potenzen die Gleichungen ftir die Koeffizienten- 
funktionen. Die Bestimmungsgleichungen fiir ® (x, y), OY (x, y) und OB (zx, y) 
sind bereits in (35) bis (37) angeschrieben. Beriicksichtigt man dabei die Glieder, 


die in erster Potenz enthalten, so entsteht fiir ®® (a, y) die Differentialgleichung 
ABD F(5 OO 6 OO GO Oa) = 
[a8 Dy) OF (Oo eee 
+ (B44 6 GO? GO? 4 5 G4) G _ 0 (55) 
mit der Randbedingung 
@®Y =0 (am Rand). 


Wir betrachten weiter den elliptischen Querschnitt. Fiihren wir (40) in (55) ein, 


so erhalten wir jetzt fiir die zweite Verbesserung erster Ordnung ©® der Spannungs- 
funktion die Differentialgleichung 


AD” = — 32.05 [(5 + 9%) a4 + 6(1 + 0%) otaty?+(1+50%) o%y"], (56) 
welche wir mit Hilfe des die Randbedingung befriedigenden Ansatzes 
OY = — 3205 (a2 + oF y? — a®) (CO, at +O, e8 ay? + 
+ C3 ety? + C,a? 2? + C, ay + C, a4) (57) 


zu erfillen versuchen. Wie man durch Hinsetzen dieses Wertes in (56) leicht nach- 
rechnet, ist 


C;= om Oe + 18 ot — 19 0? + 5 98), 

C= saz (5 +56 0% ~ 282 0! + 56 0° +58), 

C3 = sa (5 — 19." 7718 08S TL Ge oe), 

C= sng (5 + 101 g? + 434 o! + 52 0° + 155 0? — 97 01° — 10 9), 
C,= wane + 155 ot + 52 0% + 434 98 + 101 9! + 5 022), 
Co; = sng (5 + 86 g? + 251 o! — 44 0% + 251 0® + 86 01° + 5 022) 


m, = 30(1 + 14 0? + 04), 
ms = 30 (1 + 1492 +4) (1 + 6 02 +o). 
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Die zugehérigen Verbesserungen der Spannungen nehmen hiermit folgende Werte an: 
t= OY = — 648 gy [(C, + C,) vt + 
+ 2(C, + Cs) P a y? + 30, ot yt + (C1 + Cy — C4) a? 2? + 
+ 2(C, — Cs) Pa y* + (C, — C,) a4], 


Ch) 4 (1) (58) 
Ty = — O,’ = 640° x [3 C; x4 + 2(C, +C,) or a y? + 
+ (C, + Cs) of y4 + 2 (C, — O,) a? a? + 
t (C1 + Ca — Cy) Pa? y? + (Cy — C,) a4. 
Das von diesen Verbesserungen erzeugte Torsionsmoment berechnet sich zu 
pag UG V@n eB 
Mo = 8| | DO (w, y) dw dy = — (5 — 2.02 +5¢')nca’, (59) 
w=0 y=o = 


Wir spezialisieren uns jetzt auf den kreisférmigen Querschnitt, setzen 9 = 1 und 
erhalten nach (40), (44), (48) und (57) 


l 1 : 
DO = —FGo (ety —a), «OY = TH at [(e + yf — al, 


i ee 5 a? [(a + y2)3 — aS], GO) == 68 ow? [(a? + y?)> — aS); 


ferner nach (41), (45), (49) und (58) 


7 ——Go y, t) =Gox, 
7) = G3 w (a? + y?) y, oe = — G w (22 + y?) x, 
1? — — 3G o> (a? + y?)? y, ae = 365 w (2? + y?)? a, 
re = G w (a? + y?)2 y, oad = — G w (22 + y?)? a; 
und schlieBlich nach (42), (46), (50) und (59) 
Mo = S2Goad, MY = —Fa@ora’, 
M® = 3 2G 080°, MO) = —Ta@ oF ab. 


Die Naherungslésungen der Spannungen und des Torsionsmomentes sind also 
Tz, = —Goy [1 — eG? ow? (2? + y?) + (3 & — €) G4 ot (a? + y?)?], 
Ty =Gox [1 — eG? w? (2? + y*) + (3 e? — €) G4 wt (a? + y*)?*], 
M = nGo|5-—3@ ora ae e — 8) otal, 
Beachtet man noch die Stoffzahlen (53) und die Beziehungen (7), so folgt in der 


Tat die nach dem dritten Glied abgebrochene exakte Lésung (27) fiir den Kreis- 
querschnitt 


2 4 
te=—Goy [lL + F720 (8 +9) +5 rot (e+ yy, 
2 = 
oy =Goa|l +3720 (a2 + y?) Hay aw (a + v*)| 
samt dem Torsionsmoment (29) 


as 


2 1 
M=nGo| 5} + > V2 o? a8 + Zrota. 
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V. Der dreieckférmige Quersehnitt 


Als weitere Anwendung der Methode der Stoérungsrechnung behandeln wir den 
in Abb. 3 dargestellten, gleichseitigen dreieckférmigen Querschnitt. Dabei nehmen 
wir den Sonderfall (31) an und beschranken uns bei 
der Durchrechnung auf die zwei ersten Naherungs- 
funktionen. 


Wie beim elliptischen Querschnitt gehen wir auch 
hier von der bekannten strengen Lésung nullter Ord- 
nung aus. Die Spannungsfunktion @, die der Ditfe- 
rentialgleichung (35) geniigt, lautet®: 


BO) = ola® — 3.2 y® —h (a? + ¥) +r M4, (60) 
wenn man zur Abkiirzung 
Go 
/ iginee 


v 


Abb. ‘e es im  gchreibt. Hierin bezeichnet h die Dreieckhéhe. Die 
eicioeeeiogen pore: Spannungen nullter Ordnung haben folgende Werte: 


1 — —BQey(3r+h), = —c(8e2—3y —2h2). (61) 
Das Torsionsmoment nullter Ordnung betragt nach (26) 


1 
5/3 


Oe Go ht. (62) 


Die Naherungsfunktion erster Ordnung soll nun nach (36) die Gleichung 
AGO) = — 403 [27 (a? + y?) (a — 3x y?) — 
— 54h (a? + y?)? + 36 h? (23 — 3 x y*) — 8 h3 (a? 4+ y?)] (63) 


mit der Randbedingung 
@®Y) —0 (am Rand) 


befriedigen. Durch Hinsetzen bestitigt man leicht, daB die Differentialgleichung (63) 
die Losung 


3 4 
oo — 3. aly Bay? —h (a? + y) +3724] |9 (a2 + 42)? — 22h (a® — 3.x y2) + 


32 
+87? (a2 + 42) +3 hl (64) 
besitzt. GemaB (38) ergeben sich die Verbesserungsglieder erster Ordnung der 


Spannungen in der Form: 


3 
1 = Ft y|9a (at + 10.0% y? + 9 y!) — Bh (35 at — 150 a2 y? — 9 yt) + 


+ 60 A? x (a + 3 y2) +P (a? + y2) — Pt a : 
: | 
© — Ty oil9 (7 28 — 5 at y? — 15 a? yt — 3 9) — 6h w (31 24 — 70 a2 y? + 75 yA) + 


80 
+ 8018 (6 a8 — 62% y? — 3 y') — > Ia (a? + y2) — I (a? — y)]. 


; (65) 


6 Vgl. S. Timoshenko, a.a.O., S. 236. 
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Die durch die Spannungen (65) bedingte Verbesserung des Torsionsmomentes berechnet 
sich zu 


ow hs, (66) 


Als Beispiel waihlen wir wieder Stahl von Ziffer 3. Es sei: h = 12cm, wo = 
= 0°306-10%cm-. In Abb. 3 ist die Verteilung der Schubspannung 1, iiber der 
w-Achse mafistiblich aufgetragen. Die Spannungswerte 1 der ipiearen Theorie 
sind gestrichelt eingetragen. Die gré8te Schubspannung tritt in der Seitenmitte auf 


und hat nach (61) und (65) die GréBe thax = — 1389 kg/em?. Die lineare Theorie 
ergibt an dieser Stelle gem&B (61) den Wert 7, = — 1600 kg/cem?. Nach (62) und (66) 


errechnen wir fiir das Torsionsmoment M = 1986 mkg und entsprechend M — 
= 2128mkg. Die maximale Schubspannung verringert sich demnach gegeniiber 
dem Wert, den sie bei gleichem Verdrehungswinkel bei linearem Elastizitatsgesetz 
annehmen wiirde, um 13°2°% und das Torsionsmoment ebenso um 6°7%. 


VI. Der quadratische Querschnitt 


Zum Schlu8B betrachten wir kurz den quadratischen Querschnitt von der Seiten- 
lange 2a. Wir wollen hier nur die Frage nach der Gréf8e der auftretenden gréBten 
Spannung beantworten. Wie in vorigen Abschnitten, beschrinken wir uns dabei 
auf den wichtigen Sonderfall (31) und wenden als Naherungsverfahren zur Lésung 
der Randwertaufgabe (32) mit (25) die Fehlerquadratmethode an. 

Wir machen hier fiir die gesuchte Lésung ® den nur eingliedrigen Naherungs- 
ansatz 

D = ¢ (a? — a®) (y* — a), (67) 


so daB der Ansatz fiir beliebige Werte des Parameters c bereits der Randbedingung (25) 
geniigt. Wir setzen die Niherungslosung (67) in die Differentialgleichung (32) ein und 
erhalten somit die Fehlerfunktion f(z, y;c) in der Form 


f(z, y;¢) = 2Go+2c(?@+y— 207) + 8e[(32ey 4+ 10r7*y +3272 y? — 
— 20 a? x? y? (a? + y?) + at (at + 30 2% y? + y*) — 408 (2? + y?)]. (68) 


Nun verlangen wir, daB das mittlere Fehlerquadrat im ganzen Bereich méglichst 
klein wird. Die notwendige Bedingung fiir einen Kleinstwert ist 


\ | fe, 436) a dxdy = 0. (69) 


Die Integration von (69) tiber die Quadratfliche fiihrt auf die Bestimmungsgleichung 
fir den Parameter c 
78565376 
Fea — Go ee a 2 0? a’ (16c a? — 9G w) + Farge aM = 0. (70) 


Die aus dem Ansatz : 7) herzuleitenden Schubspannungen t, und t, ergeben 
sich nach (21) zu 


te=2ey (2 —¢),  t% = — 20e(y* — a). (71) 
Die in den Seitenmitten auftretende gréBte Schubspannung betrigt dann 


Tone == OC a (72) 
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Zum Vergleich fiihren wir die Naherungsrechnung mit dem Ansatz (67) auch 
fiir die lineare Theorie durch. Mit «= 0 bestimmen wir aus (70) den Wert 


\ Ilmax c* — 15 G B 
hate Mi SOE 
-G10| 2 


und somit die groéBte Schubspannung T.,,, der 
linearen Theorie 


Tax = 13644 wa. (73) 


ma. 


Dieser Naherungswert liegt nur 1:04% tiber dem 
mit Hilfe der Reihenentwicklungen zu_berech- 
nenden Wert der strengen Theorie’ 
as t= 13504 wa. 

a Q7 02 G80 %em7 max 
Abb 2) Abanlline deh pro Bree Selub: Die Anderung der maximalen Schubspannung in- 
spannung fiir quadratischen Quer- folge der nichtlinearen Zusatzglieder in (32) ist dann 

schnitt (74) 


Als Zahlenbeispiel wahlen wir noch einmal Flu8stahl von Ziffer 3 und nehmen 
a—5cman. Um den Einflu8 der Nichtlinearitat fiir diesen Fall besser zu tibersehen, 
sind in Abb. 4 die Anderungen Atpax im Verhaltnis zum Wert Tx der linearen 
Theorie in Abhangigkeit von der spezifischen Verdrehung w dargestellt. 


* 
Atnax = Umax, am Crate 


7 Vgl. S. Timoshenko, a.a.O., S. 248. 
(Hingegangen am 19. Februar 1957) 


Zur rollenden Reibung 
zwischen Scheiben mit verschiedenen Elastizitatskonstanten 


Von K. Desoyer, Wien 
Mit 6 Textabbildungen 


Zusammenfassung. Das Problem der rollenden Reibung zwischen zylindrischen Kérpern mit 
verschiedenen Elastizitaétskonstanten fiithrt unter den Voraussetzungen, daB der Rollvorgang 
stationér und elastisch verlauft und im Gleitgebiet das Coulombsche Reibungsgesetz fiir jedes 
Flachenelement gilt, auf singulare Integralgleichungen zweiter Art fiir die gesuchte Normaldruck- 
und Schubverteilung im Berihrgebiet. Diese lassen sich fiir die Grenzfalle des vollkommenen 
Gleitens und des vollkommenen Haftens geschlossen lésen. Die fiir gleiche Elastizitétskonstanten 
bzw. bei verschiedenen Elastizitatskonstanten nur fiir den reibungsfreien Fall giiltige Hertzsche 
Normaldruckverteilung tritt hier nicht mehr auf, da der Normaldruck durch die Schubspannung 
beeinflu8t wird. Das Gleichungssystem fiir den allgemeinen Fall la8t sich auf eine einzelne 
singulére Integralgleichung fiir eine Hilfsfunktion und Integraldarstellungen fiir die gesuchte 
Normaldruck- und Schubverteilung im Beritihrgebiet zuriickfihren. 


Das. Problem der rollenden Reibung zwischen deformierbaren kreiszylindrischen 
Korpern ist unter den Voraussetzungen, daB der Rollvorgang elastisch verliuft und 
in dem Teil des Beriihrgebietes, in dem Gleiten auftritt, das Coulombsche Reibungs- 
gesetz fiir jedes Flachenelement gilt, fiir gleiche Elastizitatskonstanten der beiden 
Korper geldst. In diesem Fall gliedert sich das Beriihrgebiet in ein Haftgebiet an der 
Anlaufseite und ein anschlieBendes Gleitgebiet. Im ganzen Beriihrgebiet gilt die 
Hertzsche Normaldruckverteilung. Die GréBe des Beriihrgebietes ist durch die 
Materialkonstanten, die Walzenradien und die GréBe der tibertragenen Normalkraft be- 


+ G. Heinrich : Zu L. Féppls Theorie der rollenden Reibung, Osterr. Ingenieur-Archiv IV (1950) 
S. 365—375. Vorangegangene und spitere Arbeiten siehe Literaturverzeichnis im Anhang. 
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stimmt, die Unterteilung des Beriihrgebietes in Haft- und Gleitgebiet durch das Ver- 
haltnis von tibertragener Tangential- zu Normalkraft. 

Im folgenden sollen die Ausgangsgleichungen des Problems, die ,,Hertzsche Be- 
ruhrbedingung“ und die ,,Haftbedingung“‘, auf den Fall verschiedener Elastizitats- 
konstanten der beiden Korper verallgemeinert werden. Fiir Spezialfille lassen sich 
geschlossene Lésungen dieser Gleichungen angeben. Eine Méglichkeit zur Bestimmung 
der Lésung im allgemeinen Fall kann aufgezeigt werden. 


I. Die Beriihrbedingung 


Ks seien zwei deformierbare Kérper von urspriinglich kreiszylindrischer Gestalt 
betrachtet, die mit parallelen Achsen gegeneinander gedriickt und gedreht werden, 
so da zwischen ihnen eine Normalkraft N und eine Tangentialkraft 7’ iibertragen 
wird. Der Vorgang verlaufe stationiir und elastisch. Die Deformationen in dem ge- 
samten Gebiet, in dem eine Kraftiibertragung zwischen den beiden Kérpern statt- 
findet (Beriihrgebiet), miissen von solcher Art sein, daB eine Beriihrung langs einer ge- 
meinsamen Flaiche méglich ist. Das hieraus folgende Kriterium werde Berithrbedingung 
genannt. 

Wenn man voraussetzt, daB die Ausdehnung des Beriihrgebietes in tangentialer 
Richtung sehr klein gegen beide Walzenradien ist, kann man, um die von der Normal- 
druck- und Schubspannungsverteilung im Beriihrgebiet verursachten Deformationen 
zu erhalten, jeden der beiden Korper durch eine elastische Halbebene ersetzen, die 
langs eines Stiickes ihres Randes belastet ist. 

Herrscht lings des Randes einer elastischen Halbebene im Bereich — A Sax <A 
eine Normaldruckverteilung p(x) und bezeichnet man mit 7)(x%) die dadurch ver- 
ursachte Verschiebung des Randes in Richtung der in das Innere der Halbebene 
weisenden y-Achse, so gilt? % 


A 
an”) 2 | pu) du (1) 
Ya aH’ \O zx—y° 

—A 


Fiir eine in Richtung der positiven «-Achse positiv gezihlte, im Bereich — A < 
<a <A verteilte Schubbelastung q(x) ist* 
on q(2) 1 
Se ee (1— =). (2) 


Fir den ebenen Spannungszustand ist in Gl. (1) und (2) H’ = H, fiir den ebenen 


Formanderungszustand H’ = Bae zu setzen, wobei H den Elastizitatsmodul 
und m die Poissonsche Zahl bedeuten. 


Bei Beschrankung auf eine Theorie erster Ordnung kénnen die beiden Anteile der 
Verschiebung itiberlagert werden 

y= 4?) + “®, (3) 

Der Ubersicht halber sei die Herleitung der durch die Verschiebungen ausgedriickten 

Beriihrbedingung® hier nochmals kurz angegeben. Abb. 1 zeigt die beiden Walzen 

oder Scheiben (Radien r, bzw. 7,) vor der Deformation. Mit x, sei der unbekannte 


2 L. Foppl: Drang und Zwang, Bd. III, 8. 62, Munchen 1947. 

8 (c) soll darauf aufmerksam machen, daB das Integral fiir Stellen x im Bereich — A <a@< A 
wegen der singuladren Steile w = x im Sinne des Cauchyschen Hauptwertes zu verstehen ist. 

c ohne Klammer nach dem Integralzeichen bedeutet den Cauchyschen Hauptwert. Die Ab- 
weichung von der iiblichen Bezeichnungsweise ist satztechnisch bedingt. 

4 G. Heinrich, l.c. Gl. (13). 

5 G. Heinrich, l.c. Gl. (16). 
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Abstand der Mitte des sich bei der Deformation ausbildenden Beriihrgebietes von der 
Verbindungslinie der Mittelpunkte M, M, bezeichnet. Die Punkte A, und A, sollen 
sich nach erfolgter Deformation der Korper beriihren. Sie erfahren dabei die elasti- 
schen Verschiebungen 7, bzw. 7, (jede nach dem Inneren des zugehdrigen Korpers 
positiv gezihlt). Bezeichnet man mit A die elastische Anniherung der Kreismittel- 
punkte M, und M,, dann gilt 


A=A,A, +7 + Ne- (4) 
Weiters ist in erster Naherung 
aia 1 ] 1 
Ay A, = > (t+ 2)*(>- + =). (5) 
Aus Gl. (4) und (5) folgt 
0A 1 1 é 
Ba(mt+o(-+5)taeintm)=9 (6) 


fiir alle x im Beriihrungsgebiet — A Sav <A. 

Es sei nun p(a) der auf den Korper 1 im Berihrgebiet 
wirkende Normaldruck, g(x) die auf den Korper 1 im 
Beriihrgebiet wirkende, in Richtung der gemeinsamen 2- 
Achse positiv gezihlte Schubspannung. Dann ist nach Gl. (1) 


A 
Ci” a 2 \ p(w) du 
oe =,Sti(ia HY’ SS ew’ Je 
=i 
A 
ihe, pdt Oe ty 2 p(w) du (8) 
Gio) a hg \o ew 
Abb. 1. Zur Ableitung der und nach Gl. (2) 
Berthrbedingung an, @ q(«) i 
nu — — 4) (1+), (9) 
ang a(x) 1 
a (1—=,). (10) 


Setzt man Gl. (7) bis (10) unter Beachtung von Gl. (3) in (6) ein, dann erhalt die ,,Be- 
rihrbedingung“ die Form 


d 
C g(a) + \ BM = B+ (ay + 2) (11) 
mit den Abkirzungen ‘ 
1 jy 7s 1 rg 
Babint oll Wake ee x (17a) alle. wo) 
os a 7) Hy+H wd C= 5 wa Ey an ae 


Gl. (11) muB fiir alle x im Beriihrungsgebiet — A < x < A erfillt sein. Sie zeigt, 
da8 im Falle verschiedener Elastizitatskonstanten der beiden Kérper die Normal- 
druckverteilung durch die Schubspannungsverteilung beeinfluBt wird, so daB die im 
Falle gleicher Elastizititskonstanten (C = 0) herrschende Hertzsche Normaldruck- 
verteilung hier nicht mehr auftritt. 


II. Die Hafthedingung 


In dem Teil des Beriihrgebietes, in dem Haften eintritt, miissen die Deformationen 
von solcher Art sein, da zwei einander in diesem Gebiet beriihrende Punkte beider 
Korper dieselbe Geschwindigkeit haben. Dies fiihrt zu einem weiteren Kriterium, das 
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, Haftbedingung‘ genannt werde. Zu deren Ableitung diene die Abb. 2. Es seien P, 
und P, zwei Punkte, die im Haftgebiet zusammenfallen. Ihre Geschwindigkeiten w, 
und w, sind daher gleich. Ohne Deformation wiirde P, in P,°, P, in P,° liegen. Die 
elastischen Verschiebungen zufolge der Deformation seien &, bzw. €,. (Wenn der 
Korper 1 der treibende ist, ist £, negativ.) s, und s, seien die zeitlich veriinderlichen 
Kreisbégen D, P,° und D, P,®. M,D, und M, D, seien raumfeste Richtungen, die 
Punkte D, und D, praktisch spannungsfrei. w, und w, seien die Winkelgeschwindig- 
keiten der beiden Kérper. Es gilt dann® 


d d ds ds 
reer ey Wy a (Oy 8S) ew Sg et Ope oa ty Oy 
dé dé dy de Gee des 
Ei Gnt dia dat St dt ak ag, BOR 


dé LE : F re 
Ersetzt man [~ und a durch die Dehnungen ¢, und ¢,, so erhalt man schlieBlich 
1 2 


die Bedingung 
1 @ (1 + &) = 72M (1 + €5). (13) 
Nun sei w., die Winkelgeschwindigkeit, die die Scheibe 2 hatte, wenn sie ohne Defor- 
mation und ohne Schlupf auf der Scheibe 1 abrollen wiirde: 
Te Woo = 1 W3. (14) 
Setzt man 
WO: = Woo — Aa, (15) 
so erhalt man durch Einsetzen von Gl. (14) und (15) in (13) 
unter Vernachlassigung des Produktes ¢, dw 


i) 


Ex(%) — &(%) = 7 1 (16) 


Es mu also die Differenz der Dehnungen in tangen- 
tialer Richtung im Haftgebiet konstant sein. 

Nun werde wie iiblich der Schlupf S, bezogen auf eine 
Umdrehung der Scheibe 1, eingefihrt: 


20 
S = (ry; @; — Tz ) Bie 
Ei 


(17) 


Wegen Gi. (14) und (15) kann man dafiir auch schreiben 


S=207, 2, (18) 
1 
Damit geht die Bedingung (16) tiber in 
S Abb. 2. Zur Ableit der 
Ey (x) — &,(x) = aa ee konst. (19) Haftbedingung ~ 


fiir alle x im Haftgebiet. 

Um die Dehnungen ¢, und ¢, durch die Normaldruck- und Schubverteilung im 
Berithrgebiet auszudriicken, betrachten wir wieder eine Halbebene, auf deren Rand 
im Gebiet — A< x < A eine Normaldruckverteilung p(x) und eine in Richtung der 
positiven xz-Achse positiv gezihlte Schubbelastung q(x) wirkt. Die y-Achse zeige ins 
Innere der Halbebene. Die von p(x) herriihrenden Spannungen lings des Randes 
sind’? 

(a7 pg OP ao =. P(e). (20) 


6 G. Heinrich, l.c. 8. 372. 
7 L. Foppl: Drang und Zwang, Bd. III, S. 52, Miinchen 1947. 
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Eine im Koordinatenursprung angreifende, in Richtung der positiven «-Achse wirkende 
Schubkraft Q (bezogen auf die Lingeneinheit in der z-Richtung) erzeugt in der Halb- 
ebene den Spannungszustand® 


2Q Aig 2Q xy? 


Oo, = — ——— > >? 


e=— at yp? Ou a (att 


Wegen 


1 / 1 
Ey = i’ (o.— =, 0, 


erhalt man fiir den Rand y = 0 


(ey -0 = — GW a: 


Fiir eine an der Stelle x, des Randes angreifende Schubkraft @ ist daher 
Pan ge tah 


(2x)y = 0. = — mE’ %—x,~ 


Daraus erhalt man fiir eine in — A <a <A des Randes verteilte Schubbelastung 
q(x) die Randdehnung 


A 
2 q(u) du 
yes \o wee (21) 


Uberlagerung der von p und q herriihrenden Dehnungen gibt 


A 
(¢.)y=0= — Fp a) alc gee (22) 


Ist wie vorne p(x) der auf den Korper 1 im Berihrgebiet wirkende Normaldruck, 
q(x) die auf den Korper 1 im Beriihrgebiet wirkende, in Richtung der gemeinsamen 
x-Achse positiv gezihlte Schubbelastung, so ist daher 


By 
Ss (23) 
£2) = — ( ny) Pt®) — =e ude | 
—A 
Kinsetzen von Gl. (23) in (19) liefert die ,,Haftbedingung“ in der Form 
A 
C p(x) + | wee = §* = konst., (24) 
—A 


wobei C' die Bedeutung Gl. (12) hat und S* proportional dem Schlupf 8 ist gema8 


= 1 
Up wae ay (25) 


Damit sind die Ausgangsgleichungen des vorliegenden Problems formuliert: die 
Bertihrbedingung (11) muB fiir das ganze Beriihrgebiet — A < x < A ertiillt werden. 
Im Haftgebiet muB die Haftbedingung (24) erfiillt sein. In dem Teil oder den Teilen 


§'L. Foppl, lic. 8. 29. 
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des Beriihrgebietes, wo relatives Gleiten zwischen den beiden Korpern auftritt, soll 
laut Voraussetzung | 


q(x) = fo p(x) (26) 
bzw. 


q(%) = — Mo (2) (27) 
je nach der Richtung des relativen Gleitens, gelten. ju, sei der Gleitreibungskoeffizient. 


III. Bereitstellung einiger Formeln 


Im folgenden wird die Lésung von singularen Integralgleichungen des Typs 
1 
F(@) + 2 Flu) 


loge du = Oz), —lsx<l, (28) 
<1 

benotigt. Wir entnehmen einer Arbeit von H. Séhngen?®: Ist G(x) in (— 1, 1) eine 

Funktion der Klasse LZ? (p > 1) und werden solche Lésungen von Gl. (28) gesucht, 

die in (— 1,1) einer Klasse L¢ (q > 1) angehéren, so gilt (fiir reellen Parameter A, nur 

dieser Fall wird hier bendtigt): die allgemeine Lésung der homogenen Gleichung ist 


Si 1+ 
Fy(2) ar Co n new (Rese - ? (29) 
wobei die Konstante 
1 
C, = | F,(x) dx (29’) 
=] 
beliebig vorgegeben werden kann und 
tgama = —A, —l<a«s0 (30) 


zu wihlen ist. 


Die allgemeine Losung der inhomogenen Gleichung ist 
1 


G(x) A 1 1—wuw f/1l1—u lt+aw\e Gu) 
F(«) = 7 relate 4) pee 

<i 
sin 7% % 1+ a \% ; 

See ate) eer (ey) 

wobei die Konstante 
1 

Cy = | F(x) dx (31’) 


beliebig vorgegeben werden kann. 


Die Gl. (28) bis (31’) sollen nun auf die Grenzen a und b gebracht werden, da dies 
im letzten Abschnitt bendtigt wird. Fiihrt man die Substitution 


20 ee b+a 2 b+a 


Eee SZ Oe 0 aa bee 


® H. Séhngen: Zur Theorie der endlichen Hilbert-Transformation, Math. Zs. Bd. 60 (1934), 


8. 31—5l. 
Gleichungen des allgemeineren Typs 


1 
a(x) F(x) +a\ Et ses — Gz) 


Cc ey—u 
= 
hat erstmalig T. Carleman behandelt: Sur la résolution de certaines équations intégrales, Arkiv 
for Mat. Astr. o. Fysik Nr. 26 (1922), 8.13. Vgl. auch F. Tricomi, Lezionisulle Equazioni Integrali, 
Torino 1954, S. 281—297. 
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in obige Gleichungen ein, nennt F(x) = f(x), G(«) = g(x) und 1aBt die Querstriche 
wieder weg, so erhalt man: Unter obigen Voraussetzungen hat die allgemeine Losung 
der Integralgleichung 


b 
ta) + 4\, AM du = ala), asest (32) 
die Form 
: ; 
g(x) A 1 1 b—a\-* b—u\* (b —u) g(u) 
f(x) = Se he ae REE ~ 5-5 (= =) \. (=) x—U du 
sina a« {b—a\-% x 
Or aes | ee 
mit 
teint 6 = A; —l<a S0. (34) 
Die Konstante 
“ 
O,=\fa)de (35) 


ist beliebig wahlbar. 

Die im folgenden auftretenden Stérfunktionen g(x) sind Konstante oder lineare 
Funktionen von x. In Hinblick auf eine beabsichtigte Fortsetzung des Abschnittes VI 
soll hier bereits das Integral 


—l<a«<0, (36) 


UO «v—U 


b 
VE la = Cote utcw+ cue cn 


a 


berechnet werden, das fiir Gl. (33) benotigt wird, wenn g(w) ein Polynom 2. Grades ist. 
Durch die Substitution 
b—u 
U—a 


erhalt Gl. (36) die Form 


Co (1 + 2)? + ey (1 + 2)? (b + a2) + (1 + 2) (b + a2)? + ¢3 (6 + a2) d 
ne Ze 
(z 4 (2 


Ca a 


Durch Partialbruchzerlegung des Ausdruckes in der Klammer wird 


~ zm dz za dz 
hae \o ha baat TE » | Up area: Say 
0 T—A n=l 0 
mit 
Dy = Dy =.6ai + Cpe + 0..2?.-+ 6,193 
D, = (b — a) [ce, + 6, (@ + a) +, (2? +a 2 + a?)], (39) 
Dz = (b — a) [eg + ¢3 (uw + 2 4)], 
Dy: = (b — aj C3. ; 
Die in Gl. (38) auftretenden Integrale sind bekannt!®. Es ist 
\ pope = — F metgaa (k>0, —1<a«<0), (40) 


0 


10 W. Grobner, N. Hofreiter: Integraltafel II, 421, 13c, 13¢e, 20, Wien: Springer-Verlag. 
(1950). 
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co 
ase ; | 
\eape@= bk =) (hk > 0. 1 xo <0), (41) 
0 
a (O<1+ 2 
2m gia SX oe x < OS Meme be, Bac a.6 
a ae ee ‘ ? TO 1 p) 
eat, ( el rer Dee of al ore ge 23) iA 


x 


ay 0, (—1 SE <4, 0) 


Damit ergibt sich fiir : 


= 2 3 Libs 
= 4 an ae (Cy + C1 @ + Cy x? + ©, 23) |= ae page 


7) 

(2=*)' Cy + Cy U + Cy u? + Cy U8 
G 

a 


0 o 
— (2=2) notgna|— * —{(b—a)afoy + g(x + a) + oy(2® + aa + at) + 


+ (b—ayt [oy +e (e+ 2a)] + (b—ap2=C—FE—2) oh (43) 


ies eo <0): 


b 
*/b—wu\* ey + cu + Cg U® + Cy wi 
=, 2s. = 2 BY nae 
(5) — a (Co + Cy & + Cy X* + Cy X?) a ae 
a 


fi —(F==)|— [6 0) afer + ea (@ + a) + €5 (2% + ae + a2] + 


— 


(0 == oy Y fee + ¢3(" + 2a)] + (b—a)? 


(1 — 2— 
SUSE) oh (44) 


Anwendung von Gl. (33), (34), (35) und (43) gibt: Die allgemeine Lésung aus der 
Klasse L4 (¢ > 1) in (a, b) der Integralgleichung 


b 
f(a) ++ lg HOO dy = dy + dy + de®, a<aSb, (45) 
lautet: x 
fl) = zea (sae) (O—)dot dhe + d,x*) — 
— (b —a) «[(bd, — dp) + (bd, — d,) (x + a) —d, (x? + ax + a*)] — 
a be al dee Pa) 
+ 6— apt FE 4) — 0, (Ge) (46) 
~ tg no = —<A, Sal <io: = .0. (47) 
Dabei ist ; 
Os = | f(x) da (48) 


beliebig waihlbar. Bei der Berechnung von Gl. (46) wurde die aus Gl. (47) folgende 
Identitat 


(sin 7 x)? = T+ (49) 


verwendet. 
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IV. Gleiten im ganzen Beriihrgebiet 


Wenn im ganzen Beriithrgebiet der beiden betrachteten Walzen oder Scheiben 
relatives Gleiten eintritt und vorausgesetzt wird, da8 der Korper 1 der treibende ist, 
dann geht die Beriihrbedingung Gl. (11) wegen Gl. (26) tiber in 


A 
O po p(2) + \pROM =B (ay +2). —ASxSA (50) 


Vil 


Die Haftbedingung (24) gilt in diesem Fall nicht. Gl. (50) ist eine Integralgleichung 
vom Typ (45) mit 


: Tt 
AGED ie: (51) 
und 
Je5ci B 
dy = Git SUE SIRE d, = 0. (52) 


Setzt man — « = 9 und C, = N, so lautet ihre allgzemeine Losung gemaf Gl. (46), 
(47) und (48): 


inz@O 1 A—2x\9 N 
p(a) = BAZ° —* (==) (A— 2) (a+ 240+ x)—2AO(x+ AO) + >| (53) 
mit x 
und is 
A 
N =| p(2) ax, (55) 
A 


N bedeutet die im Beriihrgebiet zwischen den 
beiden Korpern iibertragene Normalkraft, bezogen 
auf die Langeneinheit in der Richtung der 
Walzenachsen. Die Randbedingungen p(A) = 0 
bzw. p(— A) = 0 fordern 


—~2A0(e,+40) +> =0 


70 “OS 0 05 cr ae ae) 
Abb. 3. Normaldruckverteilung im  Daraus folet % +2A@—A=0. 

Falle des Gleitens gt y,=A(1— 26) (56) 

und N= BA22.0 (1 — 6). (57) 


Aus Gl. (57) laBt sich die halbe Breite A des Beriihrgebietes durch N ausdriicken: 


N 
A= 36a—o 8° (58) 


Gl. (53) ergibt schlieBlich mit (54), (56) und (57) 


Heim 7g or(i— 3h (at 


Mit Gl. (59) ist laut Voraussetzung auch g(x) bestimmt: 


q(%) = fo p(2). (60) 


Man iiberzeugt sich leicht durch direkte Integration von Gl. (59), da Gl. (55) und (57) 
tibereinstimmen. In Abb. 3 ist der dimensionslose Ausdruck p* fiir verschiedene 
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Werte von © aufgetragen!. Fiir gleiche Elastizitatskonstanten (Cr==(05.@) a=) 3) 
ergibt sich, wie bekannt?, die Hertzsche Druckverteilung 


p(x) = <= 4 2) 4aVi—(2} mit A=: Pa (61) 


V. Haften im ganzen Beriihrgebiet 


Ks soll nun untersucht werden, ob und unter welchen Bedingungen im ganzen 
Beriihrgebiet Haften eintreten kann. In diesem Fall sind Berithrbedingung (11) und 
Haftbedingung (24) fir — A <x <A zu erfiillen. Durch Addition und Subtraktion 
erhalt man daraus zwei entkoppelte Integralgleichungen des Typs Gl. (45) fiir die 
zwei unbekannten Funktionen p(x) + q(x) und p(x) — q(x): 


A 
A 
p(x) + q(x) + = \c oS du = dyy + dy, x, (62) 
—A 
A 
A Firs 
p(x) — g(a) + = \c ou = de = dey dst (63) 
=A 
mit 
Ay = —j, = = (64) 
und 
St B S* — B B 
dro aa a, dog ae ae dy, pcre dy Ole (65) 
Die Gl. (62) bzw. (63) haben gemaB Gl. (46), (47), (48) die allgemeinen Losungen 
1 A—2\~-% S* 
pa) +42) = — Bogen (Spa) [4-9 + 24am + a) + 
iS N+T 
$24 0, (F + ty —A om) + B | (66) 
bzw. Sy 
i A —x«\— %a[, 
ple) — g(a) = B52 (SSF) "(4 — a) (Ge — ey $24 we — 2) + 
N—T 
+24 o9(S — ay + A avg) — 3 | (67) 
mit 
tg 7% Oo = —A=— Fe (—l<a,<0), tgma,= =| (—1l<«,<90). (68) 
A A 
N=|pla)de baw. T=\ q(x) dz (69) 
=A A 


bedeuten hierin die auf die Lingeneinheit in Richtung der Walzenachsen bezogene 
tibertragene Normal- bzw. Tangentialkraft. 
Wegen Gl. (68) ist 
i Saito = — 1. 
Wir setzen 
een ten iS 01 (70) 


11 Dieselbe Figur ist auch fiir 0 < 0 < 1/, verwendbar, wenn man w durch — x und © durch 
1 — @ ersetzt. 


156 K. Desoyer: 


Aus Gl. (66) und (67) folgt dann: 


p(x) = smn’ (4 — ®(¢ 4+ A)-? (F +4 +249+2)— 
ee yer i taste +2+40)—7F- 
— (A — a) 9 (w+ Ay 4/5 — 9 — 2.4 (18) —a| + 
; PS N—T 
+ (A— a) 9 (2 + A124 (1 — 8) (= — 20 A(1—#)) +5 |} (71) 
und 
g(a) =F AZ? (14 — a) (w+ A 9 ( + 094248 +22) — 
—(A— 2)°-¥(a + Ay? [2.4 0(> Ly +48) we . 
4 (A — 2)? (a + Aye-1 | — ay — 24 (1-8) — | oy 
| —(A— a) * (a + AP} 2.4 (Ee 8) (> —#,—A(1—8)} + “>|. (72) 
Die Randbedingungen p(A) = 0 und q(A) = 0 fordern 
s* . N+T 
24 0( pet A8) t~=0 (73) 
und 
Se N—T 
24 (1—8)|5 eo A(t 8)| + a . (74) 
Die Randbedingungen p(— A) = 0 und q(— A) = 0 fordern weiters 
5 + 2,+2A40—A=0 (75) 
und 
a a, 2A 8) Ae (76) 
Aus den Gl. (73) bis (76) lassen sich x), S*, N und 7 durch A ausdriicken. Man erhilt: 
. Xo as 0, (17) 
S* = A B(1 — 29), (78) 
N=2A?> BS (1 — 9B), (79) | 
== .0, (80) 


Aus Gl. (79) lat sich wieder die halbe Breite A des Beriihrgebietes durch N aus- 
driicken 


4 

A =V rae (81) 
Gl. (78) gibt wegen Gl. (25) den Schlupf S, der in diesem Fall ein reiner Formanderungs- 
schlupf ist, zu 


f 1 ry on b= 29) 1/ er aa 1 
S=47,AB\— +>) 01=— 23) => 
"1 (a; rs By) | a} Vo(—9) Nala re rales + a (82) 
Als Druckverteilung bzw. Schubverteilung erhalt man unter Beachtung von Gl. (70): 
a, AB xa \o wv \l—#* ‘ / a \l—# x \e 
pa) = Se (1-4) +4) til ve) 


= LA Be ie ( +), (83) 
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a8 AB x \o x \l—é wz \l—*o w \o 
Ppeaeacly oa) Mots eal WY eal, 

AES x 

= a ta lade 
= Vato ! (zr) (82) 
Die dimensionslosen Funktionen p** bzw. q** sind in den Abb. 4 bzw. 5 fiir einige 
Werte von % aufgetragen”. Fiir gleiche Elastizitatskonstanten (7 = 0, 2 = 1/,) er- 


I 


Abb. 4. Normaldruckverteilung im Falle des Abb. 5. Schubspannungsverteilung im Berihr- 
Haftens gebiet im Falle des Haftens 


AIG 


Abb. 6. Quotient aus Schubspannung und Normaldruck im Falle des Haftens 


gibt sich wieder die Hertzsche Druckverteilung (61), q(x) = 0 und S* = 0. In Abb. 6 
ist das Verhaltnis q(x)/p(x) fiir verschiedene Werte von # aufgetragen. 


Fir #2 1/, ist lim >= +1 und lim put Ledur 


2—>—A 4—>A 


ax Seite 


Rollen mit Haften im ganzen Beriihrgebiet ist also wegen Gl. (80) nur ohne Uber- 
tragung einer resultierenden Tangentialkraft méglich und auch in diesem Fall fir 
CO +0 wegen Gl. (85) nur dann, wenn der Reibungskoeffizient uw, = 1 ist. 


VI. Vorbereitung der Lisung fiir den allgemeinen Fall der Ubertragung einer Normal- 
und Tangentialkraft 


Im Falle gleicher Elastizitatskonstanten tritt, wenn 0 <|T| <j) N ist, eine Glie- 
derung des Beriihrgebietes in ein vorderes Haftgebiet und ein anschlieBendes Gleit- 
gebiet ein'. Setzt man voraus, daB eine solche ,,Zweigebietslosung‘‘ auch im Falle 


12 Bemerkung" gilt sinngeméB auch fiir #. 
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ungleicher Elastizititskonstanten existiert® und Scheibe 1 die treibende ist, so bilden 
Beriihrbedingung, Haftbedingung und Coulombsches Reibungsgesetz das folgende 


Gleichungssystem : 
A 
C q(x) 4 \ — du = B(a)+ 2) fir —A<a2<A (ganzes Beriihrgebiet), (86) 


=A: 


A 
Gayl (xd —'§+ = konat; fir a= 7 = 4a (Hattpebion (G7) 
ay 
q(x) =u, p(x) fir —ASa<a (Gleitgebiet), (88) 


wobei die Unterteilungsstelle « = a unbekannt ist. Es ist méglich, die beiden Gl. (86) 
und (87) ahnlich wie im vorigen Abschnitt zu entkoppeln und gleichzeitig Gl. (88) zu 
erfiillen: Wir nehmen an, wir hitten die Lésung fiir p(x) und q(x) bereits und setzen 
sie in Gl. (87) ein, jedoch fiir eine Stelle im Gleitgebiet. Dann wird Gl. (87) nicht er- 
fiillt sein, sondern es wird sich auf der rechten Seite ein Ausdruck S* + h(a) ergeben. 
Um den Giiltigkeitsbereich von Gl. (87) auf das ganze Bertihrgebiet zu erweitern, 
addieren wir also rechts eine unbekannte Hilfsfunktion h(x), d. h. wir ersetzen Gl. (87) 
durch 


A 


C p(x) + \e¢ ee du = S*+h(xz) fir —A Sa SA. (87’) 


aay 
h(x) soll die Bedingungen 
h(x) = 0 fiir = SA (89) 
und 


h(— A) = endlich (90) 


erfiillen. Durch Addition und Subtraktion von Gi. (86) und (87’) erhalt man die beiden 
nunmehr fiir — A <w <A geltenden Gleichungen vom Typ (32) bzw. (45) 


A 
i 
p(x) + g(x) 4 Fs \c Ss du = dyy + dy, 2% 4 ce > (91) 
aA 
A 
ye = | 
p(2) — g(x) +22 | PHM) gy = ayy t dye +O, (92) 


res 
wobei /,, 42, dy) bis dy, die Bedeutungen Gl. (64) und (65) haben. Nach Gl. (66), (67) 
und (33) ist, unter Beachtung von Gl. (89), 


1 1 A — pers. | 
| pl) + q(x) = — BGA (23). 
S* 
(A Aes + Xo 2A x +2) +24 oy(F +e —A my) +S] + 
C 1 1 A —2\-% C A —wu)\% (A —u) h(wu) 
oo C? + 7 h(x) C? + 72 A—-2 & + i} \o(2 + a} Sa ae du, (93) 


ae! 


SIN T Oo 


ple) — g(a) = BSG (Sa). 


“S-|+ 


(4 = a) = tot 24 org — 2) 42d org — ay + A oy) — 


13 Siehe Bemerkung am Schlusse dieses Abschnittes. 
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C 1 1 A—2\~% A — wu \% (A —u) h(u) 
+ orp Me) + ee as (SS) \o(2=4) Te Say du (94) 
—A 
mit 
temo, = —4 = —s, = P< <0 tae x, = Ay eee Ay =O (90) 
Wegen Gl. (95) ist «, + «, = — 1. Wir setzen wieder 
== 8. -o,=1—%6, tgnh=F, 0< F< 1, (96) 


Bezeichnet man mit p(x) bzw. q,(x) die Druck- bzw. Schubverteilung im Haftgebiet 
a Sx SA und mit p,(x) bzw. qg(x) diese GrdBen im Gleitgebiet — A < x <a, so 
folgt aus Gl. (93) bzw. (94) mit (96) unter Beachtung von Gl. (89) fiir das Haftgebiet: 


Pul®) + dy(%) = 


a 


1 A—xz\? sin a 1 A —u\-*% (A — wu) h(u) 
Seal eee B It P;(2) ora \(iea) o—U > dl (77) 
bzw 
Py(®) — dy(*) = 
1 A—x2x\1-8 sin a 1 7 / A —u\o-1 (A — wv) h(w) 
=A 
mit 
iy eq lS N+T 
®,(x) = (A — 2)( +a. +240+2)—240(4 L a, Ad} wee: (99) 
und 
S* 
@,(x) =(A x) eet A a} — 
} N—T 
—24(1—0)(4-— 2, — A (1—0)) - 75 (100) 
Wegen Gl. (88) folgt aus Gl. (93) bzw. (94) mit (96) fiir das Gleitgebiet 
(1 + Ho) Pal&) = 
1 A—2x\? sin a 0 1 : Ae \ = 9 (Aa) hn) 
ao (Gea) B P, (x) wea (SS 4} x —U du | + 
—A 
C 
+ Gage M*) (101) 
bzw 
(1 — fo) Pe(X) = 
1 A—~2xz\l—# sin x 0 1 A, — u\e—1 (A wu) iu), 
= 4—; (452) — 8B 1 P,(%) + wae (Spa L—U dw) + 
—A 
C 
+ Gap ge MI). (102) 


Elimination von p,(x) aus Gl. (101) und (102) liefert fiir h(x) die Integralgleichung 


Cugh(x) 4 
1—w (A—« 
oe & +A 


x—U 
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mit 


,—= u+ A\?/A—u\l—? 1 u+ A\l—8/A —u\? 
Kiz0) = ea) (4=5) + he (ee a (4=5) + A104) 
Bemerkenswert ist, daB K(x, x) = 1 ist. 

Damit ist das System GI. (86), (87) auf eine einzelne Integralgleichung fiir h(x) 
zuriickgefiihrt, wobei Gl. (88) bereits erfiillt ist. Nach Lésung von Gl. (103) [mit 
den Bedingungen h(a) = 0 und h(— A) = endlich] ergeben sich py(«), Yq(%) und pg(*) 
durch Einsetzen der Funktion A(x) in Gl. (97), (98) und (101) oder (102). Sollte sich 


dabei herausstellen, daB8 nicht alle Rand- und Ubergangsbedingungen 


pel—A) = 0, py(A)= 9, dy(A)=9, pula) = Pola), da(4) = Mo Pela) (105) 
fiir alle w, und fiir alle in Betracht kommenden 7'/N erfiillbar sind, mi&te an der An- 
laufseite ein zusatzliches Gleitgebiet eingefiihrt werden. Die Beziehungen (85) im Ab- 
schnitt iiber das reine Haften zeigen, da diese Moglichkeit, zumindest fiir kleine (1 
und bestimmte Bereiche 7'/N, in Erwigung zu ziehen ist. Die angegebene Methode 
14Bt sich auch auf diesen Fall verallgemeinern. Gl. (87) ware dann durch zwei Hilts- 
funktionen h,(x) und h,(a) fiir die beiden Gleitgebiete zu erweitern, fiir die sich in 
analoger Weise zwei gekoppelte Integralgleichungen ergeben, die jedoch in punkte- 
fremden Intervallen gelten. Man kann daher, zumindest prinzipiell, eine der beiden 
lésen und diese Lésung in die andere einsetzen und diese losen. Die unbekannten 
Funktionen p(x) und q(x) in den einzelnen Gebieten lassen sich dann in analoger Weise 
durch h,(x) und h,(x) darstellen. 

Beziiglich der ,,Zweigebietsl6sung“‘ sei noch erwahnt, da® sich fiir den Spezialfall 
[to = 1 die Hilfsfunktion h(x) mit Hilfe der im Abschnitt 3 zusammengestellten Formeln 
in geschlossener Form finden léBt. Man erhalt unter Beachtung der angegebenen 
Randbedingungen h(x) = — B(a — x)?(A — x)!—*. Die Darstellungen fiir p(x) und 
q(x) enthalten hypergeometrische Reihen. Da dieser Spezialfall jedoch nicht von 
praktischem Interesse ist, soll er hier nicht gesondert behandelt werden. 
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“ Es sind Methoden bekannt, die es gestatten, Gleichungen vom Typ (103) auf Fredholmsche 
Integralgleichungen zurickzufihren (N.I. Muskhelishvili, Singular Integral Equations, Gro- 
ningen, Holland 1953, S. 323ff.) oder direkt zu lésen (W. D. Kupradse, Randwertaufgaben der 
Schwingungstheorie und Integralgleichungen, Berlin, 1956, 8S. 189ff.). Wenn die Lésung von 
Gl. (103) mit ertraéglichem Rechenaufwand médglich ist, soll diese in einer spateren Arbeit folgen. 
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Otto Lange — 70 Jahre 


Der Jubilar, der am 26. April in voller Frische und unermiidlicher Tatigkeit seinen 
70. Geburtstag feiern durfte, kann wohl mit Recht auf ein auBerordentlich erfolgreiches 
Lebenswerk zuriickblicken. Ihm ist es zu danken, daB ein so reprisentatives Sprachrohr 
wie der in der ganzen Welt bestens bekannte und geschiitzte Wiener Springer-Verlag 
den 6sterreichischen Wissenschaftlern und Technikern fiir ihre Publikationen zur 
Verfiigung steht. 


Otto Lange hat ein Leben lang seine groBe Begabung und seine nie erlahmende 
Tatkraft dem Verlagswesen gewidmet, insbesondere dem Springer-Verlag, in dessen 
Berliner Stammhaus 1910 seine Laufbahn begann. Schon 1924 wurde er Leiter des 
Wiener Springer-Verlages, dessen geschaftsfiihrender Gesellschafter er seit 1935 ist. 
Als dann im Jahre 1945 das Wiener Geschaftshaus einem Brand zum Opfer fiel und 
der Erfolg jahrelanger Arbeit vernichtet war, ging Otto Lange mit gréBtem FleiB und 
dem ihm eigenen bewundernswerten Optimismus unter schwierigsten Verhaltnissen an 
den Wiederaufbau des Wiener Verlages, und die groBe Zahl bedeutender wissenschaft- 
licher Werke, die er in rund zehn Jahren auflegen konnte, beweist seine Tatkraft. 
Schon Ende 1945 erschienen eine ganze Reihe 6sterreichischer Zeitschriften wieder 
und neue wurden gegriindet, darunter auch das ,,Osterreichische Ingenieur-Archiv‘“‘, 
dem er besondere Sorgfalt angedeihen 1i8t. Damit ist ein reprasentatives wissenschaft- 
liches Organ fiir die 'Techniker geschaffen worden, das nicht nur | ereits anerkannten 
Fachleuten, sondern auch jungen Habilitanten und Dissertanten die Verdffentlichung 
langerer Arbeiten erméglicht. 


Otto Lange, dem das Ansehen des Wiener Springer-Verlages besonders am Herzen 
liegt, hat schon viele an sich ausgezeichnete Werke verlegt, auch wenn sie von vorn- 
herein keinen grofen Verkaufserfolg versprachen. In Anerkennung seiner Verdienste 
um die technischen Wissenschaften verlieh ihm die Technische Hochschule Wien im 
Jahre 1947 und die Montanistische Hochschule in Leoben 1956 die Wiirde eines 
akademischen Ehrenbirgers. 


Die ésterreichische Technikerschaft wiinscht Herrn Lange, dessen Frische, Leb- 
haftigkeit und Optimismus sein Alter kaum erraten lassen, auch fiir die Zukunft allen 
Erfolg bei bester Gesundheit. F. Magyar, Wien 


Mitteilung 


Jahrestreffen 1957 der Verfahrens-Ingenieure 


Das diesjahrige Jahrestreffen der Verfahrens-Ingenieure wird vom 6. bis 9. Oktober in 
Stuttgart stattfinden. Es wird veranstaltet von der VDI-Fachgruppe Verfahrenstechnik, der 
Fachabteilung Apparate und Anlagen fiir die chemische und verwandte Industrien in der Fach- 
gemeinschaft Apparatebau des VDMA und der Forschungsgesellschaft Verfahrenstechnik. 


Nach einer Plenarveranstaltung am Vormittag des 7. Oktober, auf der unter anderem Prof. 
Dr. J. Bartels, Géttingen, iitber das Internationale Geophysikalische Jahr berichten wird, 
finden am Nachmittag des 7. und am 8. Oktober die Fachvortrage in zwei parallel laufenden 
Veranstaltungen statt. Fiir den 9. Oktober sind Besichtigungen und Exkursionen vorgesehen. 


In zeitlichem Anschlu8 an das Jahrestreffen, vom. 10. bis 12. Oktober, wird die Jahrestagung 
der Fédération Européenne des Constructeurs d’Equipement Pétrolier (FECEP) abgehalten. 
Hierbei werden vor allem Fragen der korrosionsfesten Baustoffe und der MeB- und Regelgerate 
behandelt werden. 

Auskiinfte iiber das Jahrestreffen erteilt die VDI-Fachgruppe Verfahrenstechnik, Frank- 


furt/M., Rheingau Allee 25. Niahere Einzelheiten iiber die FECEP-Tagung sind bei der Fach- 
gemeinschaft Apparatebau des VDMA, Diisseldorf-Oberkassel, Lueg-Allee 63, zu erfahren. 
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Vibration Analysis Tables. Von R. E. D. Bishop und D.C.Johnson. VIII, 59 S. Cambridge: 
At the University Press. 1956. 10s. 6d. 

Die Durchrechnung von zusammengesetzten linearen Schwingungssystemen kann, obwohl 
sie keine grundsatzlichen mathematischen Schwierigkeiten bietet, recht zeitraubend sein. Man 
wird es daher sehr begriiBen, da’ im vorliegenden Heft eine Zusammenstellung fertiger Formeln 
gegeben wird. Im einzelnen enthalten die Tabellen HinfluBzahlen fur Systeme von einem bis 
drei Freiheitsgraden und fiir aus Untersystemen zusammengesetzte Systeme. AnschlieBend wird 
eine Tabelle von Frequenzgleichungen gebracht. Dann folgen EinfluBzahlen, Eigenwerte und 
Eigenfunktionen fiir die schwingende Saite, fiir die Drehschwingungen von Stében mit Kreis- 
querschnitt und fiir die Langs- und Biegeschwingungen von Balken mit konstantem Querschnitt. 
Querkrafteinflu8 und Rotationstragheit werden nicht beriicksichtigt. Eine Reihe von Hilfs- 
funktionen, ebenso wie die ersten fiinf Eigenfunktionen des querschwingenden Balkens unter 
verschiedenen Lagerungsbedingungen sind tabuliert. 

Darstellung und Druck sind sehr klar und iibersichtlich. Druckfehler konnte der Besprecher 
nicht entdecken, obwohl er eine Reihe von Stichproben machte und auch einige kompliziertere 
Formeln nachprifen lief. 

Fiir den Benutzer, der die Tafeln unabhingig von dem von den Verfassern geplanten Lehrbuch 
verwenden méchte, ware vielleicht die Wiedergabe eines oder mehrerer Anwendungsbeispiele 
wiunschenswert gewesen. 

Das Buch kann jedem, der mit der praktischen Bearbeitung von Schwingungsaufgaben zu 
tun hat, durchaus empfohlen werden. H. Parkus, Wien 


Technische Strémungslehre. Von B. Hck. Fiinfte Auflage. Mit 407 Textabb., X, 422 8. Berlin- 
Géttingen-Heidelberg: Springer-Verlag. 1957. Geb. DM 29.40. 


Ks ist erfreulich zu lesen, daB bereits nach knapp drei Jahren eine neue Auflage von oben- 
stehendem Werk notwendig geworden ist, weil man daraus ersehen kann, da das Interesse an 
der Strémungslehre gestiegen ist. Anderseits ist es aber auch ein Zeichen, dai die Technische 
Stromungslehre von B. Eck, von der bereits die fiinfte Auflage nunmehr erscheint, ein sehr 
gutes Buch fiir Studium und Praxis ist, das bei jeder neuen Auflage verschiedenen Wiinschen 
durch Verbesserungen und neue Abbildungen entgegenkommt. 

Es ist zu wtinschen, dal dieses Werk auch weiterhin weiteste Verbreitung finden médge. 

R. Bruniak, Wien 


Untersuchungen iiber die dreidimensionale Potentialstré6mung durch axiale Schaufelgitter. Von 
Th. Ginsburg. (Mitteilungen aus dem Institut fiir Aerodynamik an der ETH in Ziirich: Nr. 22.) 
Mit 36 Textabb., 79S. Ziirich: Verlag Leemann. 1956. 


Das Problem der zweidimensionalen, inkompressiblen Potentialstr6mung durch ebene Schaufel- 
gitter, wie es bei der Abwicklung von koaxialen Zylinderschnitten der Axialrader auftritt, 
ist nach der Singularitétenmethode oder mit Hilfe der konformen Abbildung in grundlegenden 
Arbeiten behandelt worden. In der vorliegenden Arbeit hat sich der Verfasser die Aufgabe ge- 
stellt, die von der Dreidimensionalitat der Stroémung herrithrenden Zusatzeffekte zu ermitteln. 
Er denkt sich die Schaufeln als Schraubenflachen mit geraden radialen Erzeugenden, die als 
Stabwirbel aufgefaBt werden. Das Geschwindigkeitsfeld des Wirbels wird durch die zylindrische 
Begrenzung geandert und die sich dadurch ergebenden Zusatzgeschwindigkeiten folgen aus einer 
langeren Rechnung. Die gréBten Zusatzeffekte zeigen sich am Naben- und Au®enzylinder. 

Der Verfasser verbindet mit seiner ausfiihrlichen theoretischen Arbeit auch experimentelle 
Untersuchungen auf. Grund der elektrolytischen Analogie. Wenn auch die Korrekturen nicht 
groB sind, so ist die Berechnung derselben doch sehr interessant. R. Bruniak, Wien 


Engineering Structural Failures. The Causes and Results of Failure in Modern Structures of 


Various Types. Von R. Hammond. Mit 125 Textabb., 2248. London: Odhams Press Ltd. 
1956. 25s. 


Erfahrungen kénnen nur langsam und auf Kosten von Irrtiimern und Fehlern gesammelt 
werden. Deshalb ist es gerade fiir den Ingenieur sehr wichtig, die Ursachen und Folgen von 
Fehlkonstruktionen kennenzulernen. Der Autor unternimmt es in dankenswerter Weise, ver- 
schiedene der wichtigsten im vergangenen Jahrhundert aufgetretenen Bauschaden und Kate 
strophen zu untersuchen, die vielleicht bereits in Vergessenheit geraten sind oder als nicht restlos 
geklart galten. 

Die angefiihrten Beispiele sind aus nahezu allen Fachgebieten der Ingenieurpraxis aus- 
gewahlt und betreffen: Schittungen, Dimme, Hafenanlagen, Bauten, Briicken, Schwingungs- 
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probleme, Tunnels und SchweiSkonstruktionen. Als AbschluB zieht der Autor die Lehren aus 
den Ursachen der Versager. 

Besonders interessant ist die Untersuchung der Ursachen des Bruches der Kabinenwand 
des ,,de Havilland‘‘-Diisenflugzeuges Comet G-ALYP, der zu schweren Katastrophen geftthrt 
hat und noch allgemein in Erinnerung sein diirfte. Der Verfasser bringt die von der ,,de Havilland‘ 
in eigenen ,,Aircraft Testing Tanks“ durchgefiihrten Ermiidungsteste, durch welche die Ursache 
der Zerstérung schlieBlich festgestellt werden konnte. FF, Magyar, Wien 


Caleulo de Concreto Armado. Von 7. van Langendonck. Band I. Zweite Auflage. Mit 234 Text- 
abb. und 1 Tabellenanhang, XXXI, 5008. Sao Paulo: Associacao Brasiliera de Cimento 
Portland. 1954. 


Der erste Band von Langendoncks Werk ,,Die Berechnung des Stahlbetons‘‘ mit 500 Seiten 
Umfang wird durch einen zweiten Band mit 634 Seiten ergiinzt. Eine nicht gewohnlich breite 
Darstellung der Probleme des Stahlbetonbaus. Es ist aber nicht alles Gebrachte Lehre vom Stahl- 
beton. Dieser 1. Band bringt im 1. Kapitel auf 139 Seiten die Grundlagen der Statik, Bau- 
mechanik und Festigkeitslehre mit deutlicher Hinordnung auf die Erfordernisse des Stahlbetons, 
wobei Kraftspiel und Verformung gleichermafen Beriicksichtigung finden. Erganzt werden 
diese theoretischen Darlegungen durch ausfiihrliche Tafeln (Seite 403 bis 481). Diese vereinfachen 
die Rechenarbeiten bei der Ermittlung von Flachenmomenten und der aus diesen abgeleiteten 
statischen GréBen und dienen zum anderen Teil als Hilfsmittel fiir die Berechnung statisch un- 
bestimmter Balken und einfacher Steifrahmen. 

Die beiden Kapitel Il und III mit zusammen 258 Seiten Umfang befassen sich mit dem 
eigentlichen Thema, dem Stahlbeton. 

Das erstere ist der Beschreibung der Materialeigenschaften des Betons und des Bewehrungs- 
stahles gewidmet. Diese Beschreibung erfolgt sehr ausfiihrlich, wird durch zahlreiche, verglei- 
chende Diagramme anschaulich unterstiitzt und in zahlreichen Fu8noten wird der Leser auf das 
Schrifttum verwiesen, so da dieser reiche Literaturhinweis dieses Kapitels II fast zu einem 
Nachschlagwerk tiber alle Fragen der Beton- und Bewehrungsstahltechnologie macht. Selbst- 
verstandlich wird auch auf die in verschiedenen Kulturstaaten bestehenden Vorschriften tiber 
die im Berechnungsverfahren anzuwendenden Materialkonstanten ausfiihrlich und vergleichend 
eingegangen. 

Das III. Kapitel behandelt den Stahlbetonstab mit Druck und Zugbeanspruchung. Der 
Verfasser geht von der allgemeinsten Problemstellung aus. Er behandelt also zuerst die Aufgaben 
des mit Biege- und Schubspannungen kombinierten Lastfalles des Druck- bzw. Zugstabes, um 
dann zu den einfachen Aufgaben des reinen Druckes bzw. Zuges zuriickzukommen. Higen- 
spannungen des Stahlbetonstabes, die aus der Raumunbestindigkeit des Betons (Schwinden, 
Kriechen) folgen, sowie die Probleme des vorgespannten Stabes finden ihre Behandlung. Die 
schriftliche und rechnerische Darstellung der Probleme findet in zahlreichen Literaturhinweisen, 
guter Bebilderung und geniigender Ausstattung mit Schaubildern vergleichender Art ihre sinn- 
volle Erganzung. 

Sonderprobleme reiner Druckbeanspruchung, wie die Verankerung von Spanngliedern in 
Spannbetonbalken, Balkenlagerkérper, Fundamentgestaltung und Widerlagerkonstruktionen, 
werden schon hier behandelt. Die Behandlung dieser Sonderprobleme zwingt zu eimem weiten 
Ausholen auf Gebiete, die mit dem eigentlichen Druckstab nicht direkt zusammenhangen. Es 
finden sich z. B. zahlreiche gut bebilderte Darstellungen iiber die Spannstahlfithrung in Spann- 
betonbalken. Diese Wahl der Stoffgruppierung mag dazu beitragen, da dieses Kapitel IIIT einen 
gréBeren Umfang hat, als er uns gewohnlich in Stahlbetonbiichern entgegentritt. 

Die gute Ausstattung und theoretisch einwandfrei fundierten Berechnungsmethoden, auch fir 
nicht alltaégliche Lastfalle, sowie zahlreiche Rechenbeispiele machen dieses Buch allen Stahl- 
betoningenieuren besonders wertvoll. F. Pongratz, Wien 


Probleme der Plastizitatstheorie. Von W. Prager. (Lehr- und Handbiicher der Ingenieurwissen- 
schaften: Band 17.) Mit 52 Textabb., 100 S. Basel und Stuttgart: Birkhauser Verlag. 1955. 
Geb. sfr. 12.50. 

Einer kurzen Einfiihrung in das mechanische Verhalten plastischer Stoffe folgen einige An- 
wendungsbeispiele, so die Behandlung eines I-Tragers unter Biegung und Liangskraft, eines diinn- 
wandigen Rohres unter Zug und Torsion, ferner statisch unbestimmte Fachwerke. Den breitesten 
Raum nimmt das Traglastverfahren ein mit seinen Anwendungen auf Rahmen, Platten und 
zylindrische Schalen. Zum Abschlu8 folgen Betrachtungen uber endliche plastische Form- 
anderungen, wie sie z. B. beim Stangpressen vorkommen. 

Das Buch verdankt seine Entstehung einer Reihe von Vortragen, die der Verfasser an der 
Hidgendéssischen Technischen Hochschule in Zirich Ende 1954 hielt. Die Darstellung ist aufSer- 
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ordentlich fliissig, und manches alte Problem der Plastizitatstheorie wird in ganz neuer Beleuchtung 
behandelt. An Stelle mathematischer Ableitungen wird oft und mit groBem Geschick mit 
mechanischen Gedankenmodellen und geometrischen Darstellungen gearbeitet. 

Das Werk kann als eine willkommene Erginzung zu Prager-Hodges ,,Theorie ideal plasti- 
scher Kérper‘‘ (deutsch von F. Chmelka; Springer-Wien, 1954) betrachtet werden, kann aber 
auch unabhingig davon mit Genu8 gelesen werden, sofern man nur einigermaBen mit den Grund- 
lagen der Plastizitatstheorie vertraut ist. F.. Chmelka, Wien 


Finfiihrung in die Technische Mechanik. Nach Vorlesungen. Von J. Szabé. Z weite, verbesserte 
und erweiterte Auflage. Mit 492 Textabb., XII, 397 8S. Berlin-Géttingen-Heidelberg: Springer- 
Verlag. 1956. Geb. DM 22.50. 


Die Notwendigkeit einer Neuauflage dieses elementaren Lehrbuches innerhalb von zwei 
Jahren zeigt, wie stark die Nachfrage der Studierenden nach einer modern geschriebenen 
technischen Mechanik ist. Das ausgezeichnete Werk von Ziegler werden sich ja nur wenige an- 
schaffen kénnen. 

Der Aufbau des Buches ist gegentiber der ersten Auflage unverindert geblieben, doch wurde eine 
Reihe von Besprecherwiinschen beriicksichtigt und die erfreulich groBe Anzahl von Ubungs- 
aufgaben weiter vermehrt. Einige Flichtigkeiten sind allerdings noch geblieben, etwa bei der 
Ableitung des Schwerpunktsatzes, wo unnétigerweise vom starren Kérper ausgegangen, als An- 
wendung aber Gescho8 und Geschiitz bzw. Boot und Insassen gebracht wird. Auch iiber den 
Energiesatz wire wohl mehr und Scharferes zu sagen. 

Sehr schén sind die zahlreich eingestreuten historischen Bemerkungen. 

Man wird das Buch dem Ingenieurstudenten gerne empfehlen. H. Parkus, Wien 


Engineering Inspection, Measurement and Testing. Von H. C. Town und R. Colebourne, Mit 
176 Textabb., 192 S. London: Odhams Press. Ltd. 1956. 21s. 


Diese umfassende Darstellung des Priifens und Feinmessens zur zahlenmaBigen Bestim- 
mung von Lingen, Winkeln, Formen, Geradheit, Ebenheit und Rauheit von 
Oberflachen technischer Kérper in Werkstatten und Kontrollabteilungen neuzeitlicher Fertigungs- 
betriebe verliert sich nicht in bloBer Beschreibung von MeBgeraten und MeSmaschinen, sondern 
behandelt in eigenen Abschnitten auch die Grundlagen mechanischer, optischer, elektrischer, 
pneumatischer und elektronischer Me®verfahren. Das Buch zeigt an praktischen Beispielen, 
wie mit einfachen mechanischen MeBmitteln MeBaufgaben gelést werden, die erfahrungsgemaB 
den Studierenden Schwierigkeit bereiten im Herausfinden des oft recht einfachen mathematisch- 
trigonometrischen Problems und semer Formulierung. Hiebei auch den kleinsten Rechnungs- 
schritt anzugeben, fuhrt allerdings zu unndtiger Breite. 

Anderseits vermiBt man selbst den bescheidensten Hinweis fiir das Messen an Zahnradern 
und fast jeden Schrifttumsnachweis. Auf die unbedingte Notwendigkeit einer geordneten Mef- 
abteilung fiir die Reihenfertigung und zur Durchfithrung von Qualitétskontrollen wird eingangs 
hingewiesen. 

Den Verfassern ist es gelungen — wohl aus ihrer praktischen Tatigkeit am City and Guilds 
of London Institute bzw. am Keighley Technical College heraus —, fiir die Studierenden dieses 
Schultyps ein sehr gutes Lehr- und Vorbereitungsbuch fiir die Priifungen in dem neu in den 
Studienplan aufgenommenen Gegenstand ,,MeBtechnik‘’ zu schaffen. 

Fur den deutschsprechenden MeBtechniker liegt das ,,Taschenbuch fiir LangenmeBtechnik‘‘ 
und die Schriftenreihe ,,Technisches Messen in Hinzeldarstellungen“‘ vor. Dennoch ist das Buch 
von Town-Colebourne fir jenen von sehr beachtlichem Wert, weil es verliBliche Angaben 
uber das derzeitige englische Toleranz- und Passungssystem, iiber die gegenwirtig wichtigen 
Gewindeformen und ihre Toleranzen, Ausfithrungen iiber Mehrfach-Priifeinrichtungen, nicht so 
leicht zugadngliche Nachrichten uber ,,Messen wahrend des Arbeitens‘‘ und weil es den neuesten 
Entwicklungsstand von MeBgeraten und Maschinen spezifisch-englischer Herkunft enthalt. 

Rk. Vockenberger, Wien 


Higentiimer, Herausgeber und Verleger: Springer-Verlag in Wien I, Mélkerbastei 5. — Fiir den Inhalt verantwortlich: Prof. 
Dr. Franz Magyar, Wien IV, Technische Hochschule, Karlsplatz 13. — Druck: Manzsche Buchdruckerei, Wien IX, Lust- 
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Soeben erschien : 
Die Mikrophotographie 
Von 


Dr. Kurt Michel 


Aalen/Wiirttemberg 


Mit 550 teils farbigen Textabbildungen. XXXII, 740 Seiten. Gr.-8°. 1957 
(Die wissenschaftliche und angewandte Photographie, Band X) 


Ganzleinen S 840.—, DM 140.—, sfr. 143.40, $33.30 


Das Mikroskop ist schon lange nicht mehr nur das Handwerkszeug einer kleinen Gruppe von 
Forschern. Es ist heute in praktisch jedem Laboratorium ein unentbehrliches Hilfsmittel und 
in sehr vielen Fallen spielt es auch bereits in der industriellen Fertigung eine wichtige Rolle. 
Zum Festhalten der Beobachtungsergebnisse bedient man sich in der Regel der photographischen 
Verfahren. Die Kombination von Mikroskopie und Photographie, die Mikrophotographie, hat 
sich zu einer wichtigen Technik entwickelt, und es wurde schon lange als schwerwiegender 
Mangel empfunden, daB seit beinahe 20 Jahren eine zusammenfassende systematische Darstel- 
lung dieses Gebietes in deutscher Sprache fehlt. Im vorliegenden Band unternimmt es ein seit 
langem bekannter Fachmann, der als Leiter der Abteilung fiir Mikroskopie der Firma Car ZEIss 
mit den Geraten besonders vertraut ist, eine solche, den neuesten Stand beriicksichtigende 
Darstellung zu geben. Der Band gliedert sich in drei Teile. Im ersten wird ein konzentrierter 
Uberblick iiber die theoretischen Grundlagen der Mikroskopie und Mikrophotographie gegeben, 
ohne deren Kenntnis eine Beherrschung der Praxis immer Stickwerk bleiben mu8. Der zweite 
Teil behandelt die Gerate und der dritte schlieBlich die mikrophotographischen Methoden. Die 
Ausfihrungen werden durch eine groBe Zahl zum iiberwiegenden Teil fiir das Buch eigens neu 
angefertigter Abbildungen erlautert. 


Rentabilitat. Fehlinvestitionen, ihre Ursache und ihre Verhiitung. Von Oberbau- 
rat Ing. B. M. Gerbel, Wien. Zugleich zweite, umgearbeitete und erginzte Auflage 
der ,,Rentabilitat industrieller Anschaffungen‘‘. Mit 26 Textabbildungen. IX, 264 Seiten. 
Gr.-8°. 1955. Ganzleinen § 174.—, DM 29.—, sfr. 29.70, $6.90 


» Das tiberaus anregend geschriebene Buch von Gerbel erfiillt zwei Aufgaben; es enthalt einmal 
wichtige Hinweise fiir das Vermeiden von Fehlinvestitionen, Ratschlige, die in einer Zeit kon- 
junkturellen Aufschwungs ganz besonderen Wert haben, zumal sie von einem Wirtschaftsberater 
kommen, der sie sich in fiimfzigjéhriger Praxis erarbeitet hat; zum anderen vermittelt ein 
Ingenieur seinen Kollegen — auch den angehenden — allgemeine Wirtschaftskenntnisse 
und Verfahrensgrundsatze wirtschaftlichen Handelns... — Im ersten Kapitel werden die 
grundlegenden Unterschiede zwischen der technischen und kaufmannischen Denkweise auf- 
gezeigt; es folgt das Hauptkapitel ,Uber Rentabilitaét‘. Daran schlieBen sich die Ausfihrungen 
,Uber das Kaufen im allgemeinen und das Investieren im besonderen‘, tiber ,Fehlinvestitionen 
und ihre Ursachen‘ und schlieBlich iiber ,Heilmittel gegen die Krankheit des Wirtschafts- 
objektes*.* VDI-Zettschrift 


Die Rentabilitat der Produktivitatssteigerung. kin Beitrag 
zur Lehre der Rechenhaftigkeit der Wirtschaft. Vor Oberbaurat Dipl.-Ing. B. M. Gerbel, 
Wien. Mit 4 Textabbildungen und 1 Ausschlagtabelle. V, 48 Seiten. Gr.-8°. 1952. 

Steif geheftet S 34.—, DM 5.60, sfr. 5.80, $1.35 
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F. HECHT-Wien 
Die Zeitschrift erscheint vierteljahrlich (1957: 3. Jahrgang) 


Jahresbezugspreis 8 216.—, DM 36.—, sfr. 37.20, $ 8.60 


Ab 1957 erscheinen die ,,Astronautica Acta‘‘ bei unverandertem Bezugspreis 
im betrichtlich erweiterten Umfang von 256 Seiten pro Jahr. Trotz bedeuten- 
der Kostensteigerungen in den letzten Monaten hat sich der Verlag zu dieser 
_ MaBnahme entschlossen, um dem wachsenden Bediirfnis nach Veréffentlichung 


der immer zahlreicher einlaufenden wichtigen Beitrige entsprechen zu k6nnen. 


Mitglieder der Internationalen Astronautischen Féderation (IAF) 
bzw. deren Mitgliedsgesellschaften kénnen die Zeitschrift als das 
Organ der IAF zu einem um 20°/, ermaBigten Preis beziehen 


Band III, Heft 1 (1957): Canney, H. E., Jr., and F. I. Ordway, Ill, The 
Uses of Artificial Satellite Vehicles. Part II. - Klemperer, W. B., and R. M. 
L. Baker, Jr., Satellite Librations. — Kaeppeler, H. J., und G. Baumann, Uber 
Systeme mit chemisch reagierenden Komponenten im Gleichgewicht. — Hecht, 
F., und R. Patzak, Chemische Analyse von in Tiefseesedimenten gefundenen 
mikroskopischen Kigelchen vermutlich kosmischen Ursprunges. - Cremona, 
C. E., Metodo fotografico per la determinazione della resistenza aerodinamica 
di missili in volo. - Ordway, F. 1., TI, Project Vanguard - Earth Satellite 
Vehicle Program. Characteristics, Testing, Guidance, Control, and Tracking. 
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